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WAYS OF SAVING STEEL IN REACTOR CONSTRUCTION 


A. N. Komarovskii 


Translated from Atomnaya fnergiya, Vol. 7, No. 3, pp. 205-214 


September, 1959 
Original article submitted April 4, 1959 


The basic methods for saving steel in reactor construction are considered in this article. It is shown that, in sta- 
tionary nuclear reactors, the basic method for saving steel is its replacement with concrete. The article provides 
data indicating that the mechanical strength and the radiation resistance of concretes are not being utilized to 
the fullest extent at the present time, and that capillary cracks, which appear in concrete, do not affect its pro- 
tective properties. The problem of the application of protective explosion-proof shells is considered. The proj- 
ect of the Soviet Boiling Reactor is described as an example proving the efficiency of using a protective chamber 
(instead of a protective shell) which is designed for excess internal pressure. 


In reactor construction, a considerable saving of 
steel can be achieved by a rational and efficient 
design of the entire reactor plant. 

Let us consider a few basic design approaches 
whereby such a saving could be secured: 

1) the use of prestressed reinforced concrete 
instead of steel in the construction of the reactor 
vessel; 

2) dispensing with thermal shielding (usually 
steel); 

8) the use of the reactor chamber, which consists 
of reinforced -concrete side and top shielding, as an 
air-tight explosion-proof shielding instead of construct- 
ing a separate protective steel shell; 

4) considerable reduction in or total dispensation 
with the reinforcement of concrete biological shield 
structures. 


The Reactor Vessel 





The thickness of concrete shields (usually 2-3 m), 
which is determined with respect to the maximum 
biologically allowable radiation level, is considerably 
greater than the thickness required for the static strength 
and rigidity of the chamber which is formed by the side 
and top reactor shielding. 

In this case, the compressive stresses in protective 
concrete are several times smaller than the ultimate 
strength of the material. For instance, at the lower 
part of a 20 m high lateral shield made of ordinary 
concrete, the compressive stress is 4.6 kg/ cm’. 

If such a shield were made of concrete with scrap 
steel, which would have a specific weight of 4.5 tons/ 
m?, the compressive stress would be 9.0 kg/ cm, 

The ultimate compressive strength of ordinary 
Concrete 300 of different compositions is ~ 140 kg/ cm’, 
and for concrete of different compositions with scrap 


steel, the ultimate compressive strength is 110-150 
kg/cm* [1]. Thus, the compressive strength of concrete 
is not fully utilized in lateral shields, even in the case 
of very large reactors. For the most adverse thermal 
stress distribution, the safety factor will be reduced; 
however, it will still be equal to 3-7. 

The safety factors in compression for the reactor 
top shield are somewhat smaller although still consider- 
able. 

Tensile stresses can arise as a result of temperature 
drops only in the outside shield layers. These stresses 
are small, and, in large reactors, they are usually 
borne by a small part of the steel reinforcement. 

Ragardless of the fact that the concrete strength 
is utilized to such a small extent, a metallic (usually 
steel) vessel, which is subjected to considerable tem- 
perature stresses, is mounted inside the concrete 
chamber. This vessel is designed for possible acciden- 
tal and, in certain reactor types, for operational internal 
liquid or gas pressures. In the majority of reactors, the 
pressure in the vessel does.not exceed 15 atmos. The 
thermal stresses in the vessel can be considerable. 

While the concrete in the lateral biological shield 
is mainly subjected to compression, the material of 
the steel vessel is under tensile stresses. Naturally, in 
the case of such a difference in stresses, the question 
arises whether concrete biological shield structures 
could be additionally stressed. By changing the shape 
and the dimensions of the chamber which is formed by 
the shielding, the chamber can be used as a pressurized 
vessel. 

The maximum allowable tensile stress for concrete 
is small: for ordinary Concrete 300, it is equal to 10.5 
kg/cm’. The ultimate tensile strength of such concrete 
is ~ 30 kg/cm’, and, for instance, for heavy Concrete 
300 with scrap steel, it is 16 kg/cm* (2]. 
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Obviously, without special measures and without 
transmitting the main portion of tensile stresses to the 
steel armature, ordinary concrete biological shield 
structures cannot withstand the assigned loads. One of 
such measures is the compression of concrete by means 
of the so-called prestressed (usually by means of jacks) 
armature (rods, high-strength wire, or steel cable). In 
recent years, prestressed concrete has been used in 
practice on an ever increasing scale. 

However, in using prestressed reinforced concrete 
for reactor vessels, the following difficulties arise: 

1. If a massive concrete biological shield is 
brought close to the reactor core, the thermal stresses 
in the shield become much larger. 

2. Since reliable data on the behavior of structures 
made of prestressed reinforced concrete under conditions 
of violent temperature changes are not available, it is 
obviously necessary to secure large safety factors. 

3. If the nonuniformity in heating is taken into 
account, the design of such a structure becomes very 
complicated. In this, possible errors for an isotropic 
body can yield not only quantitatively inaccurate results, 
but also qualitative deviations from the actual values 
(changes in the sign of stresses in concrete). 

4, The great number of openings in the walls and 
the bottom of the vessel require a number of inserts, 
which complicates the construction of reinforced con- 
crete structures. 

These difficulties forced the joint Soviet-C zecho- 
slovak Commission of Experts to reject an interesting 
project of a prestressed reinforced concrete vessel, 
which has been developed by Academician Békhine of 


the Czechoslovak Academy of Sciences for the Czecho- 
slovak Atomic Power Station A-1. Nevertheless, the 
commission of experts confirmed, in principle, the 
possibility of constructing such a vessel. 

A perfectly gas-proof concrete cannot be prepared 
in practice, and therefore, the inside of the chamber 
must be lined with welded steel sheets which are fitted 
to the concrete and are not designed to withstand gas 
pressure. Such a steel casing can simultaneously serve 
as the internal lining in making a concrete vessel. 

The advantages of a vessel made of prestressed 
reinforced concrete consist in the over-all saving of 
steel (taking into account the prestressed armature of 
reinforced concrete and the internal lining) and in being 
able to dispense with the complex work in the prepa- 
ration of curved sheets, in their welding at the place 
of assembly, and in releasing thermal stresses arising 
in welding individual vessel parts, These advantages 
determined the application of such a vessel in French 
horizontal reactors G-2 and G-3 (Figs. 1 and 2). 

The main reasons for not using an ordinary metallic 
vessel are the expected difficulties in welding and in 
the stress release of a cylindrical vessel 14 m in diam- 
eter and 18 m long. Calculations show that cracks 
in prestressed reinforced concrete in such a structure 
can arise only for internal pressures which are twice as 
large as the calculated pressure, i.e., for 30 atmos. 
The authors of the project consider that the provision 
of the necessary openings in the walls and the bottom 
of vessels is simpler when the latter are made of 
reinforced concrete. On the inside, the vessel has the 
shape of a cylinder which, at the bottom, passes into 
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Fig. 1. Longitudinal cross section of the G-2 reactor vessel. 1) CO, 
inlet; 2) ionization chambers; 3) CO, outlet; 4) longitudinal cable; 
5) transverse cables. 





two inclined planes. In the authors’ opinion, this 
makes the prestressing of reinforced concrete vessels 
easier. 

Steel cables, made of strands of high-strength 
wire, serve as the armature securing the prestressing 
of concrete. Each cable was strained with a force of 
1200 tons by means of jacks. Besides stressed cables, 
band armatures (hoops) were provided in the vessel concrete. 

In ordet to provide an air-tight seal, the internal 
part of the vessel was lined with welded sheets of soft 
steel. This lining was fastened to the vessel concrete 
by means of a number of welded stays, and it permitted 
the mounting of the reactor inside the vessel simulta - 
neously with the construction of its reinforced concrete 
part. The end faces of the vessel were covered with 
spherical reinforced concrete domes, each having a 
thickness of 3 m. The longitudinal forces acting on 
the end-face domes were carried by the longitudinal 
cable armature. 

An over-all amount of 500 tons of armature of 
ordinary steel(besides the 2500 tons of steel for the 
base) and 760 tons of steel wire [2, 3] were used in 
constructing the vessel. 

According to data in [4], due to the saving 
achieved by usinga shield-vessel made of prestressed 
reinforced concrete, the construction costs were sub- 
stantially reduced. 

The above-mentioned indicates that steel can be 
saved by using reactor core chambers made of prestressed 


























Fig. 2. Transverse cross section of the G-2 reactor vessel. 


reinforced concrete, which serves as the biological 
reactor shield at the same time. However, the absence 
of convincing design estimates does not permit definite 
conclusions about the efficiency of this method. 

In reactors with a graphite moderator and a water 
coolant, which flows through pipes (i.e., in reactors 
without internal gas or liquid pressure exerted directly 
on the vessel), the construction of vessels made exclu- 
sively of reinforced concrete, which serves as the bio- 
logical shield at the same time (without prestressed 
armature) is, in a number of cases, technically possible 
and provides a saving in steel. In this, the internal 
concrete layer(~ 50 cm thick) should be heat resist- 
ant, and the internal chamber surface should be lined 
with thin steel sheets in order to provide an air-tight 
seal; in the case where a steel thermal shield is pro- 
vided, it should be combined with this lining. 


The Thermal Shield 





As a result of radiation heat release, the internal 
part of the reactor concrete shield can be heated to 
a considerable extent during operation. A thermal 
shield is used with the aim of reducing excessive heat- 
ing of concrete (which causes large temperature stresses 
and sharply reduces the content of chemically bound water 
in concrete, due to which the protective properties of 
concrete deteriorate). 

In the majority of cases, this shield consists of 
a low-carbon steel screen 15-20 cm thick, In individ- 
ual cases, it is made of lead or of boron-containing 
materials. If considerable amounts of energy are 
released, thermal and biological shields are provided 
with artificial cooling. The weight of thermal shields 
is considerable. For instance, in the reactors of the 
English atomic power stations in Calder Hall and 
Berkeley, the weight of the steel thermal shields is 1925 
and 4340 tons, respectively. 

Sometimes the design of thermal shields is 
very complicated. For instance, in the experimental 
reactor FR-2 (German Federal Republic), the alumi- 
num reactor vessel is surrounded with a thermal shield 
consisting of boral (7 mm) and steel (15 mm) layers 
which are cooled by circulating water, and of lead 
(60 mm), steel (35 mm), and iron (180 mm) layers, 
which are cooled by air. 

However, the experience gained in the use of 
reactors and the results of a number of investigations 
give rise to serious doubts as to whether the applica - 
tion of thermal shields is justified. Their efficiency 
in bringing about a considerable reduction in the 
undesirable heating of concrete is not great. 

Actually, a thermal shield can hardly affect the 
amount of heat transmitted from the reactor vessel to 
the concrete shield, and the amount of heat released 
in the concrete shield by the absorption of radiation 
in it is small. In the majority of cases, the procedure 
of constructing a thermal shield for a comparatively 
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small reduction in radiation fluxes incident on the 
concrete shield is not efficient. 

In a number of papers from abroad, based on the 
experience in the design and operation of reactors, 
the authors take a skeptical attitude toward the advis- 
ability of using thermal shields. For instance, it is 
noted in [5] that, in the experimental reactor of the 
Brookhaven laboratory (USA), which is not provided 
with thermal shielding, no noticeable damage has 
been observed in the concrete shield, which is located 
in the immediate vicinity of the reactor core. Tests 
of concrete specimens cut from this shield showed that 
the strength and the moisture content of specimens 
have not decreased (it should be noted that the concrete 


shield of this reactor has an external bitumen insulation). 


In [6-8], it is also noted that dispensing with a thermal 
shield is possible and advisable from the point of view 
of saving material. 

The use of thermal shields is mainly connected 
with the problem of stresses caused by heat release in 
the shield. Let us consider this problem in more 
detail. 


Influence of Temperature on Concrete Strength 

In a concrete shield, a relatively thin internal 
layer is subjected to highly intensive heating. 

According to data from [6], the thickness of this 
layer is only a few centimeters. The temperature 
distribution throughout the concrete shield thickness 
is shown in Fig. 3[9]. It follows from the figure that 
the use of a thermal shield shifts the temperature 
maximum toward the internal shield surface by only 
0.15-0.20 m, i.e., only a distance equal to the usual 
thickness of steel thermal shields. 

It was reported in [17] that, for a steady action of 
temperatures from 65-95° C on ordinary concrete, the 
deterioration of concrete strength is insignificant, and 
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Fig. 3. Temperature distribution throughout the biological 


shield thickness, 1) without a thermal shield; 2) with a 
thermal shield. 
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even for temperatures in the 260-315° C range, the 
deterioration of structural properties of concrete in a 
number of shielding structures is considered to be ad- 
missible in the case of uniform heating. For limonite 
concrete, which is comparatively seldom used, a maxi- 
mum allowable temperature limit of ~ 200° C is recom- 
mended regardless of the fact that, according to the 
performed experiments, the ultimate compressive 
strength of this concrete is reduced by only 5% at a 
temperature of 350° C. If the temperature is raised to 
~ 600° C, the ultimate strength decreases by 35%. 
According to data from [7], for intensive and frequent 
temperature variations, the deterioration of concrete 
strength is two to three times greater than that result- 
ing from a constant action of high temperatures. 

Some other authors recommend somewhat lower 
maximum allowable concrete temperatures. 

In [10], it is recommended to limit the temper- 
ature of ordinary concrete to 250° C. In this, it should 
be taken into account that the concrete strength can 
be reduced by 20-25% in prolonged heating in the 
temperature region of 100-250° C. It is recommended 
that the temperature of barite and limonite concretes 
be limited to 200° C, since, at temperatures from 
100 to 200° C, capillary cracks appear in these con- 
cretes, and, at temperatures above 200° C, deformations 
which are more dangerous to the safety of structures 
can appear. 

It is reported in [8] that limonite concrete can be 
used for the shielding of reactors at temperatures not 
exceeding 100° C. 

The maximum allowable compressive strength 
of concrete at high temperatures (even for limonite 
concretes which are the weakest under high-temperature 
conditions) are considerably larger than the compressive 
stresses acting in lateral reactor biological shields. 

For instance, the ultimate compressive strength of 28 - 
day old limonite concrete at a high temperature 
(540° C) is 47 kg/cn? even in the case of sharp and 
frequent temperature variations. At the same time, 

a compressive stress of 5.4 kg/cm? will be acting in 
a 20 m high shield wall made of limonite concrete 
(1). 

If temperature stresses arising in concrete shields 
are taken into account, the safety factor for ordinary 
(nonlimonite) concrete will be not less than 3 or 4. 
Thus, the ultimate thermal stresses in the outside shield 
layer (usually carried by the armature) and the phys- 
icochemical processes in the heated concrete mass 
rather than the compressive stresses should be taken as 
the strength criteria for concrete shields at higher 
temperatures. In the worst case, both these factors can 
cause only nonthrough capillary cracks in the outside 
stretched concrete layer, which do not affect the 
compactness and the efficiency of concrete structures. 

The top concrete biological shield of reactors 
represents a separate problem. In this case, the concrete 





TABLE 1. Average Moisture Loss in Heating Different Types of Concrete 
to 100° C over a period of Seven Days. 
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(usually demountable slabs) is subject to bending and, 
consequently, tension, if it does not lie on a rigid 
frame. 

The ultimate tensile strength of concrete is low, 
especially at higher temperatures. Thus, according to 
data from [7], in heating 120-day old limonite con- 
crete from 85 to 350° C, the ultimate tensile strength 
is reduced by 30.5%. Therefore, concrete in the top 
shield, which is subject to bending, should be insulated 
against intensive heat, or an adequate reinforcement 
of the extension zone must be secured. 


Effect of Capillary Cracks on the Protective Properties 
of Concrete 

Experiments and calculations [8] show that cracks 
in concrete do not have geometrically smooth surfaces 
which would permit a penetration of radiation through 
the shield parallel to these surfaces. The surface of 
these cracks is rough and broken. In concrete, even a 
through crack 1.6 mm wide has a minimum deviation of 
5 mm from a straight line [8]. The intensity of direct 
radiation through such a crack is small, and, on the 
whole, the radiation intensity at the outside shield sur- 
face in the capillary crack region is approximately the 
same as the intensity over the entire concrete surface. 

Thus, the following conclusions on thermal 
stresses in concrete can be drawn: 

1. In using ordinary high-quality concrete for 
biological protection, even in the case of powerful 
reactors operating under steady conditions, there is 
no danger that the concrete shielding would collapse 
or be damaged due to thermal stresses, including the 
case where thermal shields between the reactor vessel 
and the biological shield and a special cooling arrange- 
ment are not provided. 

2. In using limonite concrete for the internal 
biological shield layer in reactors where the vessel 
temperature is in the range from 150 to 300° C, thermal 
shields or a cooling arrangement are apparently 
necessary. In certain cases, both these measures are 
necessary. Therefore, the use of limonite concrete 
for the internal shield layer is undesirable. 

3. In the case of frequent and sharp temperature 
variations in the shield, it is advisable to apply thermal 
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shielding or cooling of concrete or to use special heat- 
resistant types of concrete for the internal portion of the 
biological shield. 

4. For the top concrete shield, which is subject 
to bending stresses, each individual case requires a 
study of the question whether it would be more effi- 
cient to lower the concrete temperature by installing 
a thermal shield or to provide suitable reinforcement 
of the concrete zone subject to tensile stresses. 


Effect of Temperature on the Composition and Protec- 





tive Properties of Concrete 





heating concrete, its content of free and 
chemically bound water decreases. For different 
types of cement and fillers, and also for different 
concrete ages, these reductions are of unequal mag- 


nitude. The resultant deterioration of the ability of 
concrete to provide protection from neutron radiation 
along with a decrease in concrete strength is the main 
reason for the tendency to reduce the concrete temper- 
ature and to use special thermal shielding. 
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Fig. 4. Water content in cement 
gravel in dependence on temper- 
ature. 1) Magnesium hydroxychloride 
cement; 2) Portland cement. 





TABLE 2. Temperature Dependence of the Amount of 
Chemically Bound Water in Cement Gravel 





Relative water 
Content, % 


Specimen tem-' Amount of 
perature, “C 


ater Yb 
re weight 


50 22,3 100 
100 17,8 80 
200 12,2 59 
300 10,4 45,5 
400 8,0 36,3 
500 3,0 13,6 
600 2,1 
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Figure 4 provides data on the water content in 
cement gravel in dependence on temperature [8], 
Table 1 provides data on the weight loss due to water 
evaporation in different concretes heated to 100° C 
(11), and Table 2 gives data characterizing the temper- 
ature dependence of the amount of chemically bound 
water in cement gravel consisting of ordinary, sand 
and Portland cement [12]. 

It is reported in[13 and 14] that limonite ores 
can be considered as a water-containing material until 
their temperature reaches ~ 190-200° C. 

Table 3 provides data obtained in investigating 
the moderation of neutrons in slabs made of concrete 
of different compositions and heated to different tem- 
peratures [15]. It follows from these data that the 
neutron relaxation length increases with an increase 
in concrete temperature. Similar data were obtained 
for ordinary and magnetite concretes [11]. Thus, with 
respect to moisture content, the use of thermal shield- 
ing or of biological shield cooling (or the simultaneous 
application of these measures) is justified only when 
the temperature of ordinary concrete exceeds 200-250° 
C or when the temperature of limonite concrete exceeds 
100-150° C (with a certain margin). 


Effect of Radiation on Concrete Strength 

No systematic research data on the effect of radia- 
tion on the mechanical properties of concrete are 
available. Data on failure or damage in the concrete 
of biological shields of nuclear reactors or accelerators 
under the influence of radiation have not been recorded 
in practical work withnuclear reactors for peaceful 
purposes. Individual cases of shielding concrete defor- 
mations (capillary cracks) are mainly the consequence 
of the usual shrinkage or thermal stresses). 

In experiments performed on the VERO reactor 
[15], it was established that ordinary concrete in shields 
can be utilized in operation over a period of at least 
ten years for a flux of 104 neutrons/cm*sec. Experi- 
ments [16] showed that also concrete with Portland 
cement can sustain a certain damage if it is subjected 
to an integral flux of over 3 x 10 neutron/cm® sec 
(or if subjected toy radiation). (These data are based 
on a very small number of experiments and are not 
conclusive). 





706 


Comparative investigations of concrete specimens 
with a water-cement ratio of 1:3 and dimensions equal 
to5 x 5 x 20 cm which were irradiated in a reactor 
(for a y -radiation dose rate of 10° r/hr) and of speci- 
mens which were not irradiated showed that concrete 
is resistant to radiation [8]. These and a number of 
other data indicate that ordinary high-quality concrete 
has a good resistance to radiation. In individual cases, 
the omission of thermal shields will necessitate the use 
of somewhat thicker concrete biological shields; this, 
as a rule, is simpler and less expensive than the con- 
struction of thermal shields, and it provides a substantial 
saving in metal. 

All the above-mentioned indicates that, in many 
cases, the necessity of providing thermal shields cannot 
be substantiated if a concrete biological shield of good 
quality is available. The application of special thermal 
shielding can be justified only in the case of high- 
temperature reactors. 

Until the present time, heat-resistant concrete 
has not been used in reactor construction, although it 
has been mentioned in individual papers [17] that its 
application would be possible and efficient. 

Considering that the maximum temperature 
point in concrete is at a distance of 0.35-0.40 m from 
the internal surface in the case where a thermal 
shield is not available, a layer of heat-resistant con- 
crete 0.5-0.6 m thick would probably be sufficient. 
Beyond this layer, the effect of temperature and of the 
temperature drop, even in high-temperature reactors, 
would probably not affect the concrete shield strength 
(except in the region subjected to tensile stresses) and 
would not cause a considerable reduction in the moisture 
content. This relatively thin layer of heat-resistant 
concrete (even if it contains admixtures of finely 
ground red brick or a coarse filler of this brick) can 
reliably withstand temperatures of up to 700° C [12]. 

The somewhat lower specific weight of heat-re- 
sistant concrete in comparison with that of ordinary 
concrete (1.7 tons/ m® versus 2.3 tons fm’) will, in the 
worst case, necessitate a slight increase in the shield 
thickness (by only 13-14 cm). Such an increase in the 
concrete shield thickness will completely compensate 
the somewhat lower (in comparison with ordinary 
concrete) amount of chemically bound water in the 
heat-resistant concrete layer. Increasing the concrete 
shield thickness is obviously much less expensive than 
the use of a metallic thermal shield. 


The Reactor Chamber 





In using reinforced air-tight chambers which are 
formed by the concrete biological shields of the reactor 
proper as well as of steam generators (in power plants) 
as a protective shell, the lateral and the top concrete 
biological shields must be adequately designed for 
accidental excess internal pressures (usually not exceeding 
3-3.5 atmos gauge). The requirement for air-tight 





TABLE 3. Absorption of Fast Neutrons by Heated Concrete Slabs 





Reinforced al Magnetite-limo- Magnetite Ordinary 
: nite concre 
Drying temperature, °C nite concrete coficrete concrete 


gem | A, cm 1.g/cm'| 4.cm +.g/cm}| A, cm 1,g/cm*| A, cm 








Initial 2, 12,06 
2, 12,80 
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seals involves certain difficulties if openings, doors, explosion-proof shield at the same time. In order to 
etc. are present, however, this requirement can be provide a better air tightness, the protective chamber 
satisfied. is lined with ordinary sheet steel on the inside. All 
If the reactor is built in a location with a mild the reinforced concrete structures of the lateral and 
climate, the reactor shell can be used as the housing top biological shields are designed for an excess 
for the entire plant, which obviously provides a saving internal pressure of 3.5 atmos. A protective steel dome 
in material. In locations with a moderate or northern is mounted immediately above the reactor. 
climate, the use of a protective chamber calculated In order to increase the volume of the space to be 
to withstand internal pressures and the simultaneous filled with radioactive vapor in case of accident (for 
building of a reactor housing of the usual type can be the purpose. of reducing the vapor pressure), all the 
more economical than the construction of a steel compartments of the first loop, which are adjacent 
shell with heat insulation. to the reactor chamber are connected to it and to each 
An example of the application of concrete other with membrane valves which operate when the 
shielding for these purposes is a boiling power reactor pressure in any of the adjacent sections rises. 
with a thermal power of 300 Mw (Fig. 5), which is It was found that the utilization of the biological 
under construction in the USSR. In this reactor, the shield as a protective shield in case of accident 
reinforced concrete biological shield is used as the necessitated an increase of 250-300 tons in the rein- 
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Fig. 5. Vertical cross section of the housing for the Soviet 
Boiling Reactor. 1) Bridge crane; 2) loading bridge; 3) reactor 
dome; 4) reactor; FC) chambers to be filled with radioactive 
vapor in case of accident. 





forcement of concrete shielding structures, which amount- 
ed to ~ 50% of the armature weight in ordinary 
nonreinforced biological shielding and to 25-30% of 

the over-all amount of the armature material for all 
reinforced concrete housing structures, which is equal 

to ~ 1000 tons. 

In order to determine the amount of metal to be 
used in the construction of a protective steel shell 
around the entire reactor plant, blueprints for the 
project of a spherical protective steel shell have been 
completed. Since the space which is filled with radio- 
active vapor in case of accident is larger than the space 
used for this purpose when only a protective concrete 
chamber is available, the theoretical pressure inside the 
shell was reduced to 2 atmos gauge. The diagram of 
such a sheet steel shell, which has a diameter of 50 m 
and a thickness of 35 mm, is shown in Fig. 6. The 
weight of sheet metal to be used in the construction of 
the shell is 2150 tons, 

By comparing the two variants, it can be seen that, 
in using the biological shield structure as an accident- 
proof chamber, a 2.5-fold saving of metal in building 
construction can be achieved. It should also be taken 
into account that the variant with a steel shell requires 
a considerable amount of sheet steel (which is less 
readily available than structural armature steel) and 
that, moreover, the work in constructing the reactor 
housing becomes more complicated, especially if one 
considers the necessity of using a considerable amount 
of heat-insulation material for the entire steel shell 
even if the reactor is built in the moderate climatic 
belt of the USSR. 


The Reinforcement of Concrete 





As a rule, steel reinforcement is provided in the 
external zone of concrete biological shields. A con- 
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Fig. 6. Vertical cross section of the housing for the Soviet 
Boiling Reactor with an outside metallic protective shell 
(a variant of the project). 
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siderable amount of reinforcement steel is required. 
For instance, in the large reactors at the English atomic 
power station at Hinckley Point, this amount is equal 
to several thousand tons. 

The reinforcement of concrete in shields is done 
in the first place because of the danger of excessive 
temperature stresses which would lead to the formation 
of capillary cracks and other concrete deformations. 

In recent years, the reinforcement of concrete was 
based on calculations of stresses acting in concrete [7]. 
However, these calculations are not reliable, since 
sufficient experimental data on the temperature distri- 
bution throughout the thickness of a monolithic concrete 
shield are not available. 

Reinforcement is often provided without any cal- 
culations, It was reported in [8] that, if the tempera- 
ture drop throughout the concrete shield thickness does 
not exceed 30-50° C, special reinforcing of concrete 
is not necessary. For temperature drops of up to 90° C, 
the amount of reinforcing material must be substantially 
increased. In this case, according to the authors’ opinion, 
the costs of additional concrete reinforcement and of 
constructing a special thermal shield should be compared 
in order to secure a most economical solution. For 
temperature drops of over 90° C, the reinforcement of 
concrete is impossible for constructional reasons, and, 
therefore, the temperature drop must be reduced by 
applying thermal shielding or cooling (or by using both 
these measures at the same time). 

In [7], the following principles are recommended 
as guidance in calculating thermal stresses in reactor 
concrete shields: 

1, If the cracking of nonreinforced concrete is to 
be prevented, the average shield temperature must not 
exceed the temperature of the shield cold side by more 
than 11° C. 

2. In reinforcing the concrete shield region which 
is subjected to tensile stresses, the over-all allowable 
temperature drop (between the hot and the cold sides) 
can be increased to 67° C. For a nonlinear temperature 
distribution throughout the concrete thickness this drop 
must not exceed 55° C, 

3. In the case of frequent temperature variations, 
itis recommended to limit the temperature drop to 
22 -28° C. 

4, If the temperature drop increases to a considerable 
extent, a substantially heavier reinforcement must be 
provided. For instance, if the temperature drop increases 
from 70 to 100° C, the reinforcement must be increased 
by 0.2-1.6% (in order to avoid the formation of 
capillary cracks). Therefore, in certain cases, the 
provision of a thicker inner steel thermal shield or an 
intensification of artificial cooling of concrete is more 
economical than the use of a heavier reinforcement. 

5. If pipes carrying a coolant at a temperature 
of 150-315° C pass through the concrete shield, an 
adequate cooling of the shield is sometimes impossible. 





In this.case, the concrete is usually enclosed in a 
carrying steel structure. 

In the USSR, concrete and reinforced concrete 
shielding elements operating at a temperature of 100° C 
are designed according to the "Norms and Technical 
Requirements for the Design of Concrete and Reinforced 
Concrete Structures." For higher temperatures, the 
design is done according to the "Provisory Standards 
for the Design of Heat-Resistant Reinforced Concrete 
Structures." Nevertheless, up to now, sufficiently 
accurate methods of determining thermal stresses in 
unilateral heating of such massive concrete and reinforced 
concrete structures as biological shield structures ~ 2-3m 
thick have not yet been developed. 

From all the above data, the following conclusions 
can be drawn: 

1. In connection with the fact that only a limited 
amount of data on the temperature distribution through- 
out the thickness of concrete biological shields is 
available, calculations of thermal stresses and of the 
necessary percentage of reinforcement in shields yield 


only provisory data. Therefore, suitable experiments 
and the development of methods for accurate calcula- 
tions of stresses and concrete reinforcement are 
necessary. 

2. A number ofreactors which are in operation, 
including high-temperature reactors, are not provided 
with reinforced concrete or thermal shielding. How- 
ever, no deformations have been observed in the 
shielding concrete in these reactors besides the appear- 
ance, in certain cases, of capillary cracks, which are 
not dangerous. Therefore, the recommendations of 
authors from abroad pertaining to limitations of 
temperature drops in concrete shielding for certain 
given amounts of reinforcement contain unjustifiably 
large safety factors, and they are not based on reliable 
calculations. 

3. The slight reinforcement of the concrete 
biological shielding in the majority of reactors with 
low and medium temperatures can be dispensed with, 
and the reinforcement of the internal compressed 
protective zone is superfluous in all cases. 
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An account is given of preliminary results in a study of the influence of a magnetic field on electron diffusion 
in a plasma. An intermittent increase was observed in the ratio of the electron current in the probe to the 

ion current at some critical value of the magnetic field intensity. According to the preliminary data, the 
critical magnetic field changes proportionally to the gas pressure. These facts evidently indicate the existence 
of two qualitatively different mechanisms for the transverse movement of the electrons, one of which is diffusion 


by collision. 


One of the methods permitting the determination 


of the velocity with which electrons travel in the plasma 


across the magnetic field is the measurement of the 
ratio of the electron current on a probe to the ion cur- 
rent in regions of the discharge where it is possible to 
neglect the ionization. On the basis of such measure- 
ments, D. Bohm [1] formulated his hypothesis on the 
anomalously large electron transverse mobility. Al- 
though the mechanism of the Bohm diffusion of elec- 
trons is not clear, nevertheless, it can be assumed that 
it does not always take place, and that under some 
conditions the velocity of the ordinary electron diffu- 
sion by collision is sufficient to compensate the ionic 
space charge. 

In order to verify this hypothesis, measurements 
were carried out, in the present work, of the ratio of 
the electron current in a probe to the ion current for 


Fig. 1. Discharge chamber. 


different magnetic field intensities in the interval of 
400-4500 oe. 


Apparatus and Measuring Circuit 





A diagram of the discharge apparatus is shown in 
Fig. 1. The discharge chamber consists of a cavity 
surrounded by a thick-walled copper cylinder C, whose 
end faces are closed by a molybdenum anode A and 
a diaphragm D. A heated cathode prepared by B. N. 
Makov was used in this work. The cathode proper 
consisted of a tungsten cylinder K, 5 mm in diameter, 
mounted on a special bracket on the cathode block B 
opposite an opening in the diaphragm 3 mm indiam- 
eter. The cathode was heated by the thermionic 
current from a tungsten filament N. In the anode, at 
a distance of 10 mm from its center, six holes were 
drilled, their diameters being 0.8, 1, 2, 3, 4, and 2 
mm. Ata distance of ~ 0.5 mm from the rear surface 
of the anode were mounted, opposite each opening, six 











Fig. 2. Measuring circuit . 





collectors S, which were used to take the probe charac- 
teristics. The collectors were fastened to the body of 
the shield G through glass insulators. Nitrogen was 
used as the working gas, and was let into the discharge 
cavity through tube F. 

This arrangement was placed inside a metal 
vacuum chamber, which was set in the gap of the 
electromagnet and was evacuated by an oil diffusion 
pump which maintained during the continuous flow of 
the nitrogen in the discharge chamber, a vacuum 
inside it of ~ 10°°mm Hg. The pressure in the discharge 
chamber was approximately two orders greater. 

Figure 2 shows the measuring circuit. Between 
the probe S and the anode A were connected in 
succession, a saw-tooth generator G and a measuring 
resistance R of 1 kohm. The dc component of the 
generator output voltage was compensated by a battery 
B, so that the potential of the probe could be changed 
symmetrically with respect to the anode. The signal 
was fed from the measuring resistance to the input of 
a direct-current amplifier UIPP-2, the output of which 
was applied directly to the plates of an £O-7 oscillo- 
graph, 

To avoid linear distortion of the probe currents 
observed when the circuit operated at a high sensitivity 
the probe potential was not measured directly, i.e., 
the horizontal deflection of the oscillograph beam was 
obtained from a sweep generator of the oscillograph. 


N 1 (H=4300e,k:=50) M2 (H=8500e k=400) 


During all measurements the amplitude of the 
probe potential with respect to the anode was set at 
~ 15 v. The measurements were made at a frequency 
of 50 cps. In the case of strong magnetic fields 
intensive oscillations always arose in the discharge ,and 
the electronic part of the probe characteristic was 
very strongly modulated by oscillations at some tens 
of kilocycles per second. Under these conditions the 
measuring resistance was shunted by a capacitance of 
0.1 uf by means of switch Ky, and thus, the cliarac- 
teristic was averaged over the high-frequency compo- 
nent of the probe current. 

In those cases when the electron current exceeded 
the ion current by more than 20-30 times a more 
correct measurement of the ratio of the electron current 
to the ion current was obtained by measuring the 
latter under a greater amplification than that of the 
electron current. In this way, the opening of switch K, 
connected in the detector D, and the electron current 
was detected to remove the overload on the amplifier. 

Switch Kg, shunting the input of the amplifier, 
was used to obtain the zero-current line during the 
photographing of the probe characteristics. 


Results of the Measurements 





Figure 3 presents a number of oscillographs of the 
probe characteristics corresponding to a discharge 
current of 1 amp, a discharge voltage of 180 v, and a 
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Fig. 3. Oscillograms of probe discharge characteristics for different magnetic -field intensitites 
H (k is the amplification factor). 





pressure in the vacuum chamber of 7 - 10° mm Hg. 
In this case, the pressure in the discharge chamber was 
(2-8) -10-* mm Hg. These characteristics were 
obtained for a probe diameter of 3 mm, the remaining 
probes were grounded. Oscillograms Nos. 1-7 were 
taken without the capacitance shunting the measuring 
resistance. 

In examining these oscillograms, one may note 
the absence of oscillations both of the ion current 


(lower branches of the curves) and of the electron cur- 
rent (upper branches of the curves) for a magnetic- 
field intensity of < 3000 oe. However, in a magnetic 
field only a little stronger than 3000 oe (oscillogram 
No. 7) there arises a very strong oscillation in the 
electron current of the probe. At the same time, os- 
cillograms Nos. 1-7 show quite distinctly that there 
is a steady decrease in the ratio of the electron cur- 
rent to the ion current as the magnetic field increases. 
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Fig. 4. Variation of electron current I_ 
and ion current I, in the probe with mag- 
netic-field intensity. Probe diamter Dg= 
3 mm; I, = 1 amp; V, = 180 vVip=T- 
10° mm Hg. 
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Fig. 6. Variation of electron current 
I. and ion current L, in the probe with 
magnetic -field intensity. Dg = 3 mm; 
I, = 0.6 amp; Vp = 210 vi p= 7- 10° 
mm Hg. 
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Fig. 5. Variation of ratio I-/1, with 
magnetic-field intensity (from data 
of Figs. 3 and 4). Dg = 3 mm; I= 

1 amp; Vp = 180 v; p= 7: 10° mm 
Hg. 
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Fig. 7. Variation of ratio I_/1, with 
magnetic-field intensity (from data 
of Fig. 6). 











The behavior of the electron current I and the 
ion current I, in the probe with changes in the mag- 
netic field is more clearly illustrated by Fig. 4, ob- 
tained on the basis of oscillograms Nos. 1-8. Of great- 
est interest is the intermittent increase in the 
electron current at a magnetic-field intensity of ~ 
3050 oe, while at the instant of the jump, the ion 
current essentially does not change. 

Figure 5 shows the variation of the ratio I /1, 
with the magnetic field intensity. At first, this ratio 
falls from 11 at 450 oe to 4 at 2500 oe. In the 
interval of 3000-3050 oe it rises from 5 to 17, if the 
electron current is measured on the lower curve of 
oscillogram No. 7 (without the capacitance); if the 
same electron current is measured at the averaged 
oscillogram No. 8 (with the capacitance), then the 
ratio I_/L, in a field of 3050 oe intensity attains ~ 10”. 

Figures 6 and 7 show a similar variation (currents 
are given in arbitrary units) for the same probe, but 
at a discharge current of 0.6 amp, discharge voltage of 
210 v, and the same pressure (p = 7 - 10°° mm Hg). 

In this case the magnetic field was varied over a wider 
range and it is seen that the intermittent increase in 
the electron current follows the decrease in the field 
more or less smoothly. 

If the pressure is decreased by half (3.5 - 10 - 
mm Hg) the general character of the probe charac- 
teristics remains approximately the same. However, 
the jump in the electron current accompanying the 
appearance of the strong oscillations takes place at a 
magnetic field intensity of ~ 1500 oe, and not at ~ 
3000 oe, as is the case at a pressure of 7 - 10° mm Hg. 


Discussion of the Results 





The basic experimental results of the present 
work is the discovery that at some critical magnetic 
field there is an intermittent increase in the electron 
current of the probe, this being accompanied by the 
appearance of strong oscillations. The magnitude of 
the critical value of the magnetic field increases 
proportionally to the gas pressure, which evidently 
indicates the important role of the quantities w, T 
or W T_ characterizing the magnetization of the 
plasma, and, possibly, predetermining the occurrence 
of a new state of discharge with strong oscillations and 
a large transverse velocity of the electrons. 

A strict quantitative analysis of the results which 
have been obtained is made difficult for a number of 
reasons. Firstly, the design of the probe used in this 
work cannot be regarded as satisfactory, since the edges 
of the probe opening in the anode are enveloped by 
a layer in which there is a negative potential drop at 
the anode, and, as a result, an indefinite collecting 
surface is formed for both electrons and ions. Secondly, 
the theory developed by Bohm [1] for the collecting 
of electrons by the probe in a discharge in a strong 


+ 


magnetic field is valid for sufficiently long discharges 
satisfying the condition 


L> 


a 
V DP 
D_ 

(L is the discharge length, a is the radius of the probe, 
DP is the coefficient of transverse electron diffusion 
in the magnetic field, D is the coefficient of electron 
diffusion without the magnetic field). In the present 
work, unfortunately, the opposite conditions were 
fulfilled. Finally, in Bohm's theory it is assumed for 
the kinetic energy of the ions in the plasma that eT, 
<< eT (T_ is the electron temperature), and therefore 
saturation of the electron current should have been 
observed at probe potentials that were positive with 
respect to the plasma by a value of ~ T,. As may be 
seen from the oscillograms of Fig. 3, under the con- 
ditions of this experiment, saturation of the electron 
current is not observed, which evidently is a conse- 
quence of the relatively high ion energy [2]. 

If we ignore all these reservations and nevertheless 
use for estimating the ratio I_/1, the formula of Bohm 





then, after substituting the values a = 0.15 cm, VM)= 
~4,T =lev,1+(T,/T_)™1,we obtain for H = 3000 
oe the value L/1, ™3.5, which is close to the 
measured ratio ~ 4. It should be noted that the ex- 
perimental value of I_/1, cannot be determined accu- 
rately also because of the absence of electron current 
saturation. 

Without attaching too much importance to this 
estimation, but keeping in mind the fact that there 
is no doubt that there is an intermittent increase in the 
electron current, one may consider that the results ob- 
tained here, along with the results of [3], confirm the 
existence, to some degree of two qualitatively different 
mechanisms for the transverse movement of the electrons. 
The first mechanism, acting in the region below the 
critical value of the magnetic field, is evidently 
ordinary diffusion resulting from collisions. The 
second mechanism is, as before, unclear, but it may 
be hoped that more detailed and careful investigation 
of the observed critical transition in the discharge state 
will permit an explanation of its fundamental properties. 

In conclusion, I wish to thank P. M. Morozov for 
his assistance and interest in this work, and also S. 
Sinotov, for taking part in the measurements. 
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The experimental data concerning the anomalously large mobility of electrons across a magnetic field is dis- 
cussed, It is shown that the distribution of the concentration of the secondary plasma of a discharge with an 
incandescent cathode is practically independent of the electron transverse diffusion coefficient, and therefore 
cannot be used to explain the mechanism of diffusion. An estimation is made of the electron diffusion coeffi- 
cient from the value of the electron current density at the anode, which confirms the presence of an anomalous- 


ly large transverse mobility. 


As a result of the study of experimental data on a 
discharge with an incandescent cathode in a strong lon- 
gitudinal magnetic field, D. Bohm [1] suggested the 
hypothesis that under the conditions of such a discharge, 
a new, uninvestigated mechanism for the transverse 
movement of the electrons is considerably more effective 
than the diffusion by collision. 

In particular, one of the arguments in favor of the 
anomalously large transverse mobility of electrons is 
the result of the analysis of probe characteristics. 
Assuming that the diffusion of electrons takes place by 
collision, Bohm obtained for the ratio of the electron 
to ion current in the probe I_/L, the value 4-5, while 
the I_-/I, ratio measured under the corresponding condi - 
tions attained 20-35. Although Bohm's calculations 
do not lay claim to great accuracy, nevertheless, it is 
not very likely that they can give a reduced value of 
I./1,, and therefore, they should not be considered. 

Another experimental fact which Bohm used as 
evidence of the anomalously large transverse mobility 
of electrons is that the decrease in the concentration 
of the plasma with the distance from the beam of 
primary electrons takes place more slowly than expected. 
An approximate theoretical analysis of the distribution 
of the secondary plasma concentration (outside the 
primary electron beam) has been carried out in the two- 
dimensional case. Using conditions of quasi-neutrality 
treating the transverse motion of electrons and ions 
as a diffusion process, and neglecting the ionization 
in the secondary plasma, Bohm obtained the following 
expression for the characteristic length x9 on which 
the concentration decreases e times as one moves away 
from the primary beam 


= 
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where T, and T_ are the temperatures of the ions and 
electrons in units of potential; pP and DP are the ion 
and electron transverse diffusion coefficients. The 
coefficients 8 and y relate the ion current density at 
the anode . and the electron current density at the 
anode j® to the plasma concentration n: 


7% = Bn; (2) 
pi=yn. 


Assuming further that 
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Bohm obtained, in place of (1), the expression 
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Therefore, according to (4), the measurement of 
X9 allows one to estimate the electron transverse mobil- 
ity coefficient. Comparing the measured value of Xo 
with the value calculated from (4) and assuming that 
ordinary diffusion takes place, Bohm discovered a lack 
of agreement between theory and experiment. It tumed 
out that the electron diffusion coefficient calculated 
from (4) exceeds the collision diffusion coefficient 
by two orders of magnitude. 

In his work A. Simon [2, 3] obtained an expression 
for Xp substantially different from (4) 
L DP 
Ly & — — 


; (5) 





where D,,, is the longitudinal ambipolar diffusion 
coefficient. 

Comparing the measured value of Xp to that cal- 
culated from (5), Simon found satisfactory agreement. 
It turned out, also, that as the magnetic field intensity 
increased, Xp decreased proportionally to 1/H. On this 
basis, Simon believed that in the experiments he was 
considering, the diffusion of both the ions and the 
electrons took place in accordance with the well-known 
ideas, and that it was not necessary to assume the 
existence of an anomalously large electron transverse 
mobility as Bohm had done. 

Simon, however, ignored the results of the above- 
mentioned probe measurements. Moreover, one can 
hardly, in general, determine the electron diffusion 
coefficient on the basis of a measurement of the charac- 
teristic length xp under conditions of a discharge of the 
type which Simon and Bohm had in mind. 

The purpose of the present article is to show that 
under certain conditions the distribution of the secondary 
plasma concentration does not depend on the value of 
the electron diffusion coefficient, but is determined 
mainly by the picture of the ion motion. Owing to this, 
it may be supposed that Bohm and Simon were in 
error when they took the data on the plasma concentra- 
tion distribution to explain the mechanism of the 
transverse diffusion of electrons. 


ON THE DISTRIBUTION OF THE PLASMA CONCENTRA- 
TION IN A DISCHARGE WITH AN INCANDESCENT 
CATHODE IN A LONGITUDINAL MAGNETIC FIELD 


Bohm's formula (4) corresponds the ambipolar 
diffusion of the plasma, to which, apparently, Simon 
first drew attention. It is interesting to note that Bohm 
did not use as his initial assumption the condition of 
ambipolarity of the diffusion, but condition (3) which 
he took as. equivalent to such an assumption. 

It may be stated, however, that in the conditions 
of the discharge being considered (Fig. 1), the diffusion 


100 


Fig. 1. Diagram of discharge chamber. 
K) cathode; A) anode; P) primary beam; 
S) probe. 


cannot be ambipolar. Consequently, if there is ambi- 
polar transverse diffusion outside the primary beam, 
then everywhere outside the primary beam P the re- 
sulting longitudinal current on the anode should be 
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equal to zero, i.e, the equation 


a v_ 


(6) 


should hold, where ¢ is the potential of the plasma 
with respect to the anode; v_ and v4 are the mean 
longitudinal velocities of the electrons and ions 
[(6) is a consequence of the equality of the currents 
at the anode; Pi = a fl At the same time, for the 
transverse diffusion to be ambipolar, the potential ¢, 
as one moves away from the primary beam, should in- 
crease to a value of the order of T,, which is incom- 
patible with (6) for an equipotential anode, since the 
temperature of the electrons T_ in a strong magnetic 
field should decrease comparatively rapidly outside the 
primary beam, owing to the cooling of the electron gas, 
For experimental evidence of the assumption one 
may use the results of the measurement of the currents 
going to the anode in the secondary plasma as obtained 
by the author and shown in Figs. 2 and 3. These figures 
show the current distribution on the anode surface, 
measured by means of a probe in the shape of a disk 
1 mm in diameter, which was placed in the plane of 
the anode and moved perpendicularly to the magnetic 
field. The primary current had a diameter of 5 mm, 


15C 


2 
-I;,ma 


Fig. 2. Current distribution in the 
cross section of the discharge (1 = 
200 ma; V, = 80 v; H = 2300 oe; 
p= 2-10" mm Hg). 





The series of curves presented corresponds to a discharge 
voltage Vp, = 80 v, a discharge current I, = 200 ma, 

a magnetic field intensity H = 2300 oe, and argon 
pressure inside the discharge chamber p = 2 - 10 “> mm 
Hg. The length of the discharge gap L was 10 cm, the 
diameter of the discharge tube was 30 mm. 

In Fig. 2 the distance from the probe to the wall 
of the discharge chamber is laid off along the abscissa 
axis, the ion current in microamperes is laid off along 
the upper part of the ordinate axis, while the electron 
current in milliamperesis laid off along the lower part 
of the ordinate axis. The region occupied by the 
primary beam is marked by a circle. The probe was at 
the potential of the anode. It is seen from the figure 
that at a distance ~ 1 mm from the edge of the primary 
beam the ion current predominates. Moreover, as 
may be seen from the voltage vs current characteristics 
(taken in this experiment and shown in Fig. 3), the ion 
current on the probe is essentially greater than the 
electron current, beginning with a distance of ~ 1 mm 
from the edge of the primary beam. At greater dis- 
tances (Fig. 3b) there is, in effect, a saturated ion current 
on the probe. For example, at the point r= 9 mm, 
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back from the center of the primary beam by ~ 8-9 
mm the electron current constitutes about 6% of the ion 
current. 

Clearly, these facts convincingly indicate that the 
diffusion in the secondary plasma is not ambipolar and 
that the electrons move across the field considerably 
more slowly than the ions. 

On the basis of these data one may consider that 
the electrons in the secondary plasma are, as it were, 
in a closed potential box, and to a first approximation, 
their escape from it may be neglected, i.e., it may be 
assumed that the electron gas is in a Boltzmann equilib- 
rium, and therefore the electron stream in the 
secondary plasma is practically absent. Then, as may 
readily be shown, the characteristic length x» will be 
given by the formula 


L ,/ DP 
neh 2. (7) 


According to this formula, the distribution of the 
concentration of the secondary plasma does not depend 
on the coefficient of electron diffusion, and the physi- 
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Fig. 3. Ion part of the probe characteristics taken at different points of the cross section 
of the discharge. a) close to the primary beam; b) at the periphery. 





cal meaning of it is understandable, if one considers 
that, in the case of ion diffusion, the ratio of the 
mean square displacements across and along the mag- 
netic field is given by 


x3 
L 2 
(z) 
One may also obtain (7) within the framework 
of Bohm's theory if in accordance with the smallness 
of the electron current on the anode, even in compar- 


ison with the ion current, one assumes that the condi- 
tion opposite to conditon (3) 
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is fulfilled. 
Taking into account (8) we obtain, instead of (4), 
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In the case of longitudinal diffusion of ions the order 
of magnitude of the coefficient 6 is 


(10) 


It is not difficult to see that when (10) is substituted 
into (9), one actually obtains (7). It should be noted 
that Simon's formula (5) gives practically the same 
quantitative results as (7), but, as it seems to us, it 
has a less clear physical meaning and was obtained 

by neglecting, in the balance of streams of charged 
particles, terms proportional to the transverse mobility 
of the ions and electrons. The correctness of neglect- 
ing these terms is not quite evident. 

Under conditions in which the distribution of the 
space charge is practically independent of the coeffi - 
cient of electron diffusion, an attempt may be made 
to estimate directly the value of the electron current 


at the anode at some point r. Its order of magnitude 
will be 


: LD?n 
xX (11) 


The ratio of the electron current to the ion current is 


ante 
J? a0, (12) 


The lower probe characteristic of Fig. 3 b(r = 9 mm) 
was obtained at a plasma concentration of ~ 10"! cm™® 
and a neutral gas concentration of ~7'10"cm™*, Atsuch 
a ratio of the concentrations and at T. “ev, the 
collisions between the electrons and the ions do not 
have to be taken into account. Then, in a magnetic 
field H = 2300 oe, the electron diffusion coefficient 
is DP ~10 cm*/ sec. Taking into account the fact 
that the argon ions move freely [4] and have a longi- 
tudinal energy ~ T_ = 1 ev, we obtain v, ™2 - 10° 
cm/sec. Inserting these data into (12) and setting 
L = 10 cm and r = 1 cm, we obtain 
a 

Few 5-404, 

The measured ratio of currents at the anode, as 
has been indicated previously, is 6 * 10~*. Although 
this estimation is very coarse, the observed discrepancy 
(by two orders of magnitude) is clearly not in favor of 
diffusion by collision. 

A more convincing proof of the existence of some 
mechanism of the transverse movement of electrons 
other than the ordinary diffusion is the intermittent 
increase in the ratio of the electron current on the 
probe to the ion current [5], and also the sharp increase 
in the coefficient of transverse diffusion [6]. 

In conclusion I wish to thank I. V. Kurchatov for 
his interest in this work and his constructive criticism. 
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The energy spectrum and the average enemy of the external proton beam of the six-meter synchrocyclotron of 
the Joint Institute for Nuclear Studies have been investigated under various operating conditions. The spectrum 
is described by a Gaussian curve with a dispersion of (2.8:0.3) Mev. The mean energy has been determined to 
within an accuracy of 0.1%, 


Introduction 20°. In front of the magnet there are several brass 


collimators which are essentially vertical slits; these 
collimators define a plane proton beam 0.1-0.7 cm 
in width and 2 cm in height. At the outputof the 


The increase in the beam intensity of the proton 
beam extracted from the synchrocyclotron chamber 
at the Joint Institute for Nuclear Studies [1] has made 
it possible to oumty out a number of enpeciinens in a, 3 ~pPD c 
which particle interaction cross sections can be meas- 
ured with an accuracy of several.percent. However, | dh 
the realization of this kind of accuracy is possible only 
if one can measure cross sections and energy of the 
proton beam at the same time. This situation is a 
result of the rapid increase in these cross sections as 
the energy is increased, especially in regions close to 
the threshold of the reactions being investigated. For 
example, the cross section for the formation of mesons 
in nucleon collisions at proton energies of 650 Mev 
varies at the rate of 0.7% / Mev, at 350 Mev this quan- 
tity varies at the rate of 3% / Mev. In view of this 
rapid increase in cross section at low proton energies 
one must know not only the mean energy but also the 
energy spectrum of the beam. It is especially desirable 
to know the latter in cases in which the energy of the 
proton beam is reduced by an absorber, a process which 
results in a considerable increase in the energy spread 
of the beam. 

In the present work we have investigated the 
energy spectrum of the proton beam at several values 
of the proton energy in the range 150-670 Mev. In 
addition, we have measured the mean energy of the 
protons for various operating conditions of the machine. 

















Proton Energy Spectrum 





In measuring the proton spectrum use is made of 
a magnet with pole pieces 100 cm in diameter which Fig. 1. Diagram of the experiment: P) proton 
provides a field of 16 koe in the gap. The magnet beam; Po) polyethylene absorber; Kg-Ks) collima- 
is located in the path of the external proton beam tors; PD) proton detector; C) automatic recording 
(Fig. 1) and can deflect the latter through angles up to potentiometer; B) unit for stretching wire. 





spectrometer the protons are recorded by a detector 
which consists of a thin (0.1-0.2 cm) plastic scintillator 
and an FEU-19M photomultiplier which is located 
behind it. The pulses from the multiplier are integrated 
by an RC network and the current (10° — 10 -Tamp) is 
measured by an automatic recording potentiometer 
(EPPV-51). By means of a synchronous motor the 
detector can be moved perpendicularly to the proton 
beam. Under these conditions the strip of the recording 
potentiometer registers the change in current density 

J with variation of the x coordinate (Fig. 1). The 
recording potentiometer also records coordinate markers 
so that the value of x can be determined with an 
accuracy of 0.02 cm. A typical J (x) curve is shown 

in Fig. 2. Only one detector such as that described 
above is used in the apparatus; this detector records 

the variation in the intensity of the proton beam with 
time. The time required for the measurement of a J 
(x) curve is 5 min. In the course of this time period 
the change in the magnetic field produced by the de- 
flection magnetis kept to less than 0.05%, 

As is apparent from Fig. 2, the function J (x) is 
described by a Gaussian curve. The dispersion of this 
curve o is determined by several factors (each of which 
is characterized by a particular component of the 
dispersion 0;): the energy spread in the beam cf, 
scattering of protons from the beam in air Og, the 
widths of the collimator and the scintillator O,, etc. 
The magnitudes of these dispersion factors vary in 
different ways as the distance Z is increased (cf. Fig. 1) 
so that the ratio 0 ,,/o passes through a broad maximum 
in the region 3< 7/R< 10. The quantity 7/R has been 
made equal to eight, providing a spectrometer resolu- 
tion of dE/Edx © 1% cm. In this case the ratio o,,/o, 
is 0.5, The quantity o, is reduced by placing in the 
path of the beam a thin-walled polyethylene tube filled 
with helium which makes it possible to reduce the 
scattering by a large factor. 

The proton energy spectrum is found by comparing 
the J (x) curves measured with the magnetic field 
switched on[J}; (x)] and with the field switched off 
[Jo(x)]. In order to measure the J,(x) curve the detector 
is placed in the path of the undeflected beam at point 
A (cf. Fig. 1). Under these conditions the square of 
the dispersion a, is determined by the difference 
O44 —(koy)*, where 0}, and Gp are the dispersion in 
the functions J,, (x) and Jo (x), respectively; the factor 
k characterizes the focusing effect of the magnetic 
field and is given approximately by 1-9 7/2(@ is 
the angle through which the beam is deflected in the 
spectrometer). For the values of the angle chosen 6 ~ 
16° this factor is close to unity: k = 0.96. In experi- 
ments carried out without helium the ratio 0 ¢/dp is 
0.5. The mean energy spread of the beam in this 
case is determined with an accuracy of approximately 
0.6 Mev but the relative shape of the energy spectrum 
cannot be ascertained. The use of helium makes it 
possible to increase the ratio Of /dg by a factor of 
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2.5 so that it then becomes possible to determine both 
the dispersion of, and the shape of the spectrum fairly 
accurately. The function F (x) which describes the 
spectrum can be found by solving the integral equation. 


Ju(a)=\Jy(k[e—E) F(a. 
0 


The energy spectrum of the extracted proton beam is 
found to be symmetric, It is described by a Gaussian 
curve: 


@ (E) > exp {—(E — E)*/2A}}. (2) 


At a proton energy E = 665 Mev the magnitude of 
Ax is found to be 3.0 + 0.6 Mev (measured without 
helium) and 2.7 + 0.3 Mev (measured with helium). 
The mean weighted value is 


Ar=(2,8 + 0,3) Mev. 


The spectrometer is calibrated (i.e., the quantities 
E and dE/dx are determined) by means of current- 
carrying wires (cf. below). 

The dispersion found in the present work Ag, differs 
from the estimate A, ~ 4 Mev given in [2] based on 
an ionization-chamber method [3]. The reason for this 
discrepancy lies in the fact that Af quoted in [3] is 
too high since the authors have not taken account of 
proton scattering. 

The method described here has been used to 
measure proton spectra at lower energies in which case 
the beam is stopped by passage through a polyethylene 
absorber (cf. Fig. 1). As the proton energy is reduced 
the energy spread of the beam increases because of the 
increased ionization losses and "straggling," due to 
fluctuations in the ionization losses. Hence the relative 
dispersion Ay /E increases rapidly as the protons are 

J(z) 
1,0}- 


0,8 
06 


04 


! 
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x,cm 
Fig. 2. The measured function J (x). The circles indicate 
values for a Gaussian function with standard deviation o. 





slowed down and this facilitates the measurement pro- 
cedure considerably. For example, at E = 185 Mev 
helium need not be used and wide collimators (0.7 cm) 
can be substituted since even under these conditions 

0, /op = 2.5. The functions J}; (x) and Jy (x), measured 
at low proton energies are found to be described by 
Gaussian curves, It follows that at low proton energies 
the proton spectrum © (E) is also given by an expression 
of the Guassian form as in (2). 

The measured dispersion Af is shown in Fig. 3. 
The results of the experiment are in agreement with 
the theoretical dependence plotted taking account of 
the increase in ionization losses and straggling [4]. 


Mean Proton Energy 





The most widely used technique for determining 
the mean energy of an external beam is the measurement 
of proton range by means of an ionization chamber 
[3]. However, as noted [3] this method always leads 
to a systematic uncertainty (2-3 Mev and E ~ 650 Mev) 
because of the fact that the proton scattering is not 
considered in the analysis of the Bragg curve. A proper 
calculation, which takes scattering into account, avoids 
this error; however, in this case, the ionization chamber 
technique does not provide an accuracy much better 
than 5 Mev at E ~ 650 Mev, because the ionization 
potential which appears in the range-energy relation 
is not known with the required accuracy. 

By measuring the emission angle of the Cerenkov 
radiation [2, 3]it is possible to determine the proton 
velocity directly. This method is most exact at low 
proton energies; as the proton energy increases the 
accuracy decreases (at E~650 Mev the accuracy is 2 
Mev). The optical system which is used to measure 
the Cerenkov angle is complicated and this requirement 
tends to limit the possibilities of practical application 
of this technique. 

In order to measure the mean proton energy with 
an accuracy better than 1 Mev we have used the current- 
carrying wire technique [5] (a diagram of the measure- 
ment system is shown in Fig. 1). A thin flexible wire 
filament is fastened at one end at the center of the 
detector (in front of the scintillator) and passes through 
the magnet and collimators, reproducing the path of 
the proton beam. At point B there is a device by means 
of which the wire is stretched by a weight T. A current 
i is passed through the wire. If T, i, and the mean 
momentum of protons in the beam P satisfy the relation 


P =cgT /i, (3) 


where c is the velocity of light and g is the acceleration 
of gravity, * the projections of the wire and the beam 
on the horizontal plane coincide. If the detector is 
located in such a way that the scintillator is placed at 
the maximum of the J (x) curve, and the wire passes 


through the center of the collimator slit of collimator 
Ko, and the current i is measured at this position of the 
wire, we can find P. 

Equation (3) applies if the wire is absolutely 
flexible and weightless. The elastic properties of the 
wire, however, give rise to a moment which opposes 
the electromagnetic forces. As a consequence, the 
value of the momentum computed from (3) differs 
(for a copper wire) from the true value by an amount 
5: 10° pd*/T (d is the diameter of the wire in cen- 
timeters). In addition to this correction it is necessary 
to take account of the sagging of the wire due to its 
own weight, which is given by the quantity d*J?/T 
( U9 is the length of the wire). The wire actually 
passes below the beam trajectory and this effect is 
important if the magnetic field is not uniform. The 
appropriate correction can be determined experimen- 
tally by moving the wire in the vertical direction. In 
the present case this correction is less than 0.05%. From 
the correction for the sagging of the wire under the 
effect of its own weight it follows that the length of 
the wire 2) must be made as small as possible. However, 
a reduction in the length of the wire leads to a reduction 
in the accuracy of the measurement since 


Ay/AP ~ L(l,—L—R)sin8/Pl, (4) 


(cf. Fig. 1). The minimum length 19 is determined by 
the accuracy Ay; jn, With which the wire can be 
located at the center of the slit of collimator Ky. In 
the present case L = 9 m, J) = 16 m, and Aypjp = 
=0.03cm. The wire is located at the center of the 
slit automatically by means of a contact device. The 
error AP/AP, given by (4) is 0.03%, As the diameter 
of the wire is reduced the size of the corrections are 
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Fig. 3. The beam dispersion Ay at various proton energies 
E. The curve shows the theoretical dependence of dis- 
persion on energy. 


* The value of cg at the latitude of Dubna is 2.9426 
if T is measured in g, i in amp, and P in Mev/c. 
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also rapidly reduced. Thus, with measurements with 
T = 100 g and d = 0.01 cm the correction for the 
elasticity of the wire is 0.4%( for P = 1300 Mev/c), 
while with d = 0.013 cm the correction is 1%. This 
correction is reduced by a factor of 2 or 3 if the wire 
is heated before the measurement by passing a high 
current through it. 

In the present work the measurements are carried 
out with d = 0.01 cm, T = 100 g. The current i varies 
from 0.2 to 1.0 amp. The quantities T and i are 
measured with an accuracy better than 0.05% so that 
it is possible to determine the mean beam energy E 
with an error smaller than 1 Mev: 

E = (667.1 + 0.7) Mev 
(at the output of the accelerator chamber). The mass 
of the proton is taken as 938.21 Mev. 

The accuracy indicated above can be realized 
in a limited range of proton energy (E > 250 Mev). At 
smaller energies it is necessary to use a small weight, 
a high current,and a thick wire; all of these factors 
reduce the accuracy of the measurements. 


Beam Energy for Various Proton Acceleration Conditions 





The mean energy of the beam depends on the 
strength of the magnetic field in the accelerator, the 
nature of the radial oscillations, and the conditions 
under which the beam is extracted; thus this energy 
depends on a large number of parameters which charac- 
terize the operation of the accelerator. For this 
reason the energy of the beam is not a fixed quantity, 
but varies over wide limits because of uncontrolled 
changes in the number of the accelerator parameters. 

J, amp 
280 290 } 300 310 520 














7000 4100 } 7200 4300 ~~ ta00 
J,, amp 
Fig. 4. Variation in mean energy of the pro- 
ton beam with variation in the magnet current 
Ip (@) and the current I, in the winding (3 ), 


which corrects the position of the median plane 
of the magnetic field. The arrows indicate the 
values of I, and I, which obtain under ordinary 
operating conditions. 


These variations in the magnitude of E amount to 
several million electron volts: 


Date of the 
Measurements 


1/5/57 
8/3/57 
8/20/57 
5/3/58 
6/5/58 
1/29/58 
8/20/58 
10/12/58 


E, Mev 


671.0+1.5 
667.84 1.0 
667.140.7 
658.84 1.0 
669.34 1.0 
666.34 1.0 
665.31 1.0 
663.22 1.0 


The value of E can be varied over certain limits 
by choosing an appropriate mode of operation. The 
necessity for control of the beam energy arises in a 
number of experiments in which one requires that the 
energy of the proton beam remain fixed. As has been 
shown by measurements, the energy E is sensitive to 
the strength of the magnetic field in the accelerator 
(Fig. 4). The control of E over a small range can be 
achieved also by varying the conditions of beam 
extraction (Fig. 5). By proper choice of these variables 
it is possible to increase the beam energy to the value 


} pon — 683 Mev. 


Under these conditions the intensity of the extracted 
proton beam is reduced by three orders of magnitude. 
Variations in the mode of acceleration and conditions 
of extraction are accompanied by small changes in the 
dispersion A;. Thus, when the current ly is increased 
from 4100 to 4250 amp the dispersion Ar is reduced 

by 0.5:0.3 Mev, while a change in the position (r) of 
the regenerator in the device by means of which the 
beam is extracted [1] causes the dispersion A, to vary 
as follows: 


for r=281.35cm Ag=3.7+0,5 Mev; 
for r=280.35c" Ag=2.8+0,3 Mev; 
for r=278.85cm Ag—1.8+0,5 Mev. 
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Fig. 5. Variation in the mean energy of the 

proton beam with variation of the position 

of the regenerator (r). 





As we have indicated earlier, measurements of the 
beam energy are especially important if an absorber is 
used to reduce proton energy (this method is used in 
many experiments). It is necessary to determine the 
energy of the degraded beam (this is inconvenient since 
the thickness of the absorber usually varies several 
times in the course of the experiment) or to determine 
the energy of the undegraded beam; in the latter case 
the energy spread relation for the absorber is used to 


5. g/cm’ 


Fig. 6. Energy spread Ay = E- EF’ in a poly- 
ethylene absorber of thickness S. E and E’ 
are the mean energies of the proton beam 
before and after the absorber. The numbers 
on the curves indicate the corresponding 
values of E, 


determine the energy of a stop beam. This relation has 
been measured by us for polyethylene (CHg), (it is also 
applicable for paraffin) in the range 150-665 Mev 
(Fig. 6). We may note that the usual range-energy 
relation can not be used for determining the energy of 
the protons instead of the indicated relation because 

it leads to large errors, especially for low energies of 
the degraded beam. Thus, if the range-energy relation 
is known with an accuracy of 2 Mev at E= 650 Mev 
(0.4%), when the beam energy is reduced to 100 Mev 
the uncertainty in the energy determination increases 
up to 6 Mev (6%), i.e., an increase of a factor of 15. 

In conclusion we wish to thank T'ang Hsiao- wei 
and A. A. Tyapkin for discussion of the results of the 
present work. 

Note added in proof. The method described here 
has been used to measure the mean energy of the 
external beam in operation of the synchrocyclotron for 
accelerating heavy particles: deuterons and 0 -particles. 
It is found that for deuterons extracted from the accel- 
eration chamber E = (405.3 + 0.5) Mev while for a- 
particles E = (811.34 1.0) Mev (for mass values of 
my = 1875.5 and mq = 3727.2Mev, respectively). The 
energy spectra of the beams are characterized by the 
following dispersions: for deuterons Ap=(1.720.5) Mev 
and for a-particles Ap = (3.5 41.5) Mev. 
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PHASE DIAGRAM OF THE ZIRCONIUM-RHENIUM SYSTEM 


E. M. Savitskii, M. A. Tylkina, and I. A. Tsyganova 


Translated from Atomnaya Energiya, Vol. 7, No. 3, pp. 231-235, 


September, 1959 
Original article submitted April 16, 1959 


The phase diagram of the zirconium — rhenium system was constructed by x-ray and microscopic methods and 
by measurements of melting points, hardness and microhardness of the phases, 

The region of solid solutions of rhenium in o-zirconium at 800° C is ~0.5 weight %, and at the temperature 
of the eutectoid transition it increases to 2-3 weight %, At 1600° C, 14.68 weight % of rhenium dissolves in 
B -zirconium and at the eutectoid transition point at 550-600° C, 8 weight % of rhenium. The 6 -phase is stabi- 
lized in alloys containing more than 4 weight % of rhenium. An eutectic is formed at 1600° C and 25 weight % 
of rhenium. The presence of a metastable w-phase in melts rich in zirconium was established. 

The solubility of zirconium in rhenium at 2500° C is less than 2 weight %, 

Three chemical compounds are formed in the system by peritectic reactions: 

1) At 2500° C, ZrsRey of the &-Mn type with a body-centered cubic lattice (a = 9.6-9.7 kX) and a micro- 


hardness of 1000 kg/ mm’; 


2) at 2450° C, ZrRe, with a hexagonal, close-packed lattice (a = 5.21-5.25 A; c= 8.5-8.56 A; c/a = 1.63) 


and a microhardness of 1200 kg/mm’, 


3) at 1900° C, Zr,Re of the o-phase type with a tetragonal lattice (a = 10.12 A; c = 5.42 A; c/a = 0.535) 


and a microhardness of 700-800 kg/mm’, 


Zirconium and rhenium are transition metals 
with an incomplete d-shell. The table gives some 
comparative data characterizing these metals. 

There is hardly any literature data on the zirconium- 
rhenium system. Only the chemical compound ZrRe, 
(Laves phase) is known and this has a lattice of the 
MgZn, type with the parameters a=5.251 kX, c=8.576 
kX, c/a=1.633 [1]. In [2] the authors of this article 
confirmed that the compound ZrReg is formed and 
established the formation of the compound Z1,Re, 
melting at 1900° C. 

As starting materials we used iodide zirconium 
and rhenium powder with 99.8 weight % purity, which 
were introduced into the charge in the form of briquets, 
first sintered at 1400° C in vacuum. 

The melts (25 samples of 30 g each) were pre- 
pared by arc melting with unconsumed tungsten elec- 
trode on a water-cooled copper bottom in a helium 
atmosphere at 200 mm Hg. The melts were fused two 
to three times for a uniform distribution of the com- 
ponents. The compositions of the alloys, determined 
from the charge and by chemical analysis, differed 
insignificantly. 

The microstructural analysis was carried out on 
sections, prepared in the usual way. An aqueous solu- 
tion of hydrofluoric and nitric acids in various propor- 
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tions, depending on the composition of the alloy, was 
used for etching. As the rhenium content of tne alloy 
increased, the nitric acid content of the solution was 
increased. X-ray patterns of powders or cylindrical 
samples, 0.6-0.8 mm in diameter, were obtained in 

a Debye camera with unfiltered chromium radiation.* 
In addition, we measured the Vickers hardness of forged 
and heat treated alloys and the microhardness on a 

PMT -3 apparatus for various phase components. 

The zirconium corner of the phase diagram was 
investigated on forged samples. The alloys were forged 
into bars with a 7x 7 mm cross section at 800-1000° C. 
Phase analysis of these alloys was carried out with 
samples as cast and also tempered after forging. Iso- 
thermal annealing with quenching in water was done 
in evacuated quartz ampoules at 1050 and 950° C for 
30 minutes, at 850, 800, and 700° C for 6 minutes, at 
500° C for 860 hours, and at 400° C for 2 months, 

The sections of the diagram from 45 to 100 weight 
% of rhenium were examined by means of cast samples 
as they shattered on forging. Phase analysis was carried 
out on cast alloys and alloys annealed to 1000° C for 


* The x-ray investigations of the zirconium— rhenium 
system were carried out by P. I. Kripyakevich and E, I. 
Gladyshevskiiat Leningrad State University. 





Comparative Data for Zirconium and Rhenium 





Characteristics 


Zirconium 


Rhenium 





Atomic radius, A: 
Q-zirconium, seer eee eeeeesn 
B - zirconium, eee er ew eeeesn 


Lattice parameters, A.......6.2.. 


Polymorphic transition point, °C... 
Melting point, °C. .....cccee5 
Specific gravity, g/cm®,....... 


WACOITIIG “VOMIING, 0 6 6.6 0:0 0.0 6 @ wee 


14,2 


ai 
1.61 
a-hexagonal, close-packed 
(a = 3,223, c = 5,123, 
c/a = 1,59) 
B -cubic, body-centered (a = 
= 3,12) 
862° + 5 
1860 
6.4 


8.9 
1,37 
a-hexagonal, close-packed 


(a = 2.755, c = 4,449, 
c/a = 1.614) 


Difference between atomic radius 
and that of rhenium, %: 
CE“SMCOMIUM «ccc ec ce te 
S-mirOohium . vs veceeece 

















560 hr with subsequent quenching in water. The phase 
composition of cast alloys is arbitrarily plotted on the 
diagram on the isotherm at ~1500° C. It can be con- 
sidered that the structure of cast alloys, rapidly cooled 
on the copper bottom of the arc’ furnace, corresponds 


approximately to alloys quenched from high tempera- 
ture. 


toid decomposition of the 6 -solid solution occurs at 
550° C by the reaction 8 -solid solution = c-solid 
solution + ZrRe. 


at. % 
60 =80 
T PR 

This solidus curve was determined in vacuum on | /1. 
samples 8 x 5 x 5 mm in an apparatus for measuring L*Z%, Rey, /” M 
the melting point of refractory metals [2]. The melt- 2500" 
ing point was determined with an optical pyrometer Sr% \ 
("drop method") as the point when a hole drilled in y” 
the sample filled with metal. The pyrometer was i. eZrRe 
calibrated with the aid of pure metals under the same i - 
conditions. The melting points were determined as the 
average of two to three measurements. 

The phase diagram (Fig. 1) was constructed from 
the data of x-ray and metallographic analysis, melting 
point determinations, and hardness measurements. 

The compounds Z1,Re, ZrRe,, and ZrsRey4, which 
were formed by peritectic reactions, were detected , 
in the system. Zirconium gives limited a- and B - | 
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solid solutions with rhenium. A eutectic is formed 
between 8 -zirconium and Zr,Re. Eutectoid decompo- 
sition of the 6 -solid solution occurs in the solid state 
with the formation of an o-solid solution and the 
compound Zr,Re. 

The region of solid solutions in 8 -zirconium ex- 
tends to 15 weight % of rhenium at 1600° C, but is 
reduced to 8 weight % at the eutectoid transition 
point. The 6 -phase is stabilized with a rhenium 
content of more than 4 weight %, The polymorphic 
transition point of zirconium falls from 862 to 550 
600° C with an increase in the rhenium content. Eutec- 
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Fig. 1. Phase diagram of the zirconium — rhenium 
system: A) single-phase structure; O) two-phase 
structure; @) solidus. 





The solubility of rhenium in «-zirconium is low: 
0.5 weight % at 700-800° C. The solubility increases 
somewhat with a fall in temperature and is less than 
3 weight % at 400°C. At 1600° C and~25 weight % 
of rhenium there is the eutectic transition liquid = 
B -solid solution + Zr1,Re. 

After quenching from 850-700° C, the alloy with 
0.79 weight % of rhenium (Fig. 2) has a two-phase 
structure ( a+ 8). 

X-ray investigation of zirconium-rich alloys, 
quenched from 1000° C, made it possible to detect 
the presence of a metastable w-phase with a rhenium 
content of “5.6 at. %, This phase has a complex, 
noncubic lattice. Analogous phases have been found in 
systems of titanium with a number of elements and are 
formed during the aging of titanium-rich alloys [3]. 


Fig. 2. Alloy with 0.79 weight % of 
rhenium, x 200. 
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Fig. 3. Change in hardness of zirconium-rhenium 
alloys during heat-treatment; x) annealing at 500° C, 
860 hr; e) quenching from 700° C; O) quenching 

from 800° C; @) quenching from 850° C; A) quench- 
ing from 950° C; 4) quenching from 1050° C. 
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Fig. 4. Microhardness of phases in the zirconium - 
rhenium system. 


Fig. 5. Alloy with 98 weight % of rhenium, 
x 340. 


Alloying of zirconium with rhenium leads to a con- 
siderable increase in hardness (Fig. 3). Thus, for 
samples annealed at 500° C for 860 hr, the hardness 
of pure zirconium (136 kg/ mm?) increased to 350 kg/ 
mm? for an alloy with 14.68 weight % of rhenium and 
to 420 kg/mm? for an alloy with 38.4 weight % of 
rhenium. A similar change in hardness occurred with 
various quenching conditions, when the numerical 
values of the hardness increased correspondingly with 
a rise in quenching temperature. In addition, the micro- 
hardness of separate phases was determined (Fig. 4). 

The solubility of zirconium in rhenium at the 
melting point was less than 2 weight %. An alloy with 
98 weight % of rhenium already had two phases (Fig. 5). 
and consisted of a y -solid solution of zirconium in 





thenium + ZrsReg,. A hypothetical line for the solubility 
of zirconium is plotted on the diagram. 

The compound ZrsRe,, with a rhenium content of 
90-92 weight % is formed by the peritectic reaction: 
liquid + y -solid solution of Zr in Re= ZrsReg, at 
2500° C and crystallizes in the form of dendrites with 
somewhat expanded faces. The lines of this phase 
on x-ray patterns indicated a cubic lattice (a=9.60— 
—9.70 kX) [4]. 

The lattice constant of the ZrsReg,4 is close to the 
constant of x -phases, i.e., compounds of the &-Mn 
type or a superlattice related to it of the TisReg, type, 
observed in systems of transition metals with rhenium 
(for example, for Nbg7Reg, a= 9.67 kX [5], and for 
TisReg, a = 9.59 kX [6, 7]). 

Accurate determination of the lattice constant 
of the compound showed that for the alloy with 80.6 
at. % of rhenium a = 9.693 kX and for the alloy with 
82 at.% of rhenium a = 9.626 kX. Consequently the 
compound has a variable composition. 


Fig. 6. Alloy with 83 weight % of rhenium, 
x 200. 


The microhardness of the ZrsRey phase was ~1000 
kg/mm*. More detailed data on this compound was 
published in [4]. The compound ZrRe, (80 weight % of 
rhenium) is formed at 2450° C by the peritectic reaction: 
liquid + ZrsRey, “ZrRez. The microhardness of the 
compound was 1200 kg/mm?, 

Between the compounds ZrRe and ZrsRey lies a 
two-phase region (Fig. 6). 


Fig. 7. Alloy with 55 weight % of rhenium, 
x 200. 


The compound ZrRe, analogous to o -phases, is 
formed with a rhenium content of 50 weight % at 
1900° C by the peritectic reaction: liquid+ ZrRe, = 
=*ZrRe. The lattice is tetragonal with the parameters: 
a=10.12 A; c=5.42 A; c/a=0.535. The micro- 
hardness is 700-800 kg/ mm and the region of homoge- 
nity of the compound is very narrow (see Fig. 4). 
Similar phases are encountered in many systems with 
rhenium, for example, rhenium — molybdenum [8]. 


Alloys with 55-88 weight % of rhenium have a 
two-phase structure, ZRe + ZrRe, (Fig. 7). 
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EXTRACTION PROPERTIES OF DIISOAMYL 


METHY LPHOSPHINITE 


V. B. Shevchenko, I. G. Slepchenko, V. S. Shmidt, 


and E. A. Nenarokomov 


Translated from Atomnaya Energiya,Vol. 7, No. 3, pp. 236-243, 


September, 1959 
Original article submitted December 11, 1958 


The properties of one of the new, efficient extractants for uranium, diisoamyl methylphosphinate (DAMP), are 
described, It was shown that uranyl nitrate is extracted from nitric acid solutions by this extractant with consi- 
derably higher distribution coefficients than tributyl phosphate (TBP). The extraction of uranyl nitrate and HNO, 
with DAMP solutions in hydrogenated kerosene was studied, It was shown that uranyl nitrate is extracted in the 
form of the complex [UO,{NO3)2 (DAMP)2], whose stability constant equals 2540 + 200, and HNO; is extracted 

in the form of the compound HNOs DAMP, whose stability constant is 0.30 + 0.03. 


In recent years extractants with both high selecti- 
vity and considerable extracting capacity have been 
investigated [ 1, 2]. One of these extractants is TBP, 
which has been described in detail in the literature 
[1-14]. The coefficients of UO NOg), distribution 
between nitric acid solutions and TBP are higher than 
for the diethyl and dibutyl esters, methyl isobutyl 
ketone, and other extractants which were previously 
used for the extraction of uranyl nitrate [15-19]. 

In the present work we give the results of 
studying the extraction properties of DAMP 


ass: See, | E 
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and solutions of it in hydrogenated kerosene. 

The physical properties of DAMP are as follows: 
specific gravity 0.95 g/cc, boiling point at 7 mm Hg 
130° C, viscosity 4.48 cp, flash point 130° C, and 
refractive index n”=1,4284, The kerosene used as a 
diluent was a mixture of saturated hydrocarbons with a 
boiling point of 170-240° C and a specific gravity of 
0.76 g/cc. Solutions of DAMP in kerosene were pre - 
pared by mixing appropriate volumes of these substances? 

For experiments with weighable amounts, the 
uranyl nitrate was purified by extraction with 20% 

TBP and subsequent reextraction with water. The 
isotope U*5 was used in experiments with microamounts 
of uranium. The U™* preparation contained an appre- 
ciable amount of Th and therefore the thorium 

was first separated from uranium by extraction of the 
latter from 6 N HCl with a 20% solution of TBP. It is 
known that thorium and protactinium are not extracted 
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under these conditions. The reextract containing 
us was repeatedly evaporated,almost to dryness with 
excess HNOsg,and the final residue diluted to 5 ml with 
water. This gave a solution with a specific activity 
of ~2 + 10° dis/ml~- min. This solution was used as 
a stock solution of U"* and ~0.01 of its volume was 
introduced into an aqueous phase of the required com- 
position. 

Equal volumes of the phases were used for extrac- 
tion. In work with weighable amounts of uranium the 
volume of each phase equaled 10 ml with microamounts 
of 2.0 ml. 

The phases were mixed with a mechanical shaker 
for 30 min ; the settling time was also 30 min. In 
work with weighable amounts of uranium, the distri- 
bution coefficients were determined as the ratio of 
the analytical concentrations of uranium in the organic 
and aqueous phases, and in work with microamounts 
of uranium, as the ratio of the a-activities of equal 
volumes of the organic and aqueous phases. Here the 
a-activity of the organic phase was calculated from 
the difference in the activity of the aqueous phase 
before and after extraction. In order to avoid a change 
in the phase volumes during the extraction of uanium 
in work with microamounts of uranium, before intro- 
duction of the U*’, the aqueous and organic solutions 
were equilibrated in a different separating funnel, from 
which the required volumes of solutions were withdrawn. 
The &-activity was measured on an o-scintillation 
counter. The experiments were at 19 4 1° C, 

The HNOsg distribution coefficients were determined 
by titration of the aqueous phase and the aqueous 
* The DAMP content of the extractant is given in 
percent by volume in every case. 





reextract from the organic phase with 0.1 N NaOH 
solution. 
Extraction of Nitric Acid 

Curves of the equilibrium distribution of HNO, 
between aqueous solutions and solutions of DAMP in 
kerosene are given in Fig. 1. The relations between 
the HNOg distribution coefficients and the HNO  con- 
centrations in the aqueous phase at various DAMP con- 
centrations in the extractant are given in Fig. 2. 
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HNOs concentration in organic 
phase, mole/ liter 
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HNOs concentration in aqueous phase, 
mole/ liter 
Fig. 1. Distribution of HNO, between aqueous solutions 
and extractant. Concentration of DAMP (vol. %): 1) 5; 
2) 10; 8) 20; 4) 40; 5) 100. 


HNO; distribution coefficient 


0 re eee 8 ee 8 
HNO; concentration in aqueous phase, 
mole/ liter 
Fig. 2. Distribution coefficients of HNO, between the 
aqueous phase and the extractant. Concentration of 
DAMP (vol. %): 1) 5; 2) 10; 3) 20; 4) 40; 5) 100. 


Figures 1 and 2 show that the distribution coeffi- 


cient a increases with an increase in the DAMP 
HNO 


concentration in the extractant. Here, as for TBP 
(11), with an increase in the acidity from 0 to 0.7 —- 


-1,0 N there was an increase in a » connected 
HNO 


with suppression of HNOg dissociation; with a further 
increase in acidity there was a fall in ONO. 2PPa- 
rently caused by saturation of the extractant. 


It can be conjectured that the extraction of HNO, 
by DAMP solutions, as with TBP [3, 5], is connected 
with the formation of a readily extractable complex 
by the reaction 


H*+NO;+ DAMP = HNO, DAMP. (4) 
Calculation of the equilibrium constant of this 
reaction 


K = __LINOs DAMP | 
[H*}[NO3}{ DAMP | ’ (2) 





based on experimental data on the equilibrium distri- 
bution of HNOs (Tables 1 and 2), showed that the value 
of K remained constant at HNOs concentrations in the 
aqueous phase of up to 3 mole /liter and equaled 0.30 
+ 0.03. 

The experimental values of a and those 
calculated by means of this value of * the constant 
are given in Tables 1 and 2. The calculation procedure 
in principle was the same as that described for the 
extraction of HNOs by tributyl phosphate [3]. The value 
of HINO; was Calculated by the formula 


G@HNo; = K[H*] [D]=KA, (D,—N,): - 
( 


TABLE 1. Distribution Coefficients of HNO, between 
Aqueous Solutions and a 20% Solution of DAMP in 
Kerosene (Dy = 0.80 M; K = 0.30). 





Equilibrium concentration of 
HNO;, mole/ liter 





in aqueous in organic 
phase pase 























TABLE 2. Distribution Coefficients of HNO, between 
Aqueous Solutions and a 10% Solution of DAMP in 
Kerosene (Dp = 40.0 M; K = 0.30). 





Equilibrium concentration of 
HNO ;, mole/ liter 


in aqueous in org anic 
phase phase 








0,116 
0,167 
0,236 
0,318 























where Dy is the DAMP concentration in the extractant 
before extraction, D is the unbound DAMP concentra- 
tion in the extractant after extraction, N, is the equi- 
librium concentration of HNOgs in the aqueous phase, 
and N, that in the organic phase. 

As the data in Tables 1 and 2 show, the calculated 
and experimental values of the HNOg distribution coeffi- 
cients are quite close, indicating the accuracy of the 
hypothesis on the mechanism of HNOs extraction by 
DAMP solutions (1). 

It is interesting to note that the equilibrium constant 
of the reaction H* + NOs + DAMP HNO, DAMP 
derived for the extraction of DAMP differs comparatively 
little from the corresponding value for the extraction 
of HNOs by tributyl phosphate (0.22 + 0.02), given in 
(3, 5]. This indicates that the difference in the stability 
of the compounds HNO, DAMP and HNOs TBP is com- 
paratively small. 

In Fig. 3, curve 1 shows the equilibrium distribu- 
tion of HNO, between an aqueous solution and a 20% 
solution of DAMP containing 60 g /liter of uranium. 

As can be seen, when the DAMP solution was saturated 
with uranium, the HNOsg distribution coefficient fell. 
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Fig. 3. Distribution of HNO, between an aqueous 

solution and a 20% solution of DAMP: 1) the organ- 

ic phase contained 60 g/liter of uranium, 2) the 

organic phase contained no uranium. 


Extraction of Uranyl Nitrate 





Experiments with weighable amounts of UO, (NOsg)s. 
Curves of the equilibrium distribution of UO, (NOs), 
between 0.05 N HNOs and solutions of DAMP in kero- 
sene are presented in Fig. 4. The relations between the 
uranium distribution coefficients and the UO, (NOs), 
concentrations in the aqueous phase for different DAMP 
concentrations in the extractant are given in Fig. 5. 
The same figures also show data for TBP [11] for com- 
parison. Figures 4 and 5 show that when DAMP solu- 
tions are used as the extractant, the distribution coeffi- 
cients of uranium between the aqueous and organic 
phases considerably exceed the corresponding values 
for TBP solutions of the same concentration. This 
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Fig. 4. Distribution of UO, (NOs), between 0.05 N 
HNOs solution and extractant. 
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Fig. 5. Distribution coefficients of UO, (NOs), 
between 0.05 N HNOs solution and extractant. 


TABLE 3. Composition of DAMP Solutions, Saturated 
with Uranyl Nitrate. 





DAMP con- funy 
tent of ex- 
tractant 


unt of uranium diss0l Ratio DAMP 


¢ in extractant sample 
‘at satura on( the taal [U] in satu- 
rated organid 


t 
Poa Tee ex solution 


|mole/ 
g mole 


vol. “Bititer 








40 48,6 
80 94,0 
6 
0 


0,202 
0,395 
0,790 
1,99 


186,2 
473,8 














* Due to the increase in volume of DAMP when 
UO,(NO3), was dissolved in it, the volume of the 
organic solution containing the given amount of 
uranium exceeded 1 liter, 











TABLE 4. Distribution Coefficients of UO,(NOs), bet- 
ween 0.05 N HNOs and a 20% Solution of DAMP in 
Kerosene (Dp = 0.80 M; Ky = 2540). 





Equilibrium conc. of uranium _ 





in aqueous 
phase 


in organic 
phase 





mole/ | g/ 
liter liter 


mole/ 
liter 





0,100 
0,123 
0,210 
0,315 
0,508 


0,224 
0,260 
0,320 
0,353 
0,375 


























difference is especially marked for low uranium concen- 
trations in the aqueous phase. With high uranium con- 
centrations in the aqueous phase, saturation of the 
organic phase with uranium occurred. As Fig. 4 shows, 
the uranium content of the organic phase remained 
practically constant at uranium concentrations in the 
aqueous phase of = 150 g/liter. 

Data characterizing the composition of the organic 
phase at saturation are presented in Table 3. 

The constancy of the ratio [DAMP] : [U] in the 
organic phase at different dilutions, saturated with 
uranium, shows that uranyl nitrate is apparently extrac- 
ted in the form of the complex UO, (NOs)  2DAMP, 
formed by the reaction: 


UO2+ + 2NO>+42 DAMP = 
= UO,(NO,),2 DAMP. (4) 


Calculation of the equilibrium constant of this 
reaction from experimental data 


___ {UO,(NO,),2 DAMP } 
~~ {003+} [NO5]?} DAMP 143 (5) 





Ky 


for a 20% solution of DAMP in kerosene shows that the 
value of K,, remains constant with changes in the 
uranium concentration in the aqueous phase over a wide 
range and equals 2540 + 200. Table 4 gives experi- 
mental values of the distribution coefficients of UO, 
(NO ), and values calculated from the value of Kj. The 
calculation was catried out by means of the formula 
given for extraction equilibria with TBP [12] with the 
DAMP concentrations substituted in the formula instead 
of the TBP concentrations: 


ay = Ky (22)? Dy’, (6) 


where K,, = 2540, x is the equilibrium concentration 
of uranium in the aqueous phase, D is the equilibrium 
concentration of free DAMP in the organic phase ( D = 
= Dy — 2 y, where Dy is the initial concentration of 


DAMP in the extractant and yis the equilibrium con- 
centration of uranium in the organic phase), and y are 
the activity coefficients of uranyl nitrate, whose values 
are taken from [20]. 

Formula (6) was deduced from (5) on the basis of 
the extraction mechanism expressed by (4). Here it 
was assumed that the whole of the uranium in the 
organic phase was in the form UO, (NOs), 2DAMP, and in 
the aqueous phase as UO}*. 

The satisfactory agreement of the calculated and 
experimental values of indicates the accuracy of the 
hypothesis made on the mechanism of uranyl nitrate 
extraction by DAMP solutions. 

As can be seen from the data presented, the stabi- 
lity constant of the complex UO, (NOs), with DAMP is 
extremely high and is approximately 39 times greater 
than the corresponding constant for TBP ( ~66) [12]. 
Such a high value of Ky indicates that NOs ions and 
DAMP molecules are apparently in the internal coor- 
dination sphere of the complex formed; the composition 
of the latter is probably expressed by the formula 
[UO, (NOs), (DAMP),]. 

A rise in the HNO, content of the aqueous phase 
over the concentration range studied produced an appre- 
ciable increase in the distribution coefficient of 
uranium. Curves of the equilibrium distribution of UO, 
(NOs), between aqueous solutions of various acidities 
and 20 and 40% solutions of DAMP are shown in Figs. 6 
and 7. As has been shown [13], with simultaneous 
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Fig. 6. Dependence of the UO,(NOs), distribution 
between an aqueous solution and extractant con- 
taining 20% of DAMP on the HNO; concentration 
in the aqueous phase: 1) 0.5 N HNOs; 2) 1.0N 
HNOs; 3) 2.0 N HNOg; 4) 3.0 N HNOs; 5) distri- 
bution of UO,(NO 3), between 0.05 N HNOsg and 
20% TBP [12]. 


extraction of UO, (NOs), and HNO; by tributyl phosphate, 
the uranium content of the organic phase may be cal- 
culated from the formula 


Bee ae V1+8FT,—1 
vod [ry 


4F (7) 
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Fig. 7. Dependence of the UO,( NOs), 
distribution between an aqueous solution 
and extractant containing 40% of DAMP 
on the HNOs concentration in the aqueous 
phase: 1) 0.05 N HNOg; 2) 1.0 N HNOg; 
3) 2.0 N HNOg; 4) 3.0 N HNOs; 5) distri- 
bution of UO,( NOs), between 0.05 N HNOs 
and 40% TBP [12]. 


where Ty is the TBP concentration in the extractant 
before extractions; 


F =fy/(1+ fu); i 
fu = Kuy’ (22+ N,)*; (8) 


fu =0,2N, (22 +N,). (9) 


Since the mechanisms of UO, (NOs), and HNO, 
extraction by DAMP and TBP solutions are analogous, 
the formulas presented should evidently also apply 
to extraction equilibria with DAMP. 

Replacing the coefficient 0.2 in(9) by 0.3 (the 
equilibrium constant of HNO extraction) and substituting 
the value Ky = 2540 in (8), we calculated the values 
of y and Oy. The experimental and calculated values 
of a,, for extraction of UO, (NOs), from 1.0 N HNO; 
with a 20% solution of DAMP are given in Table 5. 

Table 5 shows that at uranium concentrations in 
the aqueous phase of 0.10 — 0.16 mole / liter, the cal- 
culated values of o,, agree quite well with the experi- 
mental values, This also confirms the accuracy of the 
hypothesis on the extraction mechanism (4) and the 
accuracy of the value of Ky used (2540). However, 
calculations for uranium concentrations less than 0.1 
M could not be carried out as there are no appropriate 
data on the activity coefficients of uranyl nitrate. For 
UO, (NOs), concentrations in the aqueous phase greater 
than 0.16 M, the organic phase was almost saturated 
with uranium and calculation by (7) gave a large error. 


132 


TABLE 5. Distribution Coefficients of UO,( NOs), 
between 1.0 N HNO, and a 20% Solution of DAMP in 
Kerosene (Dp = 0.80 M; Kyy = 2540). 





Equilibrium conc, of uranium 


in aqueous 
phase 





in organic 
phase 





g/ mole/ g/ 
liter liter liter 


mole/ 
liter 





0,344 
0,355 
0,361 
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Fig. 8. Dependence of the coef- 
ficient of UO,( NOs), distribution 
between nitric acid solution and 
20% DAMP solution on the uranium 
content of the aqueous phase. 
Acidity of aqueous phase: 1) 0.05 N 
HNOs; 2) 2.0 N HNOsg. 


Figure 8 shows the dependence of the coefficients 
of UO, (NOs), distribution between 0.05 and 2.0 N 
aqueous HNOs solutions and a 20% DAMP solution on the 
uranium content of the aqueous phase. The figure 
shows that With 2.0 N HNOs the distribution coefficients 
are particularly great at low uranium concentrations 
in the aqueous phase (up to a value of ~5 g/liter). 
This shows the efficiency of using DAMP in cases where 


it is necessary to completely extract UO, (NOs), from 
the aqueous phase. 
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Fig. 9. Dependence of the 
coefficient of UO,(NOs), distri- 
bution between aqueous solu- 
tions and 20% DAMP solution 
on the content of salting-out 
agent in the original aqueous 
solution. 
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Fig. 10. Dependence of the distri- 
bution coefficient of microamounts 
of UO,(NOs), on the DAMP concen- 
tration in the extractant (in loga- 
rithmic coordinates). 


The introduction of salting-out agents into the 
solution considerably increased the uranium distribution 
coefficients (Fig. 9). Aluminum, iron, and magnesium 
nitrates were particularly efficient salting-out agents, 

Experiments with microamounts of UO, (NOg),. As 
an additional check on the accuracy of the UO, (NOs), 
extraction mechanism proposed, we measured the 


relation between the distribution coefficient of micro- 
amounts of UO, (NO 3), and the DAMP concentration 

in the extractant (Fig. 10). The relation lg a =f {lg 
[D]} was expressed by a straight line with a slope equal 
to 2, which would be expected for extraction of a comp- 
lex containing two molecules of extractant. A linear 
relation was observed only for DAMP concentrations 

= 10%, With a further rise in the DAMP concentration, 
the line deviated toward the abscissa axis. 

The data presented show that DAMP is a stronger 
extractant for uranium than TBP. 

The relations between the distribution coefficients 
and the composition of the aqueous and organic solutions 
studied indicate that the compound extracted is the 
complex UO, (NOs), (DAMP)p. 


SUMMARY 


1. A study was made of the dependence of the 
coefficients of HNOg and uranyl nitrate distribution 
between nitric acid solutions and kerosene solutions of 
DAMP on the DAMP content of the extractant and on 
the UO, (NOs), and HNOg contents of the aqueous phase. 
It was shown that under identical conditions DAMP 
extracts considerably more uranium than TBP, especially 
with low uranium concentrations in the aqueous phase. 

2. The rules of HNO; extraction by DAMP solutions 
may be explained on the basis of the following extraction 
mechanism: 


H*+NO;>-+ DAMP = HNO, DAMP 
where HNO; DAMP is the compound extracted completely 
into the organic phase. 

3. The rules of uranium extraction from nitric acid 
solutions by DAMP solutions may be explained satisfac- 
torily if it is assumed that extraction is determined by 
the reaction: 


U02++2NO;+2DAMP = 
= [UO,(NO,),( DAMP ),], 


where [UQ,( NO 3) DAMP) ] is the compound extracted 
completely into the organic phase. 

4. We determined the equilibrium constants of 
the reactions 


H*+NO>-+ DAMP = HNO, DAMP 
nu UO2++2NO>+42 DAMP = 
= [UO,(NO,),( DAMP ),]. 


The constant for the reaction of HNO; with 10 and 20% 
solutions of DAMP equals 0.804 0.03 . 

The equilibrium constant of the reaction of UO, 
(NOs)2 with 20% DAMP equals 2549 4 200. 
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This paper describes an experimental investigation of the accumulation factor in the region of small distances, 
of y radiation from Co™ and Cs” point-sources. The possibility of using literature data for a homogenous in- 
finite medium for practical calculations (under conditions of barrier geometry for a relatively small thickness 
of absorbing layers) is demonstrated. The possibility of equalizing the dosage field of a cylindrical irradiator 
by means of additional filters has been proven by experiment. Moreover, the authors propose to replace Co™ 


specimens with "> td 


One of the important possibilities for the use of 
radioactive isotopes is their application in y devices 
for the irradiation of different specimens for experi- 
mental and industrial purposes in the field of biology, 
agriculture, and radiation chemistry. 

In the majority of cases, the irradiation of compa- 
ratively large-size specimens with a density close to 
unity is required. In this, it is most convenient to use 
long-lived isotopes which emit high-energy y radia- 
tion. Co™ and Cs!*" are such isotopes. 

Gamma apparatuses with Co™ and, especially with 
Cs *" can be mobile, and they do not require a fre- 
quent replacement as sources. An important advantage 
of these isotopes, especially in the case they are used 
for scientific research experiments, is the fact that their 
y radiation is relatively monochromatic. 

The basic characteristics of irradiators are: 

1) the radiation dose rate P (r/min) at the center 
of the irradiator operating volume and the required 
degree of its regulation; 

2) the magnitude of the irradiator operating 
volume V (cm’); 

3) the maximum allowable degree of nonunifor- 
mity of the dosage field in the operating volume 
(6 =(P max/P min) 100%); 

4) efficiency of the y device (kg/hr); 

5) the over-allactivity Q(g-equiv Ra or C) of 
isotope sources. 

Irradiation characteristics are interdependent, and 
they are determined by the irradiator shape, the irra- 
diation method, and the object parameters. 

Depending on the purpose for which the device 
is used, certain characteristics are of greater or lesser 
importance. Thus, for instance, in radiobiological 
experiments on animals and plants, the uniformity of 
the dosage field in the operating volume and the 


specimens in irradiating objects with a thickness of up to 40 cm and a density = 1. 


possibility of changing the dosage rate in a wide range 
are especially important, but the irradiator efficiency 
and the amount of the isotope used are of no special 
importance. Conversely, the following characteristics 
determining the device efficiency are of special 
importance in industrial devices: the required amount 
of the isotope and the output; the dosage field unifor- 
mity is of lesser importance, and, in the majority of 
cases, the dosage rate does not have to be changed at 
all [1 and 2}. 

In irradiator design, the integral dose which a 
specimen of a given size and density is to receive over 
the entire period of irradiation in the device and the 
magnitudes of the first four characteristics are usually 
determined. The amount Q of the radioactive substance 
in the irradiator must be determined. The assigned 
requirements are satisfied by selecting the irradiator 
shape, its parameters, and the irradiation method. How- 
ever, because of the large number of different factors 
affecting the quantity Q, usually several variants of the 
device are designed. Therefore, the design method is 
of great importance in practice. Until now, no papers 
devoted to general problems in the design of powerful 
y devices for radiation treatment have been published. 
Only papers devoted to the development and calcula- 
tion of individual devices [ 3, 4 and 5] and to calcu- 
lations of the dosage field of extended sources of 
different shapes [6, 7, 8, 9, 10,and 11] are available. 
The majority of the theoretical solutions available 
in the latter papers are not suitable for practical cal- 
culations, where a number of specific requirements, 
which cannot be easily foreseen in general theoretical 
derivations, have to be taken into account. 

The purpose of the present paper is the following: 

1) experimental verification of the admissibility 
of certain approximate calculations by comparing the 
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calculations with measurements performed under actual 
operating conditions on a y device which is in use; 

2) determination of the most suitable irradiation 
method in using extended sources; 

3) determination of the fields of application of 
Co™ and Cs!*" for irradiation. 


Experimental Verification of Some Approximate Cal- 
culation Methods 

If a dosage field where an absorber is present is 
to be calculated, the multiple- scattered radiation in 
the given medium should be taken into account. The 
Fano theory [12] for the dose rate from a point-source 








in a homogenous medium gives the following expression: 


P =P, B,e-“F, 


where LU is the linear attenuation coefficient of a 
narrow beam in a given medium, F is the distance from 
the source to the measurement point, By is the accu- 
mulation factor, which takes into account the contri- 
bution of scattered radiation, and Pp is the dose rate 

in air, provided by the source at the measurement 
point. The accumulation factor By can be represented 
as the sum of two exponentials [9] of the following 
form 


B. = A,e—1¥0F 1 A.e—“avoF 
0 l + Age . 


For y radiation of Co™ and Cs®" in a medium equi- 
valent to water with respect to its absorption-scattering 
capacity, the expression determining the dose rate 
provided by a point-source assumes the following form: 


K. Q 
60. F2 
— 7,9¢e—9:967F\ r/min ( for 
K {2 
es —0,0753F __ 
P = : GO. , (19,5e 


— 18,5e—9.086F)+/min (for Cs'%’), (2) 


Fr (8, 9e-0,057F — (1) 


Co®), 


where Q is the point-source activity. The values of 
the coefficients in the expression for the accumulation 
factor are taken from [9]. 

For a linear source, for a point at the distance 
F from the middle of the source, we similarly obtain 


arctg 2 
i eee 8,9 \ e-0,057F -sec @ dy — 
bor” 


. (3) 


arctg ce 
arctg 


mae 7,9 ( ¢~9.067F «sec 9 dy | r/ min (for Co), 
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arcte 
2K oq r 
peel?) Sed »>~0,0753F- sec ? a 
or [ 19,5 5 ¢ dp 


arctg 4 
F 


— 18,5 ( e—0,086F -see ? dg |/ rnin 
0 
( for Cs1%7), (4) 


where L is the half-length of the linear source, Ky is 
the y constant, and q is the linear specific activity 
(g-equiv /cm or C/cm). 

If the absorbing medium fills only a part of the 
gap F, the degree of accuracy provided by the appro- 
ximate expressions (1)-(4) should be checked. 

Moreover, on the basis of an experimental check 
[13] which was performed with a Co™ point-source 
in water, it has been found that the Fano theory yields 
a somewhat higher value for the accumulation factor 
in the region of small distances, where 100 cm > F > 
45 cm. Similar investigation for y radiation of Cs¥’ 
in water has not been performed. In the above paper 
[13], the region near the source where F < 45 cm, 
which is of special interest in the design of irradiators, 
has not been investigated. Therefore, it was first 
necessary to determine the dependence of the accumu- 
lation factor on the distance F in water. This has been 
performed for a Co™ specimen with an activity of 952 
mg-equiv Ra. A comparison of the obtained results 
with data from [13] and calculation data is given in 
Fig. 1. For Co™ radiation in the F < 45 cm interval, 
in which we are interested, a good agreement between 
theoretical data from [9] and experimental results was 
observed. For the Cs!" radiation, it appeared that the 
theoretical values of the accumulation factor By, which 
were Calculated for coefficients from [9], are somewhat 
larger than the experimental values. These data must 
be obviously checked again for a larger interval of F 
values, 

The comparison of experimental and theoretical 
data for linear and a cylindrical irradiator with Co® 
radiation is shown in Figs. 2 and 3. Figure 2 provides 

8, 
12 
if 


10 
9 


8 
7 
6 
5 
4 
3 
2 
1 
0 


10 20 38 40 50 60 70 80 90 100 110 Fykem 


Fig. 1. Accumulation factor By for water and y ra- 
— of Co™ and Cs!" point-sources. -X-) cs!8"; 
@) Co®; O) data from [13]; —) calculation with res- 
pect to data from [9]. 





data for linear sources which are 40 and 80 cm long. 
It is obvious from the figure that an entirely satisfac- 
tory agreement between theoretical and experimental 
data can be observed when (3) is used for calculations. 
Figure 3 shows experimental data obtained in determin- 
ing the dose rate at the center of a cylindrical irra- 
diator consisting of eight rods, each having a length 
of 40 cm, which were arranged along the cylinder 
generatrices. The rods could be shifted simultaneously, 
whereby the irradiator diameter would be changed [3]. 
The irradiator was placed in water (in a well 3.5 m 
deep and 1.5 m in diameter). A chamber, made of 
2-mm thick iron with a diameter and height of 30 cm, 
which was filled with water, was also placed in the 
irradiator. The measurements were performed at the 
irradiator center without the chamber (directly at the 
center in water) and at the chamber center. Consi- 
dering that the irradiator diameter was changed from 
40 to 100 cm and that the radiation scattered by the 
well walls thus passed through a distance in water of 
not less than 1 m before reaching the measurement 
point, the medium extent could be considered as 
infinite. 


The measurement results obtained for the irradiator 
center were divided by 8 for a comparison with the 
results obtained in measuing the dose rate from a 
single rod. 

The good agreement between theoretical and 
experimental data for the central point of the irradiator 
in a homogenous medium is still not a sufficient 


indication that (3) and (4) are valid for all cases which 
can be encountered in practice. Actually, in the case 
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Fig. 2. Dependence of the radiation dose 
rate of Co™ linear sources in water on 
distance (F, cm) for q = 2.2 g-equiv Ra/ 
cm: 2L = 40 cm(— by calculation, x by 
experiment); 2L =80 cm (--~by calcu- 
lation, O by experiment). 


where the dosage field is calculated for air, the 
absorbing layer of water inside the chamber occupies 
only a portion of the distance from the source to the 
measurement point, and the radiation from different 
sections of linear sources will pass through different 
thicknesses of water and air layers. In the case where 
the cylindrical irradiator is in air and the chamber is 
filled with the object to be irradiated (in our case, 
water), the radiation passes through an absorbing layer 
which cannot be considered as infinite. Thus, the 
amount of multiple-scattered radiation is smaller than 
in the first case. 

Thus, we have four nonequivalent cases with 
respect to radiation scattering and absorption conditions 
which can be encountered in practice. Conventionally, 
we shall call them the water — water, water — air, air — 
water, and air — air cases. 

Expressions (3) and (4) for the water — air and air— 
water cases yield the same results in the case where 
R= 1/2F, where R is the radiation path in water. 

The dependence of the dose rate on the irradiator 
diameter was investigated for all four cases. Measure- 
ments by the ionization and the chemical methods 
were performed at two points inside the chamber (d = 
= 30 cm;/ = 30 cm): at the chamber center and at 
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Fig. 3. Comparison of dose 
rate measurements for radia- 
tion from a linear source and 
for radiation at the center of 
a cylindrical Co™ irradiator 
in water with calculation data. 
X) Linear source: 2L = 40 cm; 
©) measurements without a 
chamber (1/8Pp at the irradia - 
tor center); O)measurements 
in the chamber (1/8 Py at the 
irradiator center); —) theore- 
tical data. 





the distance of 12 cm from the center. This made it 
possible to determine the manner in which the dosage 
field changes with an increase in the irradiator diameter 
for various cases of irradiation. For the air — air case, 
the measurements were performed only for the central 
point (C = 0), which was due to the fact that in 
measurements inside a well 1.5 m in diameter a con- 
siderable amount of radiation scattered by the well 
walls cannot be avoided. 

In order to establish whether theoretical and experi - 
mental data agree under conditions where the pre- 
sence of a considerable amount of scattered radiation 
was eliminated, the dose rate from a single linear 
source (2L = 40 cm, q = 2.05 g-equiv Ra/cm) was 
measured in air outside the well, The measurement 
results are shown in Fig. 4. The measurements in air 
at the irradiator center (Fig. 5) showed that the amount 
of radiation scattered by the well walls was approxi- 
mately equal to 7-10%, Therefore, in further experi- 
ments for the air — air case, we used calculation data 
for the dose rate at the points C =0 and C=12 cm. 

The measurement results for the water — air, 
water — water, and air — water cases are shown in Figs. 
6, 7, and 8, respectively. It is obvious from these 
figures that, for practical purposes, an entirely satisfac - 
tory agreement between calculation and experimental 
data was obtained. However, for the air — water case, 
a deviation of 10-12% was observed (due to a smaller 
amount of scattered radiation than was the case for 
the water — air combination). For practical purposes, 
the indicated discrepancy is not very important. 

In the majority of cases, (3) and (4) for a linear 
source provide the possibility of designing all irradiator 
shapes which are important in practice (cylinder, 
plane, two planes, and a system of linear sources), 
since all these irradiators are usually composed from a 
larger or smaller number of rods filled with radioactive 
preparations. Therefore, calculations with respect to 
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Fig. 4. Dependence of the dose rate in 
radiation from a Co™ linear source in 
air on distance for 2L = 40 cm and q = 
2.05 g-equiv Ra/cm. —) by calcula- 
tion; x) by experiment. 


individual rods, especially in those cases where they 
are not numerous, provide more accurate results than 
calculations where equations for a continuous plane and 
a hollow cylinder are used. 


Determination of the Most Suitable Method for the 





Irradiation of Objects 





The aim in designing powerful y devices for radia- 
tion treatment is,in the first place, the determination 
of the amount of radioactive substances necessary for 
creating a dosage field of an assigned magnitude and 
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Fig. 5. Dependence of the radiation 
dose rate inside a Co™ cylindrical 
irradiator in air on the irradiator dia- 
meter. —) By calculation; +) by ex- 
periment (C= 0), 
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Fig. 6. Dependence of the radiation dose 
rate in air inside the chamber for objects 
(d = 30 cm) on the diamter of a Co™ 
irradiator which is immersed in water. 
—) By calculation, +) by experiment 
(for C= 0); —) by calculation,O) by 
experiment (for C=12). 





uniformity in a given volume. Therefore, it is 
extremely important to determine which irradiation 
method secures the greatest efficiency of isotopes. For 
this, it is necessary to compare the dependence of the 
degree of the dosage field nonuniformity on irradiator 
parameters for different cases with the dose rate 
obtained at the irradiator center. In this manner, it is 
possible to determine the cases where the amount of 
isotopes can be reduced if the assigned irradiation non- 
uniformity and the dose rate are attained. 

Figure 9 shows the dependence of the degree of 
the dosage field nonuniformity on the diameter of a 
cylindrical irradiator for different cases. Therefore, 
experimental data (see Fig. 8, the dash-dot line) are 
used in calculating the nonuniformity percentage for 
the air — water case, and calculation data are used 
for the air — air case. 

As was mentioned above, the degree of nonunifor- 
mity was defined as the ratio of the dose rate at the 
point C = 12 cm to the dose rate at the point C= 0 
For a certain given degree of the dose field nonunifor- 
mity, different irradiator diameters will be required 
for different irradiation cases, whereby different abso- 
lute magnitudes of the dose rate at the irradiator 
center will be secured. Naturally, irradiation methods 
where, for a given dose field uniformity, the maximum 
dose rate is secured are the most suitable. If we accept 
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Fig. 7. Dependence of the radia- 
tion dose rate in water inside the 
chamber for objects (d = 30 cm) 
onthe diameterof a Co™ irradia- 
tor which is immersed in water. 
—) By calculation, +) by experi- 
ment (for C=0); —) by calcula- 
tion, 0) by experiment (for C=12). 


a degree of nonuniformity which is often encountered 
in practice ( Py/P») 100 = 110%), then, for different 
cases, we obtain: 
Air- Air. . . D = 71.5 cm and Pp = 66 r/min (by cal- 
culation); 
Water — air. . . D = 48.5 cm and Py = 142 r/ min (by 
calculation and experiment); 
Air — water. . . D> 100 cm and Py< 31 r/min (by cal- 
culation; Pp < 28 by experiment); 
Water- Water. . . D = 71.5 cm and Py = 34 r/ min (by 
calculation and experiment). 
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Fig. 8. Dependence of the radiation 
dose rate in water inside the chamber 
for objects ( d=30 cm) on the diameter 
of a Co™ irradiator which is in air. 

—) By calculation, -x .-x.-) by experi- 
ment (for C=0); —) by calculation, 

O) by experiment (for C=12). 
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Fig. 9. Dependence of the dosage field non- 
uniformity in a cylindrical irradiator on its 
diameter. 





Thus, the dose field uniformity increases when 
an additional water layer outside as well as inside the 
object is present. For the same assigned dose field 
nonuniformity, this makes it possible to reduce the 
irradiator diameter, and, consequently, to secure either 
a larger dose rate or to save the radioactive substance. 
For a given device, in the cases where the objects 
slightly distort the dose field in air (as can be seen from 
a comparison of the water — air and air — air cases), a 
double saving of radioactive substances can be secured. 
The authors of [4] obtained similar results in designing 
an irradiator consisting of two flat surfaces; they noted 
the usefulness of applying either additional flat filters 
or of increasing the thickness of flat isotope sources. 

Thus, the radiation treatment of objects requiring 
a high degree of irradiation uniformity should be 
performed in the presence of specially designed addi- 
tional filters. This irradiation method would be espe- 
cially suitable for objects of small sizes and densities 
i.e., for obj ects which cause slight distortions of the 
dose field in air. 


Applicability Ranges of Co™ and Cs" Isotopes for 
Irradiation 





An important problem in irradiation practice is 
the right choice of the isotope to be used as a y-ra- 









































-3 
. 0 10 20 30 





R,cm 


Fig. 10. Attenuation of broad (full lines) 
and of narrow (dotted lines) beams of y - 
rays from point and linear sources in water. 
1) Co™ point-source, broad beam, Q = 1 
g-equiv Ra, e) by experiment; 2) Cs!®” point 
source, broad beam, Q = 1 g-equiv Ra; 

x by experiment; 3) Co™ linear source; 
broad beam; Q = 100 g-equiv Ra, 2L = 100 
cm; 4) Cs"8" linear source: broad beam, 

Q = 100 g-equiv Ra, 2L = 100 cm; — - -) 
narrow beams. 


diation source. As was indicated above, Co™ and Cs!*" 
are the most promising isotopes atthe present time, 
however, no quantitative comparisons between these 
isotopes have been made. 

Figure 10 provides a comparison of the dose rate 
dependences on distance for point and linear Co™ and 
Cs"*" sources for broad and narrow beams. The dose 
rate was calculated according to (1) - (4) for a homo- 
genous medium (water), where the activity of point— 
sources was assumed to be equal to Q= 1 g-equiv Ra, 
and the linear specific activity of linear sources was 
assumed to be q = 1 g-equiv Ra/cm. 

It is obvious from the figure that for distances 
of less than 40 cm, the dose rate values for Cs!°” which 
were obtained by calculation are somewhat larger 
than the dose rates obtained with Co™ sources. The 
experimental dose rate values in this interval are equal. 
From a comparison of the dose rate dependences on 
distance for Co™ and Cs!*" sources, we can draw the 
conclusion that the application of csi preparations 
(with an equivalent radiation yield) for objects with 
a thickness of up to ~ 40 cm is more suitable than 
the application of Co™ preparations, which is important 
in practice. 

In practice, it is most often necessar y to irradiate 
objects of a thickness less than 40 cm and of a density 
equal tol or less. Therefore, for a great majority 
of y devices, it is more convenient to use Cs!*” pre - 
parations because of the longer half-life of Cs"*7 and 
a smaller required shield thickness, especially if heavy 
material is used. 

The use of Co™ is rational only in the case where 
it is necessary to irradiate large-size objects of great 
density or if very high dose rates, which depend on the 
magnitude of the specific activity of preparations, are 
to be secured. 

SUMMARY 


1. In the design of irradiators of y devices, for 
all cases encountered in practice, the consideration 
of multiple scattering by means of the accumulation 
factor, expressed by the sum of two exponentials, 
secures a sufficient accuracy. 

2. The most suitable method of irradiating objects 
which require the presence of a uniform dose field is 
the use of additional filters which equalize the dose 
fields ofextended sources. These filters must be cal- 
culated for each individual case. 

8. For the irradiation of objects with a thickness 
of less than 40 cm and a density not exceeding 1, the 
use of Cs'*7 preparations as y - radiation sources is 
much more efficient than the use of Co™ preparations. 
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GENERAL RELATION FOR CRITICAL THERMAL LOADING 
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Translated from Atomnaya Energiya, Vol. 7, No. 3, pp. 252-253 


September, 1959 
Original article submitted July 10, 1959 


Analysis of similar conditions involved in heat- 
transfer processes in various liquids shows that for com- 
plete equivalence of the heat-transfer processes it is 
necessary that these liquids: a) be thermodynamically 
alike (i.e., satisfy the identically applied conditions); b) 
be under corresponding conditions; and c) have equal 
values for the equivalent criteria which characterize 
the conditions existing at the liquid-solid boundary. 

The requirement of thermodynamic equivalence 
is especially important in those cases where on the 
one hand, the condition parameters vary within wide 
limits, and on the other hand, a molecular mechanism 
for the effect appears more likely, undistorted by 
complex external conditions acting on the liquid. This 
applies to the boiling of a liquid, and particularly to 
the crisis arising during the changes involved in the 
boiling process. 

It is not hard to show that the basic characteristic 
of the boiling crisis — the critical heat flow level (or 
the critical thermal loading) — must, in the case of a 
liquid boiling in a large volume, be expressible as a 
function of the following type: 


_ gapt/p, To’ ( P ve] 

Ae pi/2 ce Tix’ Bie 

where R is the gas constant ; p is the molecular weight 
of the substance; p,, Vg, and T, are the critical 
parameters of the substance; a is the average radius 

of curvature of the heating surface irregularities; p/P, 
is the applied pressure; f is a certain function which 

is identical for a group of thermodynamically equiva- 
lent substances. 

The indicated relationship for q, is obtained as 
follows, For thermodynamically equivalent substances, 
the ratio of the quantity of heat received in a given 
process by one mole of liquid to the absolute temper- 
ature of the liquid turns out to be a universal function 
of the applied parameters and of the physical charac- 
teristics of the process. The average size of vapor 
bubbles formed during boiling of the liquid appears to 
be such a characteristic; it is determined by the mag- 
nitude of the heating surface irregularities, so it follows 
that a factor ¥," 8 /a is introduced into the process 
equation. 

The expression for qc can be shown in the form 
of a product of a dimensional factor times a dimen- 
sionless function of the reduced parameters (in this 
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case itwill be: either p/p,or T/T, ) and the term 

v,i/ 3/a. The dimensional factor must consist of basic 
physical parameters of the given substance; namely, 

the critical parameters pg, V., and T,,, which are the 
most important thermodynamic characteristics of the 
substance, expressing themselves in the general quali- 
tative form as the effect of molecular forces and the 
molecular masses (i.e., the relation u/g). Further- 
more, the gas constant R may enter into the dimensional 
factor. Combining the quantities Per Tc, H/g, and R 

in such a way as to obtain the dimensions of the thermal 
flux, we obtain the above expression for qc. 

It follows from this expression that for thermody - 
namically equivalent substances, boiling under iden- 
tical conditions, and having equal (or close) values of. 
Vo, the quantity (q, w1/%) / (Po T,!/%) must be the 
same for equal values of p/p,. 

This result is in good agreement with experiment. 
The figure shows the dependence of the quantity (q, 
Tihs /®, T2A) on p/p, for benzene and pentane. * 
Benzene and pentane have equal values of the critical 
coefficient (R TIMP V,) = 3.76, i.e., they are ther- 
modynamically equivalent; their critical volumes V 
(256 and 310 cc/ mole) are very close. As can be 
seen in the figure, the points for benzene and pentane 
lie on the same curve. 
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Curve showing equality of the values of q uth / 

/ PoT,' 2 with identical p/p.values for the thermo- 
dynamically equivalent substances; O) benzene; 0) 
pentane. 


* The values for q, Were taken from the work of M. 
Cichelli and C. Bonilla, Trans, Am. Inst. Chem. Eng. 
41, 745 (1946). 
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HEAT-TRANSFER COEFFICIENT 
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September, 1959 
Original article submitted December 12, 1958 


An internal heat source appears in a liquid flow 
carrying a radioactive suspension, in the flow of a radio- 
active solution, when electric current passes through a 
stream of liquid metal. and in similar cases. Let us 
consider the influence of this source on the heat-transfer 
coefficient when a liquid flows in a cylindrical tube with 


a length sufficient for neglecting the inlet section effect. 


We shall use the following notation: 

t is the temperature; 

R is the present radius; 

Ro is the tube radius; 

x is the coordinate in the tube axis direction; 

d is the thermal conductivity coefficient; 

a is the temperature conductivity Coefficient; 

Ww is the flow velocity; 

qy is the internal heat source density; 

V = (try tY(ttw- t) is the dimensionless temperature; 

Po=2W Ry /ais the Peclet number; 

w is the mean-discharge flow velocity; 

Nu = 2QRy /A is the Nusselt number; 

t is the average flow temperature with respect 
to enthalpy in a given tube cross section; 

try is the tube wall temperature; 

w = w/w is the dimensionless velocity; 

Z = qyRo/2dtw is the relative heat source density; 

— = R/R, is the relative radius; 

a is the heat-transfer coefficient; 

9 is the heat-transfer coefficient for qy = 0; 

Atw is the heat flow through the tube walls; 

X = x/Rz is the relative axial distance. 

After writing the heat propagation equation [1] 


aL C+9 )® ox ]+ 


we integrate this equation within the limits from 0 to 
R, assuming that q, = const. After transformation into 
the dimensionless form, we obtain 


+i wR ot 
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dy Ro 5\ Hb ade 
4% 
pa 


As a first approximation, let us introduce in this ex- 
pression,which is the average value with respect to the 
cross section,the value dv /0X, 


By determining the mean-discharge temperature, 
we obtain 


v=1=2\ wEvdE. (3) 
0 


By determining the heat-transfer coefficient from the 
usual equation 


(4) 


we find that 


Thus, the degree of influence of an internal heat 
source is proportional to the difference 


E 
2 
A=> wkdg—&. (6) 


For w = 1, the difference (6) is equal to zero, i.e., in 
a flow with maximum turbulence the heat-transfer 
coefficient does not depend on the internal source den- 


sity. For laminar flow, i.e., for a parabolic velocity 
distribution, 


Be Bo (7) 
G  1-+-0,272 Z° 





For turbulent flow with the velocity distributed accord- 
ing to the 1/7 law and for Pr = 0, 


EL (8) 
a, 1-+0,0834 Z’ 


For q, = 0, (5) becomes the Lyon integral ratio [2]. 
The figure shows the results of calculations where 
the above equations are used. It is obvious that the 
influence of an internal source on the heat-transfer 
coefficient is slight and that it becomes manifest for 
w#1. Actually, if heat is supplied through the walls 
to a flow where heat is generated (Z> 0), the layers 
with a lower velocity are heated more intensively. In 
this case, the temperature difference t,,, -t is greater 
than in thecase where q, = 0. Consequently, the heat- 
transfer coefficient decreases with an increase in the 
source density. When a heat sink (Z < 0) is present 
in the flow, the heat-transfer coefficient increases 
with an increase in the internal source density. 
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Influence of the internal heat source density on steady 
heat transfer in tubes, 1) Flow with maximum turbu- 
lence(w=1); 2) laminar flow; 3) turbulent flow where 
the velocities are distributed according to the 1/7 law. 


LITERATURE CITED 


1. S&S. S. Kutateladze, Principles of the Heat-Transfer 
Theory [in Russian] (Mashgiz, Moscow, 1957). 
2. R. Lyon, Chem. Eng. Progr. 47, 11 (1951). 





HEAT TRANSFER IN TRANSVERSE FLOW OF MOLTEN 


SODIUM AROUND A SINGLE CYLINDER 


A. A. Andreevskii 


Translated from Atomnaya Energiya, Vol. 7, No. 3, pp, 254-256 


September, 1959 
Original article submitted November 25, 1958. 


The heat transfer in tubes in the transverse flow of 
liquid metals has been meagerly investigated. Some 
information on this problem can be found in [1-3]. 
After the present work was completed, [4 and 5] were 
published in which some experimental data on heat 
transfer in the flow of mercury around a group of tubes 
arranged in checkerboard order were given. 

The device for investigating the heat transfer to 
sodium, which has been used in this work, consists of 
a closed heated circuit with forced liquid metal cir- 
culation. Due to the fact that the circuit was air-tight, 
it was possible to evacuate it and to fill it with an 
inert gas under a pressure of 0.1-0.2 atmos. Also, a 
bypass with a grid filter (10,000 meshes /cm?) was 
provided in the circuit for the trapping of oxides. * 

The working section of the device consisted of a 
rectangular box with an inlet step diffuser and an outlet 
convergent nozzle. The hydraulic characteristics of 
the working section inlet were investigated in special 
experiments which were performed with water. An 
electrocalorimeter, the frame of which consisted of a 
hollow Steel 20 cylinder with an outside diameter of 
22 mm and a working length of 200 mm, was installed 
in the middle portion of the working section. A steel 
tube, which was insulated from the calorimeter by 
mica, served as the heater. In order to secure a uniform 
heat flow along the calorimeter circumference, special 
attention was paid to securing a uniform thickness of 
the mica insulation. The thermal load of the calori- 
meter surface under operating conditions was 7-9 x 104 
kcal/ m?- hr, which guaranteed a reliable measurement 
of the basic quantity to be found — the temperature 
drop from the wall to the liquid. Six grooves (every 60°) 
were cut on the internal surface of the calorimeter 
frame for the placement of chromel-alumel therco- 
couples. The thermocouple hot-junctions were clamped 
under metallic blocks made of the same steel as the 
calorimeter frame. The thermocouple tips were located 
in the cross-sectional plane passing through the middle 
of the calorimeter, which permitted the simultaneous 
measurement of the heat transfer at six points along 
the calorimeter circumference. 

An evaluation of errors showed that the maximum 
total measurement error arising in determining the heat- 
transfer coefficient did not exceed 410%, 


The presence of oxygen(<1.0x 10 %), iron 
(< 5.6 x 104%), and nickel (< 4.1x 10%) in molten 
sodium was detected by chemical analysis. 

The experiments were performed for sodium tem- 
peratures from 220 to 250° C, and the velocity interval 
from 0.1 to 0.4 m/sec, which corresponded to an in- 
terval of Reynolds numbers from 4000 to 20,000 and an 
interval of Peclet numbers from 25 to 125.1 The ob- 
tained Nu=f (Pe) relation is shown in Fig. 1. This re- 
lation was obtained for ¢ = 120° on the basis of read- 
ings of thermocouples Nos. 4 and 6, which were mutually 
symmetrical with respect to the flow axis. 

A comparison of experimental data obtained for 
different angles, shows that the relations between the 
local heat transfer and the velocities at the front and 
the back cylinder regions are different. 

In the back cylinder region (from ¢ = 120 to 180°), 
the relation between heat transfer and ve locity is des- 
cribed by a logarithmic straight line with a slope close 
to 0.35. On the front cylinder surface (y= 0 and 60°), 
the dependence of heat transfer on velocity has a more 
complex character. For Peclet numbers > 50, the 
experimental points are grouped around a curve with 
an exponent close to 0.5. In the interval of Peclet 
numbers from 40 to 50 (Pe ~ 8000), the heat-transfer 
curve has a bend, and for Pe> 50, the heat-transfer 
coefficient depends only slightly on velocity. The 
causes of such a change in the character of heat transfer 
are not known at the present time. However, it can be 
assumed that the heat-transfer coefficient at the cylin- 
der front zone is more sensitive to changes in flow 
conditions than the heat-transfer coefficient at the back 
zone. 


* For more details, see [1 and 6]. 


ft In calculating the Pe and Re numbers, the upstream 
flow velocity was taken for the detérmining velocity, 


and the cylinder diameter was taken for the determining 
dimension. 


t For 9 = 60° as well as for y = 0°, the line averaging 
the experimental points has a bend in the region of 
Peclet numbers from 40 to 50. Fhe dependence of the 
heat transfer on velocity for g = 180° is similar to the 
lower line in Fig. 1. 
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The heat-transfer distribution around the cylinder 
perimeter is shown in Fig. 2. The points which were 
used in plotting the heat-transfer curves were obtained 
for the corresponding angle ¢ and Pe number (with 
respect to the averaging lines of initial Nu = f (Pe) 
curves for y = const). It is obvious from the figure that 
the heat-transfer coefficient reaches its maximum value 
at the leading cylinder generatrix (y= 0°). The heat- 
transfer minimum is observed in the trailing zone (¢ = 
=180°). 

In Fig. 3, the heat-transfer curves for sodium (Pr 
= 0.007) and air(Pr ~0.7) flow around a single cylin- 
der are compared. The heat-transfer distribution curves Fig. 2. Heat-transfer distribution along the peri- 


for transverse flow of liquid metal around a bunch of meter of a single cylinder. Measurement points 
tubes arranged in checkerboard order have a similar with thermocouples: O) No. 1; @) No. 3; A) No. 4; 
character [1, 3, and 4]. A) No. 5; @) No. 6. 
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Fig. 3. Variation of the relative heat-transfer co- 
efficient along the perimeter of a single cylinder. 
Fig. 1. Dependence of the local heat transfer on 1) Sodium (O-Pe = 0,25, Re = 4,000; A-Pe = 125, 
sodium velocity. A) Thermocouple No. 5 (¢ = 0°); Re = 20,000); 2) air (Re = 16,000), 
@) thermocouple No. 6 ( ¢ = 120°); O) thermocouple 
No. 4(¢ =120°). 
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Fig. 4. Average heat transfer in transverse flow of liquid metal around a single cylin- 
der and around a bunch of tubes arranged in checkerboard order. 1) Single cylinder in 
sodium flow (Pr = 0.007); 2) bunch of tubes (data from [5]); A)Pr = 0.0215; O) Pr= 
=0.0185; A) Pr = 0.0160; x) Pr= 0.0138. 





’ The absence of a heat-transfer maximum in the 
trailing zone of a single cylinder and of tubes arranged 
in a bunch in transverse flow of liquid metal can be 
explained by a weaker hydrodynamic influence of the 
flow on heat transfer due to the high thermal conducti- 
vity of media with small Prandtl numbers. 

The dependence of the average heat transfer of 
a single cylinder on liquid, metal velocity is given in 
Fig. 4. in Nu vs. Pe coordinates. The obtained experi- 
mental data for the 50 = Pe = 125 interval is well 
described by the dependence. 

For the sake of comparison, Fig. 4. shows also the 
experimental data for mercury, which was . obtained by 
Rickard, Dwyer and Dropkin. 

It can be stated that a qualitative agreement exists 
between the results obtained in the present paper and 
the data from [5]. Actually, it is obvious from Fig. 4 
that, in flow of liquid metal around a bunch of tubes 
as well as in the case where media with Pr ~ 1 flow 
around tubes, a higher (by 20-40%) heat-transfer inte- 
sity is observed. The dependence of the heat-transfer 
coefficient on liquid metal flow velocity is stronger 


for a bunch of tubes (the exponent is ~ 0.6-0.65). In. 
view of the fact that the results of the present investi- 
gation are confirmed by data from [5], the author 
considers that for the 50 = Pe = 5000 interval, the heat 
transfer in flow of sodium and mercury around a single 
cylinder can be calculated according to the proposed 
equation until new data becomes available. 
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The spectrum of gammas emitted from the active 
zone of the IRT reactor, described in [1], via one of 
the horizontal channels of the reactor has been mea- 
sured by means of a Compton electron magnetic spec- 
trometer. The gamma spectrometer used in an earlier 
experiment [2] was employed in the measurements. 

Results of the measurements may be seen in Fig. 1. 
The instrumental spectrum is presented there. It 
was obtained by subtracting from the measured coin- 
cidence spectrum (pair electrons for the given y spec- 
trum yield a very small contribution to the coincidence 
spectrum) the contribution from electrons of pairs 
dislodged with Compton-recoil electrons by y rays 
from the radiator section of the spectrometer. Figure 
2 gives the corrected gamma-ray spectrum. The 
value of v(Ey)Ey, where v (Ey) is the number of 
y quanta per unit energy interval, is plotted on the 
ordinate. The area of the spectral interval[v(Ey )Ey]- 
‘AEy is proportional to the energy carried off by y rays 
having energies ranging from Ey to Ey + AEy. In 
making the transition from the instrumental spectrum 
to the corrected spectrum, the spectral sensitivity of 


the instrument was taken into account, as well as y -ray 
absorption in the passage of gamma quanta from the 
reactor core to the radiator of the spectrometer. In 
particular, traversal by y rays of 8 meters of air and 

9 cm of paraffin mixed with boron carbide was taken 
into account, 

The general outline of the experiment is repre- 
sented in the upper right-hand corner of Fig. 2. It is 
clear from the diagram that the y rays emerging from 
the reactor core traverse ~ 10 cm of graphite, which 
surrounds the core. No correction was applied for 
possible absorption of y rays by graphite in that layer. 
Absorption of y rays in that layer brings about a sharp 
reduction in the y spectrum at energies < 0.7 Mev. 

The IRT reactor burns enriched uranium; the 
moderator is ordinary water, and the basic structural 
material is aluminum. y lines due to neutron capture 
in aluminum (Ey = 7.72 Mev, etc.), carbon (Ey = 4.95 
Mev), hydrogen (Ey = 2.23 Mev) and y lines of energies 
1.78 Mev, associated with 8 decay of Al”, may there- 
fore be delineated with certainty in the y spectrum 
of the reactor. At energies in the neighborhood of 
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Fig. 1. Instrumental y-ray spectrum for y rays emitted from the core of the IRT reactor. 
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Fig. 2. Adjusted spectrum of y rays emitted by core of IRT reactor. 


1.4 Mev, a diffuse peak of uncertain origin appears. It 
might possibly be caused by y rays emitted by radio- 
active elements which are fission products of uranium. 
The major portion of the unresolved part of the y 
spectrum must be determined by y rays originating in 
the reaction (n, y) involving u?®* and u?* (the binding 
energy of the last neutron in u*™ and U** is 4.63 Mev 
and 6.41 Mev, respectively), and by y rays produced 
by fission. 

Note the substantial difference in the y-ray spectra 
of the RFT reactor [3] and the IRT reactor. Since the 


basic structural material used in the RFT reactor is 
stainless steel, the contrast between the y-ray spectrum 
of this reactor and that of the IRT reactor is due to the 
strongly developed high-energy portion of the spectrum, 
due to neutron capture events in the iron, nickel and 
chromium present. 

In the table below, we present relative values of 
Jv(Ey) Ey dEy found from Fig. 2, for various intervals 
in the spectrum. These values are cited separately for 
y lines of energies 2.23, 1.78, and 7.72 Mev, while 
the value for the first y line is taken equal to 100. 
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In conclusion we should like to take this opportu- 
nity to express our gratitude to the staff which made 
possible the operation of the IRT reactor during the 
measurements. 
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The variation in the radiative capture cross sec- 
tions for neutrons reacting with gallium isotopes was 
measured by comparison with the variation observed 
in the cross section of the nuclide I*"", used as a refe- 
rence standard. Samples of the test isotope and of 
the standard were irradiated simultaneously. while 
exposed to the same fast neutron flux. The 8 activity 
induced in the samples was measured by end-window 
counters. After a drop in activity, the samples were 
irradiated in a uniform flux of thermal neutrons, and 
the induced activity was measured using the same 
counters. 

If the time of irradiation, fluorescent lifetime, 
and measurement interval are kept constant for each 
sample during irradiation by fast and thermal neutrons, 
the capture cross section may be determined from the 
relation 
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thermal and fast neutron fluxes; of, of are the cap- 
ture cross sections of fast and thermal neutrons of the 
standard; ofh _is the capture cross section of thermal 
neutrons of the nuclide under investigation. 

A fast neutron source was provided by the reaction 
T (p, n)He®, which was brought about using a Van de 
Graaf accelerator. Over the 200-1,400 kev energy range, 
the energy spread of the neutrons amounted to.+ 30 
kev, while it amounted to + 50 kev over the 1,400- 
2,000 kevrange. Irradiation by thermal neutrons took 
place in the thermal channel of the fast neutron experi- 
mental reactor. In order to eliminate resonance and 
fast neutron effects, recourse was had to the cadmium- 
difference measurement technique. 

In the calculations, the capture cross sections for 
thermal neutrons by nuclides 17, Ga® and Ga” were 
taken to be 5.6 + 0.3, 1.4 + 0.3 and 4.0 + 0.7 barns, 
respectively [1]. The capture cross sections for fast 
neutrons by I" were taken from [2]. The mean stan- 
dard deviation of the experimentally measured values 
of K did not exceed 3-4%, The results of the measure- 
ments are plotted on the accompanying diagram. The 
radiative capture cross section for Ga® reveals a smooth 
energy dependence, while that of Ga’! manifests a 
sharp drop in the region ~ 550 kev. It is evident that 
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Variation in the radiative capture cross sections of fast neutrons reacting with gallium 
nuclides, as a function of neutron energy. 





this effect is due to the presence of inelastic scattering 

of neutrons from Ga” at energy levels ranging from 510 

to 610 kev [3]. The competition introduced by the ai 

inelastic scattering leads to a fall-off in the radiative 

capture cross section graph. 2. 
The authors would like to express their gratitude 

to A. I. Baryshnikov, V. 1. Zotova, and V.G. Nedopekin, 3. 

for their kind assistance in the execution of this research 

effort. 
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This study of y emission arising in the capture 
of the thermal neutrons by chlorides was performed on 
a beam from the reactor of the First Atomic Power 
Station [1]. The chloride specimen consisted of che- 
mically pure PbCl, powder. Scintillation counters 
were placed at right angles to each other and to the 
thermal neutron beam. A layer of boron carbide, 
3 mg/cm? thick, was placed between the crystals 
and the target, to absorb any thermal neutrons scattered 
by the target. NAKT1) crystals(@ 30, h = 20 mm, and 
@ 40, h = 50 mm) were combined with an RCA - 5819 
type photomultiplier tube. The pulses of the channel 
to be analyzed were directed to a 100-channel pulse 
analyzer, described in[2], and the pulses of the con- 
trol channel were fed to a differential discriminator. 
The pulses of both channels were fed in parallel to a 
fast coincidence circuit having a resolving time T = 
0.8+ 10°7sec. The value of the resolving time was 
secured by the absence of true coincidence losses, 

The required amplitude interval is delineated 
by the discriminator circuit in which the slow coin- 
cidence circuit (T = 1078 sec) is to be found, and in the 
presence of a pulse originating in the fast coincidence 
circuit, the pulse of the control channel gates the 
analyzer input. The pulse-height spectrum recorded 
by the pulse analyzer must therefore contain peaks 
corresponding to the y photons emitted simultaneously 
with the y photons recorded in the control channel, 

In order to break the complex spectra down into 
their components, graphs were plotted of the shape 
of lines from y sources with monochromatic energies 
inthe range 0.128-2.76 Mev, for a crystal 30 mm in 
diameter and 20 mm in height. The absolute effi- 
ciency (ratio of area under the total absorption peak to 
absolute number of y quanta from the source) as a 
function of the energy of the y quanta in the energy 
interval 0, 08-2,76 Mev was determined for the same 
crystal. The absolute activity of the sources was 
determined by the y-y coincidence method in an 
energetic analysis of coincident y photons. 

The spectrum of random coincidences was computed 
on the basis of the formula 


Nran (£E) = 2T N,N2(E), 


where T is the resolving time of the fast coincidence 
circuit, Ny is the speed of light in the control channel, 
being a constant magnitude for the measurement in 
question, and N»(E) is the speed of light in the channel 
being analyzed, this speed being taken from the coin- 
cidence-less spectrum, Preliminary measurements 
showed that the random-coincidence spectrum calcu- 
lated on the basis of the above formula is inclose ag- 
reement with the spectrum measured while introducing 
into one of the channels, a large delay, sufficient to 
eliminate true coincidences, 

The background of random coincidences remained 
below the 5% level in the measured coincidences in all 
regions of the spectrum (with the exception of the 
region of energies lying below 511 kev), 

In addition to random coincidences, spurious 
coincidences were also taken into account, these 
being due to the following causal phenomena; 1) 
Compton scattering of y photons from one crystal 
into another; 2) impingement on a crystal in the test 
channel of annihilation y radiation accompanying the 
formation of an electron-positron pair in the crystal 
in the control channel, 3) the contribution to the pulse 
spectrum of pulses singled out by the window of the 
differential discriminator of the Compton distribution 


‘ attributable to y rays of greater hardness. 


For each fixed amplitude interval selected in the 
control channel, the coincidence spectrum from the 
test channel was investigated by several different 
approaches, by measurement of the gain factor of the 
corresponding amplifier, The coincidence spectra so 
obtained were corrected for random coincidences and. 
were then plotted on a common graph, after reduction 
to constant conditions, Points indicated by similar 
markings on the accompanying graphs represent mea- 
surements carried out with gain factor fixed. The 
time of measurement (5-40 hr) was chosen such that 
the number of true coincidences in the peaks would 
not be less than 100 pulses per channel. Where gain 
instability in any channel was found to exceed 1%, 
that measurement was rejected. 

The diagram shows the spectra of y photons emitted 
simultaneously with hard gammas, Diagram a depicts 
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a spectrum of gamma photons emitted simultaneously 
with y photons in the energy interval 7.4-8.5 Mev. The 
control channel was then used to select out y photons 
having energies of 8.58 and 7.79 Mev. Furthermore, due 
to the statistical dispersion, pulses stemming from the 
photopeak of y rays having an energy of 7.42 Mev can 
fall into the control channel. A peak corresponding to 
the energy 0.79 Mev stands out sharply in the spectrum, 
and a weaker peak at 1.16 Mev is also visible. This 
implies cascades of y rays, 7.79-0.79 and 7.42-1.16 
Mev. As we see from diagram (a), the y line of energy 
735 kev, reported in [3], does not appear in the spectrum. 
The shape of the photopeak at energy 790 kev shows that 
if there were y rays present with energy 735 kev, the 
intensity of these latter photons does not comprise 10% 
of the intensity of the 790 kev line. 

At energy ~ 0.5 Mev, a peak appears which is 
an instrumental artifact resulting from a pair effect 
in the crystal of the control channel when the crystal 
in the test channel was impinged upon by one of the 
annihilation quanta, 

If y rays with energies of 0,79 and 1. 16 Mev are 
emitted in a cascade by y rays with energies of 7,79 
and 7,42 Mev, respectively, the ratio of their inten- 
sities in the coincidence spectrum must be equal to 
the ratio of the intensities of hard y photons in the 
control channel. 

Using as point of departure the ratio reported in 
[4] for the intensities of the y lines emitted by cts 
the ratio of the y-line intensities in the control channel 
(taking line shapes into account) in the energy interval 
6.6-8.5 Mev (see diagram b), 17, 7 /\q.14 = 2.5. From 
the spectrum shown in diagram b, we find that Ip, »9/ 

I}, 16 = 2.6 + 1. The numerical agreement confirms 
the assumption made in [4] to the effect that cascades 
7.79-0.79 Mev and 7,42-1.16 Mev take place and 
proceed with no branching off through the 0.79 and 
1.16 Mev levels, respectively. The spectrum over 

the energy interval 6.1-8.5 Mev (see diagram c) shows 
a photopeak appearing in the control channel, and 

the beginning of a Compton peak due to y lines of 
6.96 Mev energy, with a gentle shoulder corresponding 
to the 1.60 Mev y line showing up in the spectrum of 
the test channelin the region of 1.60 Mev, the 1.60 

y line presumably [3-5] being found in the cascade 
with the y quantum of energy 6.99 Mev. In this region 
of energies, the control channel also embraces a 
portion of the spectrum of pulses belonging to the y 

line of energy 6.64 Mev; ay line of energy 1.95 

Mev appears correspondingly in the test channel spectrum 
which serves to confirm the 6,64-1.95 Mev cascade, 

It was found from the coincidence spectrum of the 
channel that I, ¢/1,, 95 = 0.18 + 0.07, while the ratio 
of intensities of hard lines, [4], with energies 6.99 Mev 
and 6,64 Mev found to accompany them in the cascade 
(with the contribution to the control window taken into 
account), must be Ig g9/ Ig, g4 = 0.13 + 0.02. The 
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numerical agreement in intensities confirms the pre- 
sence of the 6, 99-1.60 and 6.64-1.95 Mev cascades, 

When the threshold of the control channel is set 
at the position corresponding to energy 5.7 Mev (see 
diagram d), pulses corresponding to an energy of 6,11 
Mev also enter the channel, As is evident, this brings 
about a sharp increase in the peak in the energy region 
of 0.5 Mev, which may be due solely to the appea- 
rance of true coincidences of y quanta with energies 
of 6.11 and 0.5 Mev. 

Assuming that the contribution of the instrumental 
spurious peak at energy 0.51 Mev did not change, we 
find that half of the area included in this photopeak 
is due to the appearance of y quanta having energies 
of 0.5 Mev. This fits well with the existence of the 
complex cascade 6, 11-0. 51-0, 790-1. 16 Mev [4]. 
Earlier, as the width of the control channel was incre- 
ased, the appearance of a new y line in the test 
channel took place because of a decrease in the total 
intensity of the y lines, which were already present 
in the spectrum. In the case in point, the transition 
from the energy region of 6,1-8.5 Mev (see diagram c) 
to the region of energies 5,-8.5 Mev (see diagram d) 
does not result in a decrease of the total intensity of 
the y lines having energies 0,79, 1.16, 1.60 and 1,95 
Mev. This being so, a decrease, due to the appearance 
of a new y line of energy 0,51 Mev, in the total in- 
tensity of the y lines alluded to above is offset by the 
appearance, as suggested in [4], of y lines coinciding 
with it, and having energies of 0.79 and 1.16 Mev. 

In the transition to the next lower level at 4.9 
Mev (see diagram e), pulses belonging to energies of 
5.72 and 5,50 Mev also arrive in the control channel. 
The contribution attributable to y rays of energy 5. 50 
Mev may be neglected, In that case, a y line of 
energy 2.9 + 0.1 Mev will appear in the test channel 
spectrum. Its appearance is associated with the appea- 
rance of the 5,72 Mev y line. These y lines are 
obviously found in the cascade, 

The measurements reported here, as well as the 
investigations of coincidence spectra with y lines at 
0.79, 1.16, and 1,95 Mev, permit us to draw the 
following conclusions; 

1, There exist cascades corresponding to 7.79- 
0.79 and 7,42-1,16 Mev. Their intensities are 
interrelated as 7.8; 14, and the intensities of the 
cascades 6, 99-1. 60 Mev and 6,64-1.95 Mev are in 
turn related as 1.9; 14,4, in agreement with the 
supposition advanced in [4]. 

2, Cascades corresponding to 5, 72-2. 88 and 
6. 11-0. 51-0. 79-1. 16 Mev take place. The last- 
mentioned cascade is confirmed by the fact of coin- 
cidences of y quanta corresponding to 0. 79-6,11 and 
1.16-6.11 Mev. This last cascade actually passes 
through the 1.95 Mev level, which is borne out by 
the coincidence of y quanta having energies 1. 95 and 
6.1 Mev. The coincidences of y quanta having ener- 





gies 1.16 and 6.64 Mev would indicate that, starting 
with the 1,95 Mev, a transition is possible not only 
to the ground state, but that even a transition through 
the 0.79 Mev level or through the 1.16 Mev level 

is possible. 

The cascading transitions obtained show complete 
agreement with the scheme of y transitions compiled 
in (4). 

In conclusion, the author would like to take this 
opportunity to express his gratitude to Doctor of 
Physical and Mathematical Sciences A. K. Krasin 
for his kind and unflagging interest, valuable com- 
ments, and cooperation in carrying out the research, 
and to V, P. Radchenko for assistance kindly rendered 
in performing the necessary measurements, 
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The fundamental method used in the direct study 
of the fine structure of a decay still remains the mea- 
surement of the spectra of a particles by means of mag- 
netic spectrometers. The most promising spectrometers 
are those incorporating spatial focusing in an inhomo- 
geneous axisymmetric magnetic field [1-3]. Two 
existing & spectrometers of that type are known; one 
which focuses @ particles at an angle of a v¥2(4], and 
one which focuses alphas at an angle of 180°, using the 
fringe effect [5]. The first variant was selected by 
the authors in elaborating the design of our own & spec- 
trometer. 

The principal drawback in the instruments referred 
to is the rather small aperture (product of solid angle 
by source area), With a view toward improving this 
feature while retaining the high resolution of the spec- 
trometer, the radius of the instrument's central orbit 
was relatively enlarged (py) = 155 cm). 

Figure 1 shows a photograph of the © spectrometer 
so described. Its principal details are the magnet with 
field coils and the vacuum chamber of nonmagnetic 


stainless steel. The magnet enclosure is devised in the 


. 3, pp. 262-264 


form of a double mushroom-shaped assembly consisting 
of three parts; the core, the cylindrical shell, and two 
**hats". The width of the pole pieces is ~ 70 cm, 

and the pole gap is 35 cm. The magnet parts are made 
to relatively high tolerances (0.2mm). Specially 
machined cover plates dimensioned to analytically com- 
puted specifications are fitted to the pole pieces. The 
total weight of the facility runs to ~ 90 tons. The 
chamber volume is ~ 1000 liters. The chamber is 
pumped out by a high-vacuum VA-54 assembly and 

a forevacuum pump of the VN-2 type. The extent of 


evacuation of the vacuum chamber reaches several 
units of 10°° mm Hg. 


Sources of o-radioactive isotopes are placed in a 
special device, making it possible to perform alternate 
exposures of four targets without jeopardizing the 
vacuum. The maximum size of usable alpha sources 
is 100 X 10 mm. 


*Submitted at the Ninth All-Union Congress on Nuclear 
Spectroscopy (Khar*kov, January, 1959). 


Fig. 1. General view of the alpha spectrometer. 





Alphas emitted are recorded by a proportional 
counter or by thick-layered phtographic plates. A set 
of plates 500 mm long and 90 mm thick may be 
exposed simultaneously. The plates are placed in a 
special receptacle which may be used to permit the 
successive exposure of four sets of plates. Replacing 
of one set of plates by another and closing the recep- 
tacle is accomplished without affecting the vacuum. 

The magnet winding is fed via a selenium rectifier 
hooked up to a three-phase bridge circuit. The voltage 
across the rectifier is supplied from a synchronous 
motor-generator operating on alternating current at a 
35 kva power rating, and passes through a DN-35 satura- 
tion choke coil. The current is 700-1300 amp over the 
operational range of the instrument, which corresponds 
to a magnetic field intensity of 2, 0-3.5 koe, 

The basic circuit arrangement found in [6] is 
employed to stabilize the magnetic field. A 5LO-38 
type CRT placed between two fixed "reference" mag- 
nets is utilized as a component,sensing variations in 
magnetic field intensity. This system is located in the 
scattered magnetic field of the spectrometer. 

In the case of spectrographic measurements, the 
peak deviation of the magnetic field from a given value 
Hp did not exceed 2+ 10% in 8 hr of continuous opera- 
tion of the instrument. Measurement of the value of the 
magnetic field intensity is based on the phenomenon of 
proton paramagnetic resonance. The measuring device 
enables determinations to be made of the value of the 
magnetic field over the effective operational range, 
to an accuracy of =1: 10°. 

Much attention has been devoted to obtaining the 
required field distribution in the magnet pole gap, and 


E,=6,066 Mev 





to a study of its topography. In particular, measure- 
ments were performed on the azimuthal symmetry of 
the magnetic field distribution, and the position of the 
median surface was also ascertained. An analysis of 
the measurements showed that azimuthal distortions 
and deviations of the median surface from the mid- 
plane within the effective volume remained within 
permissiblelimits. | To remove the fringe effect, we 
used special magnetic shims which functioned as mag- 
netic potential dividers [7]. These were a set of iron 
strips between which slabs of cardboard were sandwiched. 
The arrangement of these strips along the vertical was 
arrived at by experimental trial and error. Assemblies 
containing sets of the strips were positioned over the 
entire length along the inner and outer boundaries of 
the pole gap. The use of the shims enabled us to 
secure a somewhat improved field distribution over the 
entire volume. 

Determination of the topography made it possible 
to select the proper geometry of the diaphragms colli- 
mating the beam of particles, and to ascertain the 
maximum solid angle of the instrument (8 - 104 x 4 m), 

Investigations of the ion optical properties of the 
instrument enabled us to find the position of the focal 
plane, and to carry out a calibration of the spectrog- 
raph with respect to energies. The energy range of 
the a particles recorded simultaneously in the photo- 
graphic plates proved to be equal to AE/E = 10%. 

The half-width AHp/Hp of the line in that range did 
not exceed several hundredths of a percent. The disper- 
sion of the instrument was close to predicted values, 
being AE/Ax = 2.3- 10“ per mm, which corresponds 
to 1.2 kev/mm for an & particle emitted by Po*2®, 
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Fig. 2, Fundamental groups of & particles of Cm”, 
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One of the most important problems in nuclear 
power engineering is the reduction of weight going 
into shielding. The design of a minimum-weight 
shielding depends on many factors [1]. Without pro- 
ceeding into an analysis of all of these factors, let us 
merely indicate that the geometry of the radiation 
shielding barrier plays a role not to be underestimated 
in achieving optimum weight design. The tolerable 
radiation levels may be different at different points in 
space, and the shielding thickness may,in consequence 
of that,be a variable quantity, but one such that will 
assure the maintenance of a specified level of radiation 
at that particular point in space, 

The solution of the problem of the optimum radia- 
tion shielding geometry in the general form for radia- 
tion sources which are complex in composition and 
shape involves formidable mathematical difficulties. 
For isotropic sources of simple geometry (e.g. linear, 
disk, cylindrical-surface or spherical-surface geomet- 
ries), the solution of the variational problem may be 
used to derive an equation describing the optimum~ 
weight shielding geometry. 

Consider the problem of determining the geometry 
of a radiation shield of optimum weight for a linear 
source extending 2L in length and having an activity 
M (mg-equiv Ra) uniformly distributed along the length 
of the source (Fig. 1). 
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Fig. 1. Geometric representation for design 
purposes, 
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The dose rate at point A, without taking multiple 
scattering within the shield itself into account, is 
determined from the equation 


0 
, MK our \ MK -vy sec ¢ 
Ri’ dz—\ Li ° dq, (1) 
0 


P= \ or 


where J is the linear attenuation coefficient in the 
shielding material (for a narrow beam); M/ 2L is the 
specific activity of the source; K is a proportionality 
factor; 9p is arctan L/H. 

The weight of the shielding barrier is 


¢ 
G= arcs | 2ptydx = \ 2ptyH sec? pd, (2) 
0 


where p is the density of the shielding material, and 
t is the transverse dimension of the barrier. 

It is imperative to find the geometry of a shielding 
barrier, i.e., y = f (x) or r= f (Y) under the condition 
where the dose rate at point A is of some specified 
value. The geometry must be such that the weight of 
the radiation shielding will be at a minimum. This is 
a variational problem having a conditional extremum 
[2]. 

We now find the minimum of the functional 

0 
G= \ QetyH sec? @ dp 
0 
for the isoperimetric condition 
0 
FS enV 80 & dep — Pp, (3) 
0 


where Pp» is some specified value of the dose rate at 
point A. 

The solution of such problems will be aided by 
setting up the auxiliary functional 


0 


G** = \ (2ptytt sec? Kv sec *) dq, 
0 





where \ is the Lagrangian constant multiplier. The 
functional G** is studied on the minimum by means 
of the usual methods of the variational calculus, i.e. , 
by the solution of Euler*s equation 


Fy 5 Fy’ =0, 


where F stands for the integrand. 


In our case, the function is not dependent on y, 
and, accordingly, 


Fy = O, 
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see pe“"Usec?_9, (4) 
From (4), we may determine y: 
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To determine X, we substitute (5) into (3): 
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so that 
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Finally, we have 
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y= in | —THp, 





Tene | 
COS Do : 
r=In [at al : 


(6) 





LHP, COS Po 


We thus obtain the equation r = f ( ¢) describing 
the shape of the radiation shield having optimum-weight 
characteristics. The curve arrived at is seen in Fig. 2. 
Similar reasoning is used to solve the problem of a 
disk-shaped radiation source. 

The equation (6) derived for the geometry of the 
radiation shielding barrier fails to take into account 
multiple scattering within the shielding. The varia- 





Fig. 2, 
shield. 


Geometry of optimum-weight radiation 


tion in shielding thickness due to multiple scattering 
effects may be taken into account by what is known 
as the method of successive (iterated) approximations 
[3] for values of r found by means of (6). 

It must also be noted that the variation in the dose 
rate at the point in question due to scattered radiation 
(e.g., scattered in air) is not taken into account in the 
calculations. The contribution of scattered radiation 
to the dose rate for a radiation.shielding geometry 
such as described by (6) will obviously be larger than 
that for a radiation shield of constant thickness. 

The experiments which we performed using a 
linear Co™ source,1 meter in length showed that, with 
the shielding geometry calculated on the basis of (6), 
the measured dose rate agrees, within the limits of 
experimental error, with the dose rate specified for 
this design, and that the weight of the shielding barrier 
is 20% less than that designed on the basis of multiple 
attenuations, 
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In the work discussed below, an experimental 
material was produced which may prove of value in 
the design of shielding against y radiation, and in 
calculations of multiple scattering of y photons. 

The total radiation intensity at distances up to 
600 m removed from the source (r), and the angular 
distribution of radiation at distances of 150, 250, and 
400 m were derived in the course of the work. Mea- 
surements were carried out with-the source and detector 
placed at a height of 1.5 m above the earth. At 
some points, the measurements were conducted using 
a source supported at the top of a slender duralumin 
mast, 25 m in height. 
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Fig. 1. Counter sensitivity plotted against energy 
of y radiation. 











Fig. 2, Diagram of the experiment. 


For measurements at large distances from the 
source, a powerful gamma emitter was made of Au, 
having an activity of ~ 10,000 C, The total radiation 
intensity was measured with a Geiger counter of the 
type described in reference [1]. The sensitivity res- 
ponse of this counter as a function of y-radiation 
energy is shown in Fig. 1. Note that multiply scatte- 
red y photons from the Au" source, at energies of 
120-410 kev, are registered by the counter with 
uniform sensitivity. Each event registered in the 
counter corresponded to a y-ray energy flux of 8 + 
£0.4 Mev/cm? passing through the counter. Radia- 
tion at energies of 60-120 kev was registered by the 
counter with slightly enhanced sensitivity. However, 
the maximum increase in sensitivity did not exceed 
20%. Radiation at energies below 50 kev completely 
escaped detection by the counter, In referring to 
measurements of the radiation intensity in what follows, 
we accordingly have in mind only radiation Ey > 50 kev. 

The angular distribution was measured by a radia- 
tion detector assembly consisting of an array of four 
counters placed in series. Theeffective volume of the 
detector setup was 4.5 X 4.5X 0.5cm. This detector 
was placed in back of a thick leaden disk measuring 
21 cm in diameter, in a manner such that the centers 
of the disk and of the detector coincided on the same 
straight line with the center of the source. The 
arrangement of the disk, source, and detector may be 
seen in Fig. 2. It is clear that the leaden disk shields 
the detector from any radiation arriving at the point 
of measurement at angles smaller than ¢y. This angle 
¢ is read off from the straight line joining the source 
at the point of measurement. By varying the distance 
between disk and detector, it was possible to obtain 
information on the angular distribution of the radiation 
fur angles y from 0 to 90°. The results of the mea- 
surements performed are shown in Figs. 3 and 4. 

As we see from Fig. 3, when the source and 
detector are placed close to the ground, the intensity 
of the radiation at great distances from the source 
proved to be roughly half that established for a homo- 
geneous air medium. When the source was placed in 
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elevation above the gound, this difference was lessened, 
e.g., at a height of 25 m, the radiation intensities 

at the same distances differed by only 1.5 times. The 
data fit well with the considerations voiced in [4] as 

to the effects of having the source at a distance from 
ground level.on the intensity of radiation in air. 

Data on the angular distribution of scattered ra- 
diation are given in Fig. 4, where we have indicated 
the percentage contribution owing to y photons arriv- 
ing at the point at angles smaller than the specified 
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Fig. 3. Absorption coefficients and inten- 
sities of nonscattered radiation on the dis- 
tance between source and detector: a) ab- 
sorption coefficients K (at T=0°C and H= 
740 mm Hg), showing the number of times 
by which the intensity of radiation in the air 
above the ground is weakened as a result of 
interaction with airborne matter; b) the same 
for the source placed 25 meters above ground 
level; c) exponential law of attenuation of 
the intensity of primary radiation from the 
source, with scattered radiation effects neg- 
lected; d) law of attenuation of radiation 
intensity of emissions of Au ina homoge- 
neous air medium [2, 3]. 
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Fig. 4. Angular distribution of the scattered 
radiation. Distance from radiation source 
(in m); ©) 150; X) 250; Ol) 400. 


angle y. As we see, within the limits of the statis- 
tical scatter of the radiation counter counting rate 
(+ 10% for r = 400 m), the angular distribution of the 
radiation at substantial distances was independent of 
the distance from the source, i.e., the angular distri- 
bution represented by the curve plotted in Fig. 4, while 
stable over the first 100-150 m, varied very slightly 
as r was increased further. 

The authors take this opportunity to express their 
acknowledgement to O. I. Leipunskii,with whom 
the idea underlying the experiment originated ,and 
to V. A. Rogachkov, V. A. Shabashov, V. N. Rodi- 
onov, whose kind assistance was of great value in the 
work with the powerful y-ray emitter under field 
conditions, 
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The inelastic scattering of fast neutrons on nuclei 
is of considerable interest both for the study of the 
structure of the nucleus and for the pratical design of 
fast reactors, biological shielding in nuclear facilities, 
etc. 


The present paper reports studies of the energy 
levels in nuclei,based on measurements of the energy 
of y rays emitted in the inelastic scattering of neut- 
rons of energy Ey = 3 Mev. The energy spectra of 
gamma rays emitted from titanium, chromium, stron- 
tium, iodine, barium, tungsten, iridum, and iridium, 
and bismuth nuclei were measured by means of a 
scintillation y spectrometer using Nai(T1) crystal 
scintillators, 


Experimental arrangement. The D(d, n)He® 
reaction was utilized as a neutron source. A beam 
of deuterons accelerated to an energy of 200 kev 
impinged on a target of deuterium adsorbed on zir- 
conium. This neutron source, with a power of 107 
neutrons/ sec, was brought to approximately 2 m from 
the massive accelerator parts. 





A lead cone 36 cm in height with a base diameter 
of 5. 8 cm was used to shield the crystal from direct 
impingement of radiation emanating from the target. 
The crystal with accompanying photomultiplier tube 
and preamplifier, the test samples, and the cone 
were suspended from the ceiling of the room on 
caprone filaments 0.7 mm in diameter. 


The elements to be tested were prepared in the 
form of rings. The height of the he avy-element rings 
was 30 mm, and that of the light-element rings was 
40mm. The rings had an I, D. of 60 mm and an 
O. D. of 85 and 95 mm for heavy and light elements, 
respectively. While measurements were in progress, 
the rings were fitted around the crystal. 


The scintillation y spectrometer consisted of a 
single Nal(T1) crystal 40 mm in diameter and height, 
connected to a photomultiplier via a light-pipe of 
plexiglass 6 mm thick, a nonoverdriven amplifier and 
a 128 channel kicksorter with ferrite memory storage 
unit. The energy resolution of the spectrometer was 
10% for gammas emitted by Zn®, 


Monitoring for constancy of neutron flux was 
handled by a full-wave boron counter set at a distance 
of 2 m from the accelerator target, at right angles 
to the direction of deuteron flux. 

Results of the measurements, In order to find out 
what was the contribution of either y rays or of 
neutrons scattered from surrounding objects and 
scattering agents to the cyrstal count, measurements 
of the pulse spectrum for pulses due to a carbon scat- 
terer were taken after each run of measurements with 
the scatterer being tested. All of the results were 
reduced to a uniform monitor count, and the back- 
ground measured with the carbon scatterer was sub- 
tracted. However, in the low-energy region of the 
energy spectrum (down to 0.3 Mev), the background 
was not subtracted, since this would have been asso- 
ciated with difficulties arising in the scaling process, 
owing to nonuniform absorption coefficients of the 
y rays of those energies, for carbon and heavier 
elements being tested. The y-ray spectra in that 
energy region were therefore arrived at, for certain 
elements, without deduction of background. 

Figures 1-4 (see following pages) give the spectra 
of the y rays accompanying the inelastic scattering 
of 3-Mev neutrons on nuclei of strontium, barium, 
tungsten, and iridium. The Table presents energies of 
Y quanta as measured in the course of our work and 
as reported by other authors. 

In conclusion, the authors view it their welcome 
duty to thank A. K. Krasin for his constant interest 
in the work and his cooperation in setting up the 
experiment, A. N. Serbinov for his kind assistance, 
and I. I. Bondarenko and V. V. Stavinskii for their 
fruitful discussions. 
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The question of the effect of the space charge of 
a beam on the free oscillations of particles acquires 
much significance in relation to the problem of obtain- 
ing high currents in accelerator machines. This 
problem has been discussed in the literature [1, 2], but 
the screening effect of the vacuum chamber walls and 
of the electromagnet pieces has not been taken into 
due account. The work reported by the authors takes 
into account both these factors, with omly the horizon- 
tal chamber walls being considered in the problem, 
since the screening effect contributed by the vertical 
walls is substantially less, 

Let us visualize a beam of particles as an infinitely 
long cylinder of radius p. Now, if we assume the 
walls of the vacuum chamber to be ideally conducting 
and take the permeability of the iron to be infinitely 
large, we shall find the following expressions to hold, 
outside of the beam, for the vector and scalar poten- 
tials [3]; 
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J 
where A is the component of the vector potential 
along the beam axis; B = v/c; N is the number of 
particles in a unit volume; x and z are deviations from 
the beam axis in the horizontal and vertical direc- 
tions; o is the ratio of the vertical aperture of the 
chamber to the spacing h between the electromagnetic 
pole pieces. Assuming p/h << 1 and solving the 
Maxwell equations for the inner region of the beam, 
we get: 
xt? : 
6x =2neNx [ 1s (p/h)? | ; 
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6x,2,Hx,, are components of the electric and mag- 
netic fields. The calculation of those fields leads to 
the following equations for free oscillations; 


x" +k? [1—(n—6n,)] x=0, 
2” -|-k® [n+-8n,] z=0, 


where differentiation is carried out along the length 
of an arc of the equilibrium orbit 0; & is the curvature 
of the orbit; n is the index of the accelerator guide 
magnetic field; 
2 2 
— Farge (1B) EE (erm (1/02 2p) | 
(4) 
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(here the upper sign refers to Sn, and the lower sign 
refers to 5n,; E is the total energy of a particle). 
From (3), it is not difficult to derive the formula 
for the frequency variation of free oscillations [4, 5]. 
ihe 
dys.2= Seg | UBF 
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The quantity l¥,, z |? is the mean value of the 
square of the modulus of Floquet's function, which 
in our case is normalized such that the Wronskian 


¥,.,: V2.2 — Wr2Vx,2 = —2j 


where L is the perimeter of the orbit. If the operating 
point falls at the center of a cell bounded by reso- 
nance bands, then assuming that | Avy 2| = 4 
(i.e., is equal to the distance to the nearest reso- 
nance), we find the limitation on the beam current 
for a beam injected into the accelerator chamber, 
From (5), we see that the screening effect 

(second term within the brackets) yields a contribu- 
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Fig. 1. Spectrum of pulses excited by y rays accompanying inelas- 
tic scattering of 3-Mev neutrons in strontium. 
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Fig. 2. Spectrum of pulses excited by y rays accompanying inelas- 
tic scattering of 3-Mev neutrons in barium. 
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Fig. 3. Spectrum of pulses excited by y rays accompanying in- 
elastic scattering of 3-Mev neutrons in tungsten. 
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Fig. 4. Spectrum of pulses excited by y rays accompanying inelastic 
scattering of 3-Mev neutrons in iridium. 
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tion of the order of 10-20% in the nonrelativistic 

case, Contrariwise, in the relativistic case (B ~ 1), 

it may play a paramount role, leading to an additional 
limitation on the number of sectors accumulated as 
the particles accumulate in the accelerators. 
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The properties of carbon materials which are of 
practical engineering significance, in particular its 
thermal properties, are related to the atomic and mo- 
lecular structure of those materials, 

X-ray diffraction analysis data have shown the 
coexistence of amorphous and crystalline phases in 
the process of homogeneous graphitization of carbon 
(at 2000°-2800° C) [1]. This process must therefore 
be assigned to the category of orientation processes, 
which are similar to processes involved in the forma- 
tion of oriented polymers [2]. 

Homogeneous graphitization of carbon is character- 
ized by a continuous improvement in the spatial order- 
edness of the crystal lattice of graphite, through azi- 
muthal orientation [3] of parallel-layer elementary 
atomic lattices of carbon in groups (packets) formed 
prior to the crystallization stage. 

This letter reports results obtained in research 
on the specific heat, thermal diffusivity, and thermal 
conductivity, as well as x-ray powder determina- 
tions of the degree of graphitization of petroleum 
coke under isothermal conditions, for a range of diffe- 
rent temperatures and furnace times. Petroleum 
cracking coke of 1.405 g/cm* density with an ash 
content of 0.08% and a 5.13% yield of volatiles was 
used in the experiment. The dry weight of the sam- 
ples ranged 0, 50-0, 60 g/ cm*, and the grain compo- 
sition ranged from 0 to 0.5 mm. 

Graphitization was accomplished in a furnace 
with a graphite heating unit, in an atmosphere of 
nitrogen and argon. The maximum departure of the 
temperature from set point were + 25°C, The degree 
of graphitization was determined[4] on the basis of 
measurements of the interlayer spacing d (002) from the 
formula 


y= dmax —4 , 


dinax 4min 


where d is the measured interlayer spacing; di, = 

= 3.356 A is the interlayer spacing in the crystalline 
lattice of critically graphitized carbon; diya = 3.425 
A is the interlayer spacing of nongraphitized carbon. 
X-ray patterns of the graphitization products were 
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obtained by filtered radiation of copper and iron in 
cylindrical 70- and 43- millimeter diffraction cameras, 

The specific heat and the thermal diffusivity 
constant of the samples were determined at 20-35° C 
according to the method of the regular thermal regime 
[5] modified for the study of trace quantities of 
matter [6]. The results of parallel experiments differed 
by not more than 1%, 

On the basis of experimental data on the determi- 
nation of the specific heat c,, the thermal diffusivity 
a, and the weight by volume 6, the value of the 
thermal conductivity constant ) was calculated as 
XN=a cp (6 ). 

Figure 1 shows the degree of graphitization as 
a function of temperature and of the isothermal 
holding time, corresponding to the equation 
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Fig. 1. Interlayer spacings d, as a function of 
the time of isothermal graphitization of cracking 
coke. Temperature, °C; O) 2000; X) 2150; O) 
2300; A) 2420; 0) 2800; 
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where E = 95 kcal/ g-atom is the energy of activa- 
tion of the graphitization process; Ky = 10, 

The results of studies of the specific heat c,, 
thermal diffusivity @, and thermal conductivity A» 
of the samples as temperature functions and as func~- 
tions of the isothermal holding time in the furnace, 
as well as a function of the degree of graphitization 
y are plotted in Figs. 2-4, 

The regular patterns observed in the variation 
of the quantitites y, cp, 4, and A may be explained 
by the kinetics and the mechanism operating in the 
structural-chemical transformations, The original 
coke had the structure of a three-dimensional polymer 
in the form of an assemblage of relatively small 
planar monolayers of carbon interconnected by side- 
chain atomic groups containing, in addition to carbon, 
hydrogen, oxygen, and many other elements [3]. In 
the precrystallization stage of heat treatment (up to 
~ 1800° C), thermal destruction of the side-chain 
radicals took place, accompanied by a growth of 
carbon monolayers due to the coalescence of the latter 
and the formation of packets of parallel monolayers 
connected by side chains of carbon with monolayers 
of the neighboring packets. These processes were 
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Fig. 2, Measurements of the specific heat (0) and thermal 


diffusivity (O) as a function of the time of isothermal 
graphitization of cracking coke. Temperature, in °C; 
1) 2000; 2) 2300; 3) 2400; 4) 2600; 5) 2800. 
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accompanied by the escape of volatile products with 
the attendant loss of a large quantity of the terminal 
atomic groups, and led to an increase in the fraction 
of carbon in the monolayers being included into the 
structure of the three-dimensional polymeric carbon 
skeleton. The structural-chemical transformation in 
this stage of heat treatment found their reflection in 

a reduction in specific heat c, (reduction by ~ 50%) 
and an increase in the thermal diffusivity constant 

a (Fig. 3). During the graphitization stage proper, 
orientation processes occurred in the stacking of carbon 
monolayers into place in the crystal lattice of graphite, 
those processes being linked to the chemical processes 
involved in the destruction of the side carbon chains, 
as indicated by the high value of the energy of acti- 
vation of the graphitization process (~95 kcal/g-atom). 

A lowering by ~ 1/3 in the specific heat during 
the graphitization stage pointed to loss of two degrees 
of freedom in the overall vibrational spectrum of 
carbon atoms, 

The physical causal agent for that effect is rela- 
ted to the redistribution of the frequency of thermal 
lattice vibrations (as indicated earlier for lamellar 
structures [7] in the monolayers of carbon freed from 
the side valency chains during the process of destruc- 
tion of the side carbon chains, 

The linear behavior of the dependence of the 
specific heat Cp thermal diffusivity constant 4, and 
thermal conductivity on the degree of graphitiza- 
tion y (Fig. 4), if the initial period of graphitization 
(~ 0.1), during which stacking processes in the 
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thermal conductivity ( A) as functions of the degree of graphiti- 


zation. 


precrystallization stage are still possible, is excluded, 
points to structural and chemical transformations of the 
polymeric carbon backbone. 

Relatively small increases in the thermal conduc- 
tivity A, the basic changes . which occurred during 
the precrystallization stage, were noted during the 
graphitization stage. In order to obtain carbonaceous 
material exhibiting relatively high thermal conducti- 
vity combined with rather high mechanical strength 
properties, it was not necessary, therefore, to effect 
the complete graphitization of the carbon. Heat 
treatment of the original petroleum coke stock need 
be carried through no further than the initial stages 
of graphitization. In order to enhance heat stability, 
which is related to the thermal diffusivity, efforts 
must be made to reach critical graphitization of the 
carbon, Data obtained enable us to also explain the 
small degree of variability in the specific heat of 
neutron-irradiated graphite [8, 9]. When graphite is 
subjected to neutron bombardment, the three-dimen- 
sional ordering of the carbon, y, suffers considerably, 
and the interlayer spacing d between the periodic 
carbon monolayers increases, with accompanying 
swelling of the graphite [8]. 

The graphite monolayers in samples having a 
high degree of graphitization of the carbon [y ~ 1) 
are practically without side bonding (carbon chains), 
The internal changes in the monolayers in response 
to neutron bombardment, the azimuthal disorientation 
(and consequently the lowering of y and increase in 
d) are not associated with the transformations in the 
spatial polymeric carbon backbone, and are therefore 


without effect on changes in specific heat. We may 
add to this fact that the swelling of the graphite in 
itself testifies to the absence of side links between 
the monolayers of neighboring crystallites. 

The results of the research may be used to produce 
carbonaceous materials from petroleum coke having 
preengineered thermodynamic properties, 


LITERATURE CITED 


R, Franklin, Acta Cryst. 4, 253 (1951); R. 
Franklin, Proc. Roy. Soc. 209, 196 (1951). 

V. A. Kargin and G, L. Slonimskii, Uspekhi 
Khim. 24, 785 (1955), 

V. I. Kasatochkin, Izvest. Akad. Nauk SSSR, 
OTD 9, 32 (1951). 

V. I, Kasatochkin and A, T, Kavernov, Doklady 
Akad. Nauk SSSR 117, No. 5 (1957). 

G, M. Kondrat’ev, The Regular Thermodynamic 
Regime [in Russian] (Gostekhizdat, Moscow, 
1954). 

V. K. Zamoluev and V. I. Kasatochkin, Izvest. 
Akad. Nauk SSSR, OTN 11, 190 (1957). 

V. K. Zamoluev and V. I, Kasatochkin, Koks i 
Khimiya 6, 21 (1957). 

V. V. Tarasov, Doklady Akad. Nauk SSSR 46 
117 (1945); V. V. Tarasov, Trudy Inst. Krist. 
Akad, Nauk SSSR 10, 304 (1954), 

Hennig and Howe, Report 751, Geneva Confe- 
rence (1955) Vol. 8, 

W. De Sorbo and W. Tyler, J. Chem, Phys. 
26, 244 (1957). 
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In this contribution, an equation describing the 
transport of isotopes in separatory equipment is derived. 
No assumptions are made in the derivation as to the 
method of separation used, the design of the equip- 
ment, or the character of the fluid dynamic flows, 
which enables us to utilize the equation for the study 
of a wide variety of separatory devices, The equa- 
tion was derived by a method first used in the study 
of a thermal diffusion column [1, 2], and is based 
on an assumption that the variation in concentration 
over the transverse cross section of the equipment 
is small, which is a restraint well satisfied in most 
cases of isotope separation, 

A considerable number of the processes used for 
isotope separation are described by the equation 


t1—te= —B [Ve+-8c (1—c)], (1) 


where 7; _is the convection velocity of the light 
isotope; T is the convection velocity of the isotope 
mixture; c is the molar fraction of the light isotope 

in the isotope mixture. Coefficients B and 6 are, 
generally speaking, functions of the coordinates. Their 
value is dependent on the method of separation in 

use, For example, in the thermal diffusion process, 

B is the diffusion coefficient, 6 = a Vin T(a 

being the thermal diffusion constant, and T the tem- 
perature); in centrifugation 


6 ae an wr 
(w is the angular velocity, m, is the mass of a 
molecule of the ith isotope, and k is Boltzmann's 
constant). 

In what follows, we shall restrict ourselves to 
the steady-state two-dimensional problem, for the sake 
of simplicity. For process equipment possessing 
cylindrical symmetry, the equation is derived in simi- 
lar fashion. 

A diagram of the separatory process equipment 
is shown here. The flow of the mixture of isotopes, 
admitted to the bottom of the equipment and dis- 
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lan 


Diagram of separatory 
equipment. 


charging from the top, is equalto I, The correspond- 
ing flow of light isotope is lj. Then, obviously, 


a 
I= \ t,dz; 
0 
a 
= \ T%,d2. 
0 
Let us now write out (3.) for the components: 
dc 
tiz—txe= — Bi {5-4 8:0(1—e) ; (2 


Tiz—T:c= —B { 52-10. (1—e)}. (3) 


We introduce the notation 


x 
G(z,z)= ( t, dz; @(x,3)= —p{ 24 6c (1—o)} 
0 


Next, we take the integration of (3) with respect to 
X from 0 to a; 


a a 
1,—Ieq= gaz —\ G (2, 2) 32 ae, (4) 
0 0 





where Cq is the value of c(x,z) atx =@. We then 
eliminate 0c/@x from (4). For this, we integrate 
the equation of continuity of the light isotope, div 
T= 0: 
a x 
Tix = — \ TW, dx 

i (5) 
and transform (5), with the aid of (3) and the equa- 
tion div 7 = 0: 


ac c 8 
. c 
t= —5, \ g dx— -\ Tz 5g ot 
0 1) 


Oc 
+ ct, — 3 


Tx Ox dz . 
Substituting (6) into (2) and introducing the notation 


eet 
p (2, oes \ rP \ par 
0 (7) 


x 
+ \ tT: Z dz )—b,0 (1—c), 


du Tx 
a ve 


x 


u=\ te 


0 


We then integrate this equation, and keeping in mind 
the fact that dc/ dx = 1/%* du/ dx, we find 


x” 
° S54 \ 
dc B 


jie (8) 
a=5° 5 pre ° 


dx’ +p. 


Substituting (8) into (4) yields: 


dx! v | dz. 


When G(x, Z) may be represented in the form G(x, z) 
= X(x)Z(z), it is easily demonstrated, with the aid 
of the equation of continuity div 7 = 0, and the 
obvious condition T, = 0 at x= 0, that G(x, z,)/T 

is independent of x. This enables us to transform 
(9) to the form 


a 


a 
I,—Icg= \ pdz—e” 
0 


a -| “X dy 
x \ G (x, z) e dz. (10) 
0 


Put Ac = c(Z, 0)—c(z, 2). In accord with (2), 4c ™ 
~aé,. Inthe separation of isotopes, a6, << 1. This 
atone the penny of replacing the allie of c, 
de/ dz and d*c/ dz” acating into the integrands by 
Cq dceg/dz, andd 2ca/ dz”, respectively, retaining 
accuracy to terms of the order of (a6,)’, and to then 
remove them from under the sign of the integral 
— —— to x. In addition, terms proportional to 

d’ca/ dz" and a, d/d,[ca(1-ca)], also of the order 
of (a6, y » May be deacanied. Equations (9) and (10) 
acquire the form 


I,—Ica=H (2) eg (1—ca)—K (2) “8 . Gy 


If G(x, z), which must be found by solving the hydro- 
dynamic problem, is known, the task of calculating 
H(z) and K(z) reduces to quadratures. When (10) 
and the constraint 6, << 6, are valid, as is usually 
the case in isotope separation equipment, then 


t tsa { tar 
8 rs, 29) 


, 


=. ee (13) 
K (z)=e B i+) 8 dz. 


Equation (11) may be put to use in studies of 
both elementary separation equipment and of column- 
type devices with a complex configuration of hydro- 
dynamic flow patterns (cf., e.g. [3]). 
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A very familiar property of tributylphosphate 
(TBP) which sets it off from many other industrial 
solvents is its capacity to yield relatively stable chem- 
ical compounds with extractable components [1], 
which results in a competition between various substan- 
ces in the organic phase and, as a consequence, in the 
exclusion of some compounds in favor of others. 
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Uranyl] nitrate displayed a particularly vigorous 
ability to crowd out competitors. Owing to the high 
value of the partition ratio (Kj * 40 [2]), uranyl 
nitrate was capable, even at moderate concentrations, 
of displacing HNOs (Kj, * 0, 2) and other compounds 
from the organic phase, In this connection, at the 
bottom of the tower, where uranyl nitrate is absent 


ntration radient of HNO. 
amas moles iter . 
uj -02 -0f 0 Of 02 03 


\ 
\ 





1 





NN 


* 





| 








T 
i 


T 


Reextraction 





L 


———number of stage 


Exttaction 


"A 





























Cc 


P (2 "= 3g/ liter 
2\r"=4 mole/ liter 


7 


"= 110 g/liter 


p(t 
r"=3,2 mole/ liter 


1 








co 














oy [2 





z 
ai z™ % from feed 


Fig. 1. x-y plot of the extraction and nature of variation in concentration gradient; a) uranyl nitrate 
(1, 2, 3 equilibrium lines for acidities, xH = 2, 3, 4 moles/ liter, A, operating line, Il, critical 
point); b) HNO, (1. . .7. numbers of stages); c) variation in concentration gradient of HNO, along 


height of tower; d) plutonium. 





and the partition coefficients of HNOg and of trace 
elements present are comparatively large, those 
substances are extracted efficiently. At the top of the 
extraction tower, on the other hand, uranyl nitrate 
being present there in a large yield, forces out the 
remaining elements from the organic phase, i.e., 
reextraction of the latter takes place under the effect 
of uranium. An internal circulation appears, leading 
to a buildupof some compounds in the midportion of 
the extraction tower; the peak concentration may be 
several times in excess of input (cf. Fig. 2a). In this 
letter, we take up briefly the question of this internal 
circulation on the basis of x-y diagrams plotted in 
somewhat unusual form. The further solution of the 
problem (whether by graphic or analytic means) pre- 
sents no great difficulty. 

Consider first the distribution of bulk elements, 
and then the distribution of trace elements governed 
by the former. A graphic analysis of the extraction in 
multicomponent systems, where the partition coeffi- 
cient o depends on the concentration, of several 
different components (e.g. , oY on x and x ), may 
be carried out by the method of successive iterations, 


r'g/1 x % 


(approximations), plotting individual equilibrium curves 
for each plate. Since the equilibrium curves for 
uranium do not vary much asacidityis increased, the 
x-y diagram for uranyl nitrate has the usual form (Fig. 
la). Equilibrium curves may be based on data from 
the literature [2] (e.g., for 20% TBP in kerosene), the 
operating line may be based on assigned terminal 
points or one such point and the average ratio of flow 
of phases n = V/L.* The stages are plotted, to a 
first approximation, on the equilibrium curve corres- 
ponding to an initial acid concentration x9 (dashed 
line), and to a second approximation by taking the 
variations in acid concentration on each plate. The 
corresponding relationship between concentrations 
and stage number of tower height, at HETS = const and 
n = 4, 0,is seen in Figs. 2a, 2b. 

The x-y diagram of HNOs (Fig. 1b) differs in the 
plotting of individual equilibrium curves for each plate 
on which uranyl nitrate is to be found (usually involv- 
* In practice, sufficient accuracy may be secured by 
assuming n equal to the ratio of flows, at entrance to 
the extraction column, of the aqueous phase. 
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Fig. 2. Distribution of concentrations along the height of the extraction tower; 

a) design distribution of concentrations by stages, in the aqueous phase; 1) uranyl 
nitrate; 2) HNO; at number of stages N = 7(2a, N = 10; 2b, at N = 3); 3) pluto- 
nium (1a and 3a uranyl nitrate and plutonium at n = 3,85, 9 = 0.95); b) distri- 
bution along height of packed column: 1) uranyl nitrate, 2) HNO3; c) dependence 
of concentration of plutonium in discard on 7; d) diagram of tower in the form of 
a cascade of stages. 





mi two or three of the top plates, e.g., in Fig. 1 

= 110 g/ liter, x,U = 3 g/liter, xy * 0; the numer- 
ical symbols are given to correspond to Fig. 2d); for 
the remaining plates, the equilibrium curve for the 
pure acid is used. In the extraction process, a tran- 
sition occurs from the equilibrium curves P, to P, and 
Ps, while the total equilibrium curve acquires the 
form of a loop (Fig. 1b). A peculiar form is taken 
by the operating line. It rotates through 180° at the 
point of transition from extraction to back-extraction, 
at that point, the concentration gradient $= y, —y 
and the transverse flow of material j 1 = kd change 
signs(Fig. 1c), while the acid concentration runs 
through a maximum. The operating line is not bounded 
by the inflow and exit points (Fig. 1b); furthermore, 
those points are not fully specified: only the initial 
es anecoien of the acid in the aqueous (x9 ) and organic 
cy 4 1) phases and the ratio of flows n = ee (i,e., 
the slope of the operating line tan & = \ = 1/n) are 
given, while the correspondi ng values a the second 
set of coardinates (xj and y;’) are unknown, The 
position of the operating line is therefore chosen by 
assigning some value of y; or x,, and dividing the 
series of trial plots, in order to obtain a specified 
value of yN+1 oF x f for each "anand number of 
plates considered [ xy and yy". may be found to an 
approximation by assuming that the partition coeffi- 
cients a" (x"’, x~) are dependent solely on xU, and 
using (1) and (2)]. 

The maximum buildup of acid, xH ox: is bounded 
by the intersection of the operating line with the 
equilibrium curve, and increases as n = V/L; atn = 
= 2, x4 = 3 moles/ liter, and at n = 4, xH x =4 
moles/ liter, Since Xeoias is bounded, then, as the 
number of plates are increased, a plateau appears in 
the concentrations, as seen in Fig. 2a. As is evident 
from inspection of Figs. 2a and 2b, the graphs plotted 
may be extended, to some degree of approximation, 
to packed columns.f However, the HETS in the 
uranium zone is somewhat higher than in the uranium- 
free region, 

Extraction of trace elements. The partition 
coefficient of element M is determined by the concen- 
tration of bulk elements (uranyl nitrate and HNOs), 
aM - of! cx U x), Accordingly, it is not difficult 
to construct the ee curve for each plate: 
= aM a xH) xM x; Where i is the number 





Yequilb. i 


of the plate, and xP and 3 are as determined above, 


All of the constructions may be carried out for NHO, 
graphically (as an example, we have in Fig. 1d the 
x-y diagram for Pu (IV); %,, is taken from data 


reported in [3]). However, in some cases it is simpler 
to perform an analytical calculation, Using the equa- 
tions for the operating line Yi = Ynga + A O24 - Xn) 

and for equilibrium y; ~ %j X; 48 recurrent expressions, 
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we find the concentrations of the element of interest 
in the discard (X,)) and in the organic phase at exit 


(ys): 


— Tot nye (1—Ta; + Z1la)) . 1 
1+ 2Ma, () 





yv.i1la; + Ax,d Ta; 


ee 





(2) 


where A = 1/n; CA = Oj/A =n; TOj= Oj, d,... 
ay jC 1a; = Oly + Oly Ol, Hes ot (OQ, . . Xp). The 
degree of extraction of the element, T, at Yue = 0 
is 


_ 4 (1-+-0,-++a,03+.. » Ogi . . A, ‘ 
ag 1+ Ma; (3) 





Analyzing the x-y diagrams and (1)-(3), w¢ may 
break the trace elements down into two groups, When 
the partition coefficient of the element is comparatively 
small in the absence of uranium (% < 1/n, a < 1), 
the buildup (as in the case of HNOs) will be bounded 
by the intersection of the operating line with the 
equilibrium curve, and the degree of buildup Xmax/Xo 5 
=/(1-O); when O&» << 1, then that element will be 
extracted poorly, and will pass almost in its entirety, 
for pr actical purposes, into the discard. In that case, 
we find, from (3), T =aT << 1, xy = x9 (1—- @ ) 
® Xp, the separation factor K = 1/T 8 1/a,.% Elements 
for which Gp» = on > 1 may be considered readily 
extractable (TC “1, xy ™ 0), In that case, the build- 
up is not bounded and will be dependent on the number 
of plates on which uranyl] nitrate is found in amounts 
sufficient to secure reextraction, i,e., for the cons- 
traint 3, < 1 to be observed. If this constraint 
is not satisfied on at least one plate (n > 1/a,,:,), 
then there will be no buildup of the element at all, 
the element will be extracted monotonically [e.g., 
Pu (IV)] atn > 5, The degree of internal buildup 
in a readily extractable element comprises (provided 
the uranium barrier is acting on m plates); 


_Xmax_ . : 
Imax Xo ( %4n)(%n). . .(%,n) 





tT Measurements performed by S. Medvedev and E, 
Rodionov. 


t > It is clear that for similar trace elements, satura- 
tion by uranium of the exiting organic phase entails 


a decrease in oO, and will lead to improved separa- 
tion. 





Thus, for Pu(IV), m=1, a * 0,18, and the degree In schematizing the distribution and visualizing 
of buildup q,,,, = 1/(4° 0.18) = 1.4 (Fig. 2a). The & = % = const in the uranium-free region, while a = 
buildup will be all the greater, the smaller the value =, = const in the uranium-loaded region, we find 
of 4, i,e., the greater the saturation of TBP with 
uranium at exit from the extraction column. Ig q aw 
The buildup will rise sharply with an expansion of N=No+Ny+AN= “ee 
the uranium zone, e.g.,asn decreases (curves la 4N leq wigs (4) 
and 3a in Fig. 2a), etc., right up to the critical UT Igagn ’ 
region corresponding to the intersection of the operat- 
ing line with the equilibrium curve for uranium, At where Np is the number of plates required to extract 
the transition beyond the critical region, the inade- the element in the absence of uranium; Nyy is the 
quacy of the volume of solvent leads to discard of the number of plates in the uranium zone; AN is the 
basic components. number of plates required to reduce the concentration 
Accordingly, the degree of approximation to the of the element from the maximum to the initial 
critical regime n = 1 (n, X9-, To) = yi3/yP (x97) is concentration; q = 1/(1-T’) is the degree of separa- 
_an important characteristic of the tower; on 9 depends tion with respect to the particular element. 
on the buildup and the losses in valuable components The increase in the required number of plates 
(see Fig. 2c), and the required accuracy in control. associated with internal circulation is rather high, 
In the presence of a uranium barrier and buildup, e. g., for Pu(IV) atl = 99.99%, Ny = 4.5 and Ny + 
additional plates are required for the complete extrac- + AN = 2(a 40% increase). In addition, there is also 
tion of valuable components (such as plutonium). a certain reserve margin of plates AN, required to 
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centrations by stages in extraction with scrubbing operation (1. uranyl nitrate, 
2. HNOg, 3, Pu); b) distribution of concentrations by stages in the case of reco- 


very in aqueous phase (1. uranyl nitrate, 2, HNOs); c) diagram of extraction and 
scrubbing tower. 





avoid waste during excursions in the direction of the 
critical mode of operation, in the course of process 
control (Fig. 2a). 

More complicated cases, such as extraction with 
scrubbing, or solvent recovery (Fig. 3), are approached in 
similar fashion. In the first case, the concentrations of 
the elements in the aqueous phase at the feed point of 
the basic solution vary stepwise in line with the displace- 
ment equation L"X" 4) , LoXqq = (L" + L)Xo; the concen- 
tration in the organic phase varies smoothly, on the other 
hand, The concentration of bulk elements are found by 
graphical means, and the concentrations of trace elements 
are found by either graphical or analytical means (in the 
latter case, the solution of the system of equations for the 
scrubbing and extraction sections of the tower is required, 
using the displacement equations as the conditions which 
must be obeyed, and assuming y, = Y "N+p: By assum- 
ing yyy4 1 = 9, we find the concentrations of the discharg- 
ing streams; 
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where 


A=Nla’, B=XIia’, C= XMla, D=(1+B-+4 AC), 


For most compounds, x'y = 0, and the equations 
required for the calculations are simplified. In parti- 


cular, the concentration of the elements remaining 


in the organic phase after the scrubbing operation will 
be 


Ila’ Sa (1-4, )*22 


1 =, A/(A-+ B) = Se 
n=nAl+ = — sna’ + le’ Sila (7) 





From (7) and (3), we find that the separation 
factor in scrubbing, K = y,/ y;' = 1/1 a =1/( a,'n) 
(Q2'n). . . The scrubbing step therefore effectively 
removes only those elements having rather small parti- 
tion coefficients (if «;" = 1/n), Scrubbing, by setting 
up additional circulation of the uranyl nitrate, brings 
the extraction section of the transfer unit closer to the 
critical mode of operation; the ratio of process flows 
of phases no = V/ Ly must be increased slightly compared 
to extraction without scrubbing: np = f (n", 7). 

The performance of a tower for solvent recovery 
of uranium and plutonium (Fig. 3b) is viewed similar- 
ly. It is interesting that, at large flow ratios (n> 10), 
a sharp increase in acid concentration on the last plates 
leads to the uranium being displaced by the acid, i.e., 
to internal circulation of the uranyl nitrate, 
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In the solution of a number of problems encoun- 
tered in the design of powerful y facilities for radia- 
tion processing of objects and in the design ot radia- 
tion shielding, the dose field of sources of different 
(and sometimes complex) configurations must be 
determined. 

Exact solutions are available only for the simplest 
types of extended sources [1, 2],some of which are of 
no practical interest. For several other cases of prac- 
tical interest, nomograms have been developed or 
auxiliary coefficients have been determined through 
numerical integration of the equations derived for the 
dose rate [3-5]. 

Up to the present, there is no generally applica- 
ble method of calculating the dose field of extended 
sources which would make it possible to carry out 
the analysis in any desired case (in a homogeneous 
medium or in its absence) to some exactly determined 
appr oximation. 

The method proposed for calculating the dose field 
of extended sources makes it possible to dispense with 
numerical integration, and to obtain approximate 
expressions, with a desired degree of accuracy, in 
elementary or tabulated functions, We shall consider 
this method in the light of exarnples of calculations 
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Fig. 1. Transition from a rectangular planar source 
LS with constant surface activity q to a fictitious 
linear source O d = J with variable specific linear 


activity over different segments: Oy, =L.. qy3 L S 
<yp =S...q9; S Sys = br Fie 
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of source configurations of practical interest which 
have no generalized solutions up to the present writing. 

Infinitely thin rectangular plane of finite dimen- 
sions, In order to calculate the dose rate at any point 
P found over a planar rectangular source, outlined 
in Fig. 1 by a dashed line, it is sufficient to find the 
dose rate due to a rectangular segment Oabc of the 
source, having dimensions L and S for any relationship 
of the sides, 

Let the distance of the point p from the radiating 
plane be R, and let the specific surface activity be 
q, the width of segment Oabc be S= nL(n 21), All 
elements of the plane segment which are found at some 
equal distance p from the point p make the same 
contribution at that point to the dose rate, and may 
be replaced by a source dl found at the distance p, 
but exhibiting an activity equal to the sum of the 
activities of the replaced elements dS, 

In other words, the dose rate which is established 
at point p by the segment Oabc may also be derived 
from. a linear source Od of length J, having an activity 
distributed in obedience to some law along the length 
of the source (cf, Fig. 2). The total activity of the 
linear source J will then be equal to the total activity 
of the segment Oabc. 

We shall find out how the activity is distributed 
along the length of the fictitious linear source J. 

The elements of the plane which are found at 
some uniform distance from the point p lie on arcs of 
circles of radius r , described from a point O(Fig. 1). 
The linear specific activity of the source J at each 
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Fig. 2. Distribution of specific linear activity along fic- 
titious linear source, for different ratios of the sides 
S/L= n of the planar rectangular source replaced. 
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point of the source is proportional to the length of the 
corresponding part of the arc sweeping over the area 
of the segment Oabc of the planar source. 

At first, the linear specific activity increases 
over the interval of the linear source OSy SL, in 
obedience to the rule 


n=159- (1) 


The linear specific activity over the segment 
L Sy SS then drops off according to the law 


q2=Qy arctg VP—D * (2) 


From the point y = S to the point y = is? +L? =J, the 
linear specific activity drops off to zero, obeying the 
law 


L ‘inti 
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Figure 2 depicts the distribution of the specific 
linear activity along the fictitious linear source, whose 
action at point p is equivalent to the action of a 
segment Oabc of the planar source equivalent. The 
diagram shows how the distribution of the activity 
varies as the ratio S/L =n is increased. 

For dz and qg, we may assume a linear distribu- 
tion pattern to hold. We then obtain 


m 
"=95 y=ay for OC y<L. 


In place of qz and qg, let us take 
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The mean error in a linear approximation follow- 
ing (4) is 2-5% depending onthe source parameters 
and source geometry. 
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The dose rate from the segment Oabc of the planar 
finite source in air will be 
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A check was carried out on (5) by comparing 
results with a nomogram generalizing the results of 
analysis by numerical integration, which was in turn 
performed for several values of the sides ratio of the 
segment Oabc [3]. Agreement was found to be with- 
in the experimental limits of error in the determina- 
tion of the values on the nomogram for all values of 
R/S and L/S for which the nomogram was constructed. 

When the planar source was placed in an absorb- 
ing medium, with multiple scattering taken into 
account, the dose rate due to segment Oabc was deter- 
mined from the equation 
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Figure 3 shows the variation in radiation dose rate 
at a distance R = 50 cm along the axes of a rectangu- 
lar planar source 200 X 100 cm, having a specific 
surface activity a = 1 mg-equiv/ cm? in air and in an 
absorbing medium (water), The presence of an absorb- 
ing and scattering medium not only lowers the dose 
rate, but also equalizes the dose field at the same 
time over the planar source. 

Planar source in elliptic form. A calculation of 
the dose rate in air and in an absorbing medium for a 
radiation source in the form of a disk, for a point 
located above the center of the disk, may be found in 
the literature [2]. We now derive the equation for the 
more general case, where the source has the form of 
an ellipse with any assigned eccentricity. Let there 
be given an ellipse with semimajor axis a and semi- 
minor axis b (see Fig. 4), The distribution pattern of 








the activity along the fictitious linear source over the 
intervalO= p Sbis 


q1 = 2ngp (1) 
(where q is the specific activity associated with a 


unit area of the ellipse). 
For the interval b = p =a, we get 


q2=4paq, 


¥y 6 f/f. 
a=arc tan’ =arc tan| ° 6 
a? — p2 
m d2=40q arctan! = = } (8) 


Figure 5 shows the distribution of activity along a 
fictitious linear source of length a = 1 for the three 
cases; a= b, a= 2b anda=10b, Whena # b, a linear 
approximation is possible in the form 


where 


(9) 


The total activity of the linear source then 
remains equal to the total activity of the ellipse. 
The distribution of activity at a remote part of the 
source is the only activity distribution to vary some- 


what. Utilizing (7) and (9), we obtain, for a source 
in air 
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Fig. 3. Distribution of dose rate over a planar 
rectangular source (Co™) in air and in an absorb- 
ing medium, 


For a source placed in an absorbing medium 
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Linear sources, The method of calculations pro- 
posed here may be used to find the equivalent sources 
which are placed differently with respect to the point 
where measurements are taken. Let there be a given 
linear source L with a constant specific activity q over 
its entire length. We shall now find the linear source 
L, placed perpendicular to the first (Fig. 6a, b), which 
is equivalent to it with respect to the dose rate estab- 
lished at the point p. The linear specific activity 
q, of this new source is a variable. 

From Figure 6, we have 
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Fig. 4. Change from plane elliptic source with 
with a constant specific activity to a fictitious 
linear source, with the variable distribution of 
activity along the source. 
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Fig. 5. Distribution of specifit linear activity 
along a fictitious linear source for different ratios 
b/a=m(m =1). 
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Fig. 6. Transition from a linear source with 
uniform activity distribution to a linear source 
with irregular activity distribution, the second 
placed perpendicular to the first: a) transi- 
tion from a source placed along the x axis to 
a source placed along the y axis; b) transition 
from a source placed along the y axis to a 
source placed along the x axis, 
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Equation (12) yields the distribution pattern for 
activity along the length of source J. 

Let there now be given a linear source of length 
L laid out along the x axis, and having a uniform 
activity distribution pattem. 

The source of length 2 equivalent to it, and 
placed along the y axis, will exhibit a specific linear 
activity obeying a distribution law 


"n=4 (13) 
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Conclusion, The method proposed makes it 
possible to find approximate equations for the dose rate 
due to planar sources of different configurations, and 
the equivalent linear sources, with respect to dose rate, 
at a given point, and placed in different arrangements 
with respect to the first source. It would apparently 
be of some interest to extend this method to apply 
to three-dimensional sources, 
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News of Science and Technology 


SIXTH SESSION OF THE LEARNED COUNCIL OF 
THE JOINT INSTITUTE OF NUCLEAR STUDIES 


M. Lebedenko 


The sixth session of the Learned Council of the 
Joint Institute. for Nuclear Studies met at Dubna from 
May 27 through June 2, 1959. Leading specialists 
in the field of nuclear physics came from 12 socialist 
countries to participate in the work of the session. 

The delegations of the member nations of the 
Institute were headed by Petrik Pilika (Albania), 
Georgi Nadzhakov (Bulgaria), Albert Konya (Hungary), 
Le Van Thiem (Viet Nam), Heinz Barwich (German 
Democratic Republic), Tsien San-tsiang (China), Kim 
Hung Bong (Korea), Sodnom Namsrain (Mongolia), 
Andrzej Soltan (Poland), Horia Hulubei (Rumania), 
Cestmir Simane (Czechoslovakia), V. I. Veksler, 

N. N. Bogolyubov, and V. P. Dzhelepov (USSR). The 
chair was shared in succession by D, I. Blokhintsev, 
Director of the Joint Institute, and Vice-Directors 
Wang Hang-chang and E, Dzhakov, 

Blokhintsev reported to the Council on the most 
significant research efforts brought to completion in 
the Institute's laboratories during the past half-year, 
Some of this work, in line with the decisions reached 
at the session, will be reported at the international 
conferences scheduled to meet at Kiev (on high-energy 
physics) and Geneva (on particle- accelerator engineer- 
ing and related experimental techniques). Wang Hang- 
chang delivered a report on the development of inter- 
national liaisons between the Institute and scientists 
in other lands, 

The opening of the sixth session of the Learned 
Council coincided with the introduction into service 
of the new and splendidly equipped building housing 
the Laboratory of Theoretical Physics, where indeed 
the session was convened for a good part of the time. 
Members of the Learned Council heard reports on 
particular experiments of highest interest performed 
by a scientific staff of the Institute which is drawn 
from many nationalities, Some of those experiments 
produced new data on the structure of nucleons and 


on the production of mesons and antiparticles. Interest- 
ing experimental data were obtained using the proton 
synchrotron, in the work of a team of physicists headed 
by Wang Hang-chang and I. V. Chuvilo. Tens of 

of thousands of plates which they took will be useful 

in the further study of phenomena observed in inter- 
actions between pions and nucleons, and in the gene- 
ration of "strange particles." New phenomena in 

the field of the decay of ™-zero mesons was reported 
by a team of experimental researchers working in the 
Laboratory of Nuclear Problems. 

The members of the Learned Council noted that 
many of the research efforts completed at the Labora- 
tory of Theoretical Physics were in close corresponr 
dence with experimental observation. 

In the, resolution adopted by the sixth session of 
the Learned Council, it was noted that the achievement 
of round-the-clock and reliable operation of the proton 
synchrotron under the conditions for which it was 
designed is a signal accomplishment to be credited to 
the staff of the High-Energies Laboratory, and that it 
has made it possible to develop experimental work 
on a broad scope within a short period. The Learned 
Council singled out as one of the most noteworthy 
achievements in the field of advanced accelerator 
design the development and commissioning, at the 
Laboratory of Nuclear Problems, of an accelerator 
incorporating spatial variation of the magnetic field.* 

The sixth session of the Learned Council pointed 
out that the Joint Institute, which was founded only 
three years ago, has already become a powerful 
research center occuping one of the foremost places 
in the vanguard of world science, 


* D. P. Vasilevskaya et al., Atomnaya Energiya 
6, 657(1959) [Original Russian pagination. See 
C. B. translation]. 


FOURTH SCIENTIFIC AND TECHNICAL CONFERENCE 


ON NUCLEAR ELECTRONICS 


The IV Scientific and Technical Conference on 
Nuclear Radioelectronics, organized by the Main 
Control Board on the Uses of Atomic Energy, met in 
Moscow on April 20-25, 1959. Over 800 specialists 
from scientific research and industrial organizations 
of the Soviet Union participated in the conference, 


Over 100 reports and communications were heard, 
covering research work in the field of design of sophis- 
ticated electronic equipment for experimental nuclear 
physics and nuclear engineering, and also covering 
instrumentation already designed and made available 
by industry. 
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The conference was opened by P, S, Savitskii. 

Panel sessions were convened on particular topics 
such as nuclear radiation detectors, pulse-height 
analyzers and time analyzers, pulse equipment, auto- 
mation of measurments and data processing, spectro- 
meters, voltage and current amplifiers, and voltage 
and current stabilizers, 

The major attention of the panel session on nuclear 
radiation detectors was devoted to scintillation techni- 
ques in radiation research, and to the design of appro- 
priate instrumentation. In particular, questions touch- 
ing upon the design of new electronic phototube mul- 
tipliers for scintillation spectrometers and for record- 
ing pulses in the decinanosecond (10° 1 sec) range 
were discussed, as were problems associated with the 
design of detectors of ionizing radiations with high 
efficiency in recording and high resolution, including 
plastic, gaseous, liquid, and hybrid scintillators. 

Several reports dealt with research on the para- 
meters of scintillators and phototube multipliers (sta- 
bility and pulse-height discrimination in photomultip- 
liers), the effect of crystal parameters on pulse-height 
discriminations, the choice of crystal-photomultiplier 
tube systems for a scintillation spectrometer, etc. A 
special subpanel session took up questions regarding 
fast: neutron spectrometers and neutron detectors. 

The panel session on pulse-height analyzers and 
time analyzers heard papers on the design of multi- 
channel amplitude and time analyzers and their 
subassemblies and associated circuitry (readout instu- 
ments, memory storage, etc.), methods for decreasing 
resolving time in pulse-height analyzers, etc. 

A significant number of the reports were devoted 
to the development and design of multichannel systems 
( a 1000-channel amplitude-time analyzer, a 1024— 
channel time analyzer with magnetic-drum memory, 

a 50-channel analyzer, etc.). 

The panel session on automation of measurements 
and data processing heard reports on the extraction 
of data from multichannel analyzers with the aid of 
automatic printout devices, on the design of output 
digital printout devices for multichannel analyzers, 
automatic facilities for measuring activity in prepa- 
rations, etc. One paper was devoted to an arrangement 
for semiautomatic processing of photographic plates 
of tracks in bubble chambers, a design worked out at 
the Institute of Theoretical and Experimental Physics 
of the Academy of Sciences of the USSR, and a design 
project for an automatic arrangement for the same 
purposes, 

The panel on pulsed equipment devoted its sessions 
to papers dealing with the design of circuits and 
assembly units in pulsed equipment. 

The spectrometer panel considered reports dealing 
with several developed spectrometer facilities. 

The panel on voltage and current amplifiers and 
voltage and current stabilizers heard papers on stabilized 
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power supplies which had been developed (using 
transistors), electrometer amplifiers in rod-lamps, etc, 

The concluding plenary session was devoted to 
problems such as the radiation stability of semiconduc~ 
tors and standard radioactive sources for the monitoring 
and calibration of radiometric equipment. 

The participants in the conference adopted a reso- 
lution envisaging the development of research and 
experimental work on the design of electronic equip- 
ment for nuclear physics studies, . 

An exposition on display during the conference 
demonstrated samples of electronic equipment designed 
at research establishments and produced by domestic 
industry, and also several specimens of foreign equip- 
ment, 

On display were pulse-height analyzers, time 
analyzers, scintillation spectrometers, scalers, coin- 
cidence circuits, phototube multipliers, scintillators, 
and other equipment. 

Below, we present a brief rundown of the charac- 
teristics of some of the specimens on display. 

Automatic AZ-1 single-channel self-recording ana- 
lyzer for automatic measurements of the energy spectra 
of radioactive radiation, and for automatic measurements 
of a-, B-, y-radiation in large quantity from standard 
smaples (targets), The instrument records pulses of ne- 
gative polarity with a rise time of from 0.01 to 1 psec. 
The range of input pulses is from 0,01 to 25 v. The 
window has a width of 2.4. Irregularities in window 
width are ~ 1%, The resolving time is 5 psec. The ana- 
lyzer is equipped with a printing recorder. The time for 
recording a single result is 2.5 sec, The analyzer may 
be operated in two modes; the number of pulses per 
exposure may be measured, or the time corresponding to 
a given number of pulses may be measured. An auto- 
matic sample shift device is built to handle 40 samples. 

Automatic AZA-2Ts single-channel pulse-height 
analyzer for measuring the spectra of radioactive 
radiation and for the automatic measurement of radia- 
tion intensity (assembled with an automatic sample 
changer capable of accomodating 50 samples). Re- 
cording of the measurement data is handled by a 
printing recorder. The amplitudes of input pulses 
range from 0,01 to 40 volts (11 range levels). Win- 
dow width is 2,4, Resolving time is 5 usec, Time 
of exposure lasts from 0.4 to 1,000 sec. There is 
a high-voltage power supply unit for the photomultip- 
lier part, The instrument is made in the form of a 
stand, Power intake is 230 va, and the instrument 
weighs 75 kg. 

Differential single-channel ADOI-3 analyzer 
with rate metert for measurements of energy spectra. 

A scintillation detector is used in the instrument. The 
resolving time is 8 sec. The amplitudes of input 
pulses range from 0, 25 to 25 v (11 ranges). Window 
widths are 1, 2, 5. The thresholds vary stepwise from 

1 to 10%, 




















AIMA-2- pulse-height analyzer with 50 or 80 
channels, The capacity of each channel is 26 pulses 
in operation using 50 channels, and Qi pulses in 
operation using 80 channels, Resolving time is 1 .3 
Hsec, The input pulse amplitudes range 0-50 v. 
Peak duty factor of the analyzer is 5000 pulses per 
sec. Memory storage of the number of pulses in each 
channel is handled by a magnetic drum. The mea- 


surement data are displayed on a dual-beam screen, 
in binary and linear scales, There exists the possibi- 


lity of analyzing only those pulses which are time- 
coincident with the control pulse, Resolving time for 
the coincidence circuit is 1.5 psec. 

BMA-50-channel pulse-height analyzer for auto- 
matic measuring of energy spectra of elementary 
particles and other products of nuclear interactions 
(Fig. 1). Channel widths are 0.2, 0,3, 0.5, 1.0, and 
2.0 vz; Range of pulse amplitudes 0,01-100 v. 
Resolving time 0,182 psec: Capacity of each channel 
16, 000 pulses: Maximum counting rate 50, 000 pulses/ 
sect Range of durations of input pulses 0.1-1000 usec, 
The analyzer is equipped with a rate meter, dead-time 
indicator, monitoring pulse generator for calibration 
purposes, high-voltage rectifier for power delivery 
to the phototube multipliers, and a detachable cathode- 
follower circuit, as well as a preamplifier unit for 
reducing channel width and threshold by 20 times. 
Provisions are made for automatic switching of the 
analyzer to handle 8000 or 16, 000 pulses in either 
channel, as desired, The analyzer circuitry incor- 
porates about 1500 cold-cathode tubes and about 200 
thermionic tubes, Analyzer operation may be gated 
by an external selector or without it. 

AI-100 100-channel pulse-height analyzer using a 
ferrite-core memory for determining the amplitude 
spectrum of pulses distributed statistically in time 
(see Fig. 2.), Channel width instability is 4 1.5%; 
Channel width 1 v; Capacity of each channel 15, 999 
or 65,535 pulses; Range of input pulse duration 1-100 
Hsec; Average counting rate 10, 000 pulses/ sec? 
Average dead-time 70 sec. The analyzer contains 
an amplifier (for amplifying pulses of 0.01-1v), a 
rate meter with measurement ranges of 1000; 3000; 

10, 000; and 30, 000 pulses per sec with an accuracy 
of +10%, When any one channel is fully loaded, 
the analyzer is blocked; pulses beyond the hundredth 
channel are not recorded, The counting system uses 
a binary-decimal arrangement, 

1024-channel time analyzer for operation in a 
neutron selector system over a broad range of energies, 
with high-intensity pulsed neutron flux. The analyzer 
is capable of handling an input duty cycle of up to 
10, 000 pulses per sec. Channel-width steps are 1, 2, 
4, 8, 32 psec. All the channels can be broken down 
into four groups (of 256 channels each), Each group 
may be assigned a different channel width. To secure 
a dead time of 1 psec, a special intermediate memory 

















unit is introduced into the analyzer circuitry, to keep 
counting losses at peak duty factor from exceeding 
1.7%, The basic memory storage in the analyzer is 
based on magnetic cores, The maximum capacity of 
each channel is 2" pulses, Measurement data may be 
presented on a linear scale in a CRT display either 
during the measurements or afterward. Definitive 
measurement data are printed out. The data are 
punched out on paper tape, with channel number 
indicated. Printing speed is 12 characters per second. 
Numerical values in binary coding may be obtained 
by reference to the display devices of the adding cir- 
cuit in the basic memory unit, 

GSP table-model y spectrometer for measuring 
radioactive radiation spectra, The spectrometer 
incorporates a pulsed amplifier, a single-channel 
differential and integral analyzer, scaler, rate meter, 
monitoring integral discriminator, and a stabilized 








Fig. 1. 50-channel pulse-height analyzer, BMA-50, 
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positive voltage source rated at 0,7-2 kev. Input 
pulse heights range 0,01 to 25 v. Resolving time of 
the spectrometer is 5 sec, and the capacity is 16, 000 
pulses, 

PS-5 scaler for counting pulses with peak pulse 
repetition rate of 80,000 pulses per sec (Fig. 3). The 
scaler assembly includes a scaling unit, high-voltage 
rectifier, detachable cathode-follower circuit, gas- 
filled counters section and photoelectronic multiplier 
unit. The resolving time for binary pulses is 10sec, 
and channel capacity 10% pulses. The scaler is opera- 
ted by gating with pulses of 0, 5-10 sec duration and 
3 to 100 voltsamplitude. Ranges of measurement. 
are 1000, 10,000, and 100, 000 pulses per second, 
The instrument reads to an accuracy of + 10%, Data 
are read by neon glow-tube and dekatron indication. 
An electronic second-timer allows time intervals of 
up to 1000 sec to be read off. Counting may be halted 
and the electronic second-timer inserted after accu- 
mulating 10°, 10°, or 10° pulses, or after a time 
interval of 1, 10, 100, or 1000 sec has elapsed. The 
stabilized high-voltage rectifier provides a voltage 





Fig. 2, 100-channel pulse-height analyzer 
with ferrite-core memory units, AI-100. 


of 400 to 2000 v across the output at a current of up 

to 3 ma. In response to a variation of + 10% in network 

voltage, the output voltage varies by not more than 

+ 0.2%, The section containing gas-filled counters may 

be operated with any industrially available counter. The 
photoelectronic multiplier unit is designed for operation 

with the FEU-19M photomultiplier tube. 

PS-0.1 scaler for counting pulses exhibiting a 
random time distribution (Fig. 4). This instrument 
is one of the units in a general-purpose "line" of 
instruments designed for nuclear research, and having 
a resolving time of 0.1 psec. Channel capacity is 
10 pulses. Sensitivity is from 5 to 30 v (regulated in 
steps of 5 v each), and the scaling system used is a 
scale-of-ten, 

BDS-1 double-coincidence instrument (Fig. 5). 
This instrument, with two channels at its disposal, 
operates with two general-purpose scintillation trans- 
ducers and places a pulse across the output if the pulses 
received from the transducers fall within a specified 
time interval. The instrument is gated from negative 
pulses of 3 v amplitude and higher, and of 107° sec 
duration, The minimum resolving time is 2 sec, 
and the resolving time may be varied from 10 ® to 
10°" sec. One of the channels is provided with a 
variable delay line of 0-107’ sec. The instrument can 
handle a duty factor of 10° pulses per sec and may be 
operated as a time analyzer. 

UIS-1 spectrometric non-overloading pulse 
amplifier for amplifying pulses arriving from scintilla- 
tion transducers with NaI(T1) crystals. A preampli- 
fier (K= 10), mounted in the form of a separate unit 
housed inside the instrument, may be coupled to the 
main amplifier, if it is required to amplify pulses of 
~10mv. The amplifier consists of a power supply 
unit, a main amplifier, and the preamplifier. The 
maximum gain factor, with the preamplifier hooked 
in, is 10,000, and is 1000 without the amplifier. 

Pulse rise time is 0.2 usec, Peak linear output signal 
is 100 v. Duty factor is 10° pulses per sec. The 
integral nonlinearity is ~ 1%(for a capacitive load of 








Fig. 3. PS-5 scaler, 





Fig. 4. PS-0.1 scaler. 


Fig. 5. BDS-1 double-coincidence instrument. 


50°60 wmf). Gain instability is less than 1% for 8 hr 
of continuous operation, when the network voltage 
varies by + 10%, Input and output pulses are of 
positive and negative polarity. 

The RV timer is used for setting exact time 
intervals, and as a time-pulse oscillator. The timer 
is capable of setting 22 exposure times lasting anywhere 


from 10 sec to 12 hrs, with exposure-time error not 
in excess of 0.5%, Duration of time pulses 5-10 Lsecy 
Amplitude 25-50 v. 

SVV-9 miniaturized stabilized high-voltage recti- 
fier as power supplies for from one to three photomul- 
tiplier tubes. The output stabilized voltage is negative. 





Provisions are made for stepwise voltage variations 





from 0.8 to 2 kv, using a stepping switch (11-position 
switch) and slide-wire potentiometer, Load current 
may reach 2 mag 0. 04% instability in response to 
network voltage oscillations of + 10% or load current 
swings from 0 to 2ma, The ripple factor is 0. 004- 
0.01%, 

Photoelectronic multiplier tubes of various designs 
and types (spectrometric phototubes, timing phototubes) 
for working with directed light flux, and other photo- 
multiplier types: Photomultipliers with phtocathodes 
ranging up to 250 mm in diameter were put on display. 

Scintillation detectors using thallium-activated 
Nal, measuring up to 150 mm in diameter and up to 
60 mm in height, CsIl(T1) crystals 80 mm in diameter 
and 25 mm in height, anthracene crystals 40 mm in 
diameter and 20 mm in height, plastic scintillation 
phosphors based on polystyrene with PTP+POPOP addi- 
tive, measuring up to 200 mm in diameter and up to 
150 mm in height, etc. 

Instruments for testing and tuning electronic 
equipment and for taking their characteristics (oscilla- 
tors, amplifiers, quartz-crystal calibrators), vacuum- 
tube voltmeters, etc), were put on display, including 
the GKI-5 millimicrosecond-range square-wave 
generator, which yields pulses lasting 7-500 milli- 
microseconds (17 values) with amplitudes ranging from 
0 to 250 v and rising wavefront lasting 5-6 m psec; the 
GITA-3 pulse generator, which produces pulses ranging 
0 to 200 v in amplitude, with a pulse repetition rate 
of 50 cps, square-wave pulse duration of 2 sec, mini- 
mum rise time 107° sec, and pulse height stability of 
107° to 10, the GMI-1 millimicrosecond-range pulse 
generator, producing pulses of positive and negative 
polarity having amplitudes up to 45 v and durations 
ranging from 3 - 107% to 1.5 1078 sec, rising wave- 
front 1.5 + 107* sec, and pulse repetition rate of 
50-6000 cps; the 4E-2 electronic frequency meter 
designed for automatically measuring frequencies of 
electrical oscillations, pulse repetition rates in the 
10 cps to 10 Mc range of frequencies, measurements of 
pulse cycle duration (starting from 100 psec of higher), 
counting the number of electrical oscillations, mea- 
suring the duration of pulses and time intervals bet- 




















NEW ITEMS AT THE EXPOSITION IN 


ween pulses (longer than 10 psec), etc.; the OS-4 high- 
speed oscillograph, with a sweep duration of 10 to 
1000 m psec (in 5 discrete values) having a time- 
marker pulse lasting 2,10 and 50 m psec, with verti- 
cal-deflection amplifier bandwidth of 1000 Mc. 

Components and modular units used in the 
assembly and design of industrially available radio 
engineering and electronic equipment were displayed 
widely, these including: 10 standardized electronic 
circuit units, normalizers, signaling devices, constant 
voltage converters, trochotron starting stage, bipolar 
cathode-follower circuit, standard scaling circuit, elec- 
tronic stabilizers, gas stabilizers, semiconductor 
rectifiers, etc. ; 2) multipurpose components in the 
MTSsIS digital measuring system(100-256 character 
memiory units and matrices); 3) MnO,-Zn storage battery 
systems, including the GB-120, GB-300, GB-5400, 
and GB-900 batteries, which provide voltages up to 
1000 v and currents up to 100 Ha (service life extending 
to 250 hr); 4) silver-zinc storage battery system 
includings the STsS-15, STsS-12, STsS-3, STsS-1.5, 
STs-0.5, and STs-0.1 batteries, all of which function 
as new electrochemical reversible current sources (speci* 
fic capacity rating attains 100 w-hr per kg weight, peak 
capacity ranges from 0,12 to 16 amp-ht,with the weights 
of the batteries ranging respectively from 4.0 to 275.0 g); 
5) alkali manganese~zinc batteries including the MTs-1K, 
MTs-2K, MTs-3K, and MTs-4K battery types, serving 
as power supplies for semiconductor equipment and other 
instruments (the capacity ratings of these batteries, for 
a voltage of 0.8 v and discharge current of 5 ma, extends 
from 0.13 to 1.15 amp-hr; as to the size of the batteries: 
height ranges from 6.6 to 9.6 mm, diameters from 15,8 
to 30.3 mm, and weights from 4 to 20 g). 

TheFourth Scientific and Tecknical Conference 
on Nuclear Electronics will contribute to the successful 
development of work in the USSR on devising sophis- 
ticated electronic equipment in the service of nuclear 
physics and nuclear engineering. 

The proceedings of the conference will be published 
by Atomizdat in the form of a separate symposium in 
the book format. 





L. P. 


THE "ATOMIC 


ENERGY FOR PEACEFUL PURPOSES" PAVILION 


V. Orlov 


(From the exposition of the achievements of the national economy of the 


USSR) 


The " Atomic Energy for Peaceful Purposes" 
exhibition has been greatly expanded for 1959, and 
now occupies two full buildings. 
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The new building houses the following section of 
the exposition: “International Cooperation by the 
Soviet Union in the Field of the Peaceful Utilization 





Fig. 1. Scale model of 400 Mw(e) nuclear power station. 


of Atomic Energy, and the Joint Institute for Nuclear 
Studies," "Nuclear Power Engineering," and ™ Ther- 
monuclear Research, ” 

The "International Cooperation” section informs 
the visitor, by means of photographs, figures, and 
repeating newsreel shows, of the many aspects in which 
the Soviet Union is rendering technical aid and tender- 
ing scientific collaboration to other countries, in the 
training of specialists, in the building of nuclear 
research centers. Here, we see on display models of 
several reactors built by the Soviet Union in other 
countries, The activities of the Joint Institute for 
Nuclear Studies are illustrated ,in particular ,by the 
theme "Charged-Particle Accelerators," Iluminated 
moving-light display panels explain the operating 
principles of the synchrocyclotron and proton synchrot- 
ron, and a model of the 680 Mev synchrocyclotron is 
demonstrated. Some examples of nuclear research 
instrumentation are on display in the same hall. 

Models and live display panels are also featured 
in the demonstrations of nuclear-electric power gene- 


rating stations and particular nuclear reactors being 
planned, now under construction, or presently in 
operation; diagrams, a model,and the technological 
channel of the world's first nuclear-electric power 
station; a live panel and a scale model illustrating 

a 420 Mw atomic-electric power station powered by 
a water-cooled, water-moderated reactor; a live- 
panel and scale model of the 400 Mw uranium-gra- 
phite superheated-pr essurized-steam reactor power 
station (Fig. 1); an electrified moving display panel 
of the 50 Mw (e) boiling-water-reactor nuclear power 
plant now being built; a diagram of a test reactor 
burning gaseous fuel; a moving panel of a fast reactor, 
and photographs, drawings, and diagrams of fast 
reactors BR-1, BR-2, BR-3, BR-5, and BN-50, 

A scale model of the 150 Mw(e) nuclear power 
plant incorporating a gas-cooled heavy-water mode- 
rated reactor, located in the middle of the hall, was 
the center of much interest. The reactor will burn 
natural uranium fuel, The design of this power 
station was carried through by the joint efforts of 





Fig. 2, Scale model of the * Alpha" toroidal thermonuclear facility. 


Soviet and Czechoslovak specialists. A power station 
along these lines will be erected in Czechoslovakia. 

A new section, " Thermonuclear Research in the 
USSR,” has been added to the exhibits, 

The exhibits in this section explain fusion reac- 
tions of light nuclei, and tell of the conditions required 
for achieving controlled thermonuclear reactions, and 
of the history of the development of research on 
controlled fusion reactions in the USSR. 

The display illustrates three basic methods for 
obtaining a high-temperature plasma: pulsed pro- 
cesses, quasi-stabilized processes (charge stabiliza- 
tion by external magnetic field), and magnetic traps. 

An operating pulse-discharge experimental 
facility with a ceramic discharge chamber is demons- 
trated in this section. A facility such as this was the 
one in which neutrons emitted in the gaseous discharge 
process were first produced, at the Institute of Atomic 
Energy of the Academy of Sciences of the USSR, in 
1952. A porcelain discharge chamber for a pulsed 
thermonuclear facility is also demonstrated here. Also 
demonstrated are scale models of the first adiabatic 
magnetic trap (scale 1: 6), toroidal pinch chambers, 
arrangements in which a plasma ring is obtained in a 
transverse magnetic field, a magnetic adiabatic trap 
(scale 1; 10), the " Alpha” toroidal thermonuclear 
machine (Fig. 2). 
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The publishing and publicity activities relating to 
exhibitions on the peaceful uses of atomic energy are 
reflected by an exhibit set up inside the lightweight- 
glass structure of the new building. The stands display 
books published in the USSR and abroad, trophies, 
scrolls, medals awarded to exhibits of the USSR; colored 
posters advertise periodicals and books put out by 
Atomizdat press, and Soviet * Atoms for Peace" exhi- 
bits in foreign countries, Literature published by 
Atomizdat is also placed on sale in that section. 

Other sections of the pavilion have also been 
amply refurbished with new exhibits. 

The section "The Atomic Nucleus and Related 
Research Techniques” has on display a facility which 
provides a graphic demonstration, to the visitors, of 
the penetrating power of various types of radioactive 
radiations, and the shielding properties of many mate- 
rials, 

The section "Raw Materials for Nuclear Power, 
and Materials Used in Reactor Design” has fathered 
a new topical subsection; ™ Working of Uraniferous 
Ore Deposits." The high level of techniques achieved 
in safety practices in underground,highly productive, 
sublevel methods, sophisticated methods for grading 
and classifying ores (radiometric beneficiation), shaft- 
head facilities, techniques for hydraulic dust suppres- 
sion, etc. are shown. 





Also on display are radiometers used in ore grading: 
the RSR-3 "Strela" ore-grading radiometer, designed for 
assay of ores in large yield (in hoppers and automatic- 
loading and hauling devices) by means of y-ray emissions, 
and automatic classification of hoppers into groups depend- 
ing on the content of radioactive metal; the RSU-T radio- 
metric scintillation device, designed for express-analysis 
of ores by means of y radiation, while the ores are in 
motion along the conveyor belt. The efficiency of the 
scintillation detector is 5-8 pulses/ sec forl ur/hr. The 
resolving time is not longer than 15 psec. 

Several new logging and prospecting radiometers 
are also on display in this section: 

Portable y logging scintillation radiometer, PRKS, 
designed for measuring the intensity of radioactive 
emissions and evaluating the relative uranium and 
thorium content in ores directly at the ore deposits. 

The device may be operated under high-humidity 
conditions; 

KRS logging scintillation radiometer, capable of 
logging cores down to 1000 m below ground. Besides 
measuring the magnitude of the intensity of y radiations, 
the radiometer is capable of determining the relative 
uranium and thorium content in ores directly at the 
deposit site. The instrumental error over the 0-50°C 
temperature range is not in excess of 15%, The scin- 
tillation transducer is manufactured in hermetically 
sealed encapsulation; 

"Rupor™ ore-prospecting radiometer using gas- 
filled counters, available in two modifications; 

RPR-1 is designed for measuring the intensity of y 
emissions under field or mine conditions and for 
logging boreholes to depths of 2 m, and RPP-1 is de- 
signed for the detection of the presence of radioactive 
ores, and also for measurements of the intensity of y 
and 6 radiations from rocks under field conditions. 
The instrument is power-supplied by batteries. The 
weight of the entire portable assembly amounts to 

3 kg; 

"Kristall" SPR-2 prospecting scintillation radio- 
meter, for recording y radiations of energies above 
50 kev, and hard B emissions. Sensitivity threshold 
is not higher than 2,2 yr/hr for a gamma background 
level of 8 wr/hr. The instrument is capable of opera- 
tion over a range of temperatures from -20° to + 50°C 
and in 98% humidity. The entire weight of the 
assembly is 3 kg. 

" Sputnik-1" companion portable prospecting 
radiometer, for detecting radioactive substances, and 
for qualitative and quantitative analysis with respect 
to y and B emissions. The sensitivity threshold is not 
higher than 8 yr/ hr for a gamma background of 15 
Hr/hr, The instrument: is capable of operation over 
a range of temperatures from -10° to +40°C, and in 
95 4 3% humidity. Weight is 1.6 kg. 

“Pioneer"pocket y-and 8-indicator for detec- 
tion of radioactive emissions. 





























The section " Radiation Shielding and Health 
Physics” has been supplemented by a coverall of 
"lavsan" textile (distinguished by its heat stability, 
high mechanical strength, and stability to quite a 
few corrosive acids), the air-inflated overalls LG-4 
made from polyvinyl chloride by a high-frequency 
welding technique, and used in work under hot condi- 
tions where the surrounding air and equipment is hea- 
vily contaminated with radioactive matter. 

Some new instruments are on display in this 
section; 

SU-1 automatically actuated alarm annunciator for 
monitoring of radioactive contamination (Fig. 3). This 
instrument is designed to sound an automatic alarm when- 
ever the contamination on surfaces of the body or clothing 
by 8- and y-active substances exceeds safety levels. 
The sensitivity threshold is cohtrolled within the range 
0.03-3 pr/sec (with respect to y emissions from Co®), 
The alarm annunciator operates in the presence of a gamma 
background of up to 0.4 ur/sec, automatically compen- 
sating for the background; 





RK-01 pocket dosimeter for measuring y-and 6~- 
radiations in shops and laboratories, and under field 
conditions. The dosimeter circuitry is transistorized. 
The power is fed from a single 1 KSU-3 "Saturn" 
battery. It weighs 600 g. Error is within +10%, 

The dosimeter can operate without recharging for 
200 hr; 

KPN-2 portable fast and thermal neutron radio- 
meter, The instrument is equipped with three replace- 
able detectors, is capable of determining the intensity 
of fast neutron flux of energy ranging from 0.5 to 14 
Mev, and is also capable of measuring thermal neutron 








Fig. 3. Automatic alarm signaling device SU-1 
for monitoring radioactive contamination. 
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flux intensity. The instrumental error does not exceed 
+ 20%, of scale value, Power may be fed from the 
lines or from a STsS-5 battery; 


IFK-1 assembly for determining the total y-ex- 
posure dose obtained by personnel during the working 
day; 





DFE-10 photoelectric densitometer for measuring 
the diffuse optical density of segments of a developed 
photographic layer on a transparent substrate; 





IV-1 facility for measuring degree of contami- 
nation of water by 8-active substances. The instru- 
ment is provided with an alarm system indicating 
whenever water contamination exceeds a previously 
set level, It is capable of exercising continuous 








BRIEF COMMUNICATIONS 


USSR. In the Physics Institute of the Academy of Scien- 


ces of Armenian SSR an analyzer of pulse amplitudes 
has been developed, which employs logarithmic trans- 
formation for conversion of pulse amplitude into dura- 
tion; owing to this, the measuring of pulses in all chan- 
nels is realized with the same relative accuracy. 

The analyzer was designed to be used with a five- 
layer proportional counter in the field of cosmic-ray study. 


monitoring control over contamination of water by 
radioactive substances, Instrumental error stays with- 
in + 30% of nominal scale readings. 

The section "Isotope Production" is attracting 
increased attention because of the graphically illustra- 
ted “Table of Radioactive and Stable Isotopes Pro- 
duced in the USSR," Also on display are types of 
radioactive emission sources, internal packings, remote- 
control appurtenances for handling hot preparations, 
and a mock-up of a hot cell, 

Other sections in the exhibit, “ The Uses of Radio- 
active Isotopes in Industry," “The Uses of Radioac- 
tive Isotopes in Biology and Agriculture”, and " The 
Uses of Radi oactive Isotopes in Medicine " have also 
been supplemented with new displays. 


USSR. Tallinsk Experimental Plant of Instrumentation 
(Estonian SSR) has taken over manufacture of 12 types 
of instruments which are based on utilization of radio~ 
active isotopes; the instruments are intended for the 
automatization of various processes in the metallurgi- 
cal, chemical,mining, food and other branches of in- 
dustry and transportation. 

The first group of automatic regulators of level of 
petroleum products has been issued. 





NEW LITERATURE 


BOOKS , SYMPOSIA, PERIODICALS 





Verkhovskii, B, I, The Use of Radioactive 
Isotopes in Process Control, Academy of Sciences 
Press, 1959, 83 pages, 1 ruble, 25 kopeks. 

This book, in popular style, deals with the phe- 
nomenon of radioactive decay, the properties of radio- 
active radiation and methods for detecting radiations. 
Process control instruments based on the use of the 
penetrating power of radioactive radiation scattered 
by 68- and y-sources, or the ionizing power of radio- 
active radiation, and instruments utilizing neutron 
sources, are described. 

D' yachenko, P. E, Moscow, Mashgiz, 1958, 
215 pp. 8 rubles, 








D' yachenko, P. E. Use of Radioactive Isotopes 
in Industry. Moscow, Mashgiz, 1958, 215 pages, 8 
rubles, 

General information is provided on radioisotopes, 
methods for detecting them and methods used in intro- 
ducing them into the objects to be tested. Exam- 
ples of the use of radioactive isotopes in technical 
laboratories and in industry are cited, 

Weiss, K. F. Standard Radioactive Preparations. 
Translated from the German, Moscow, Fizmatgiz, 
1958. 244 pages, 8 rubles, 45 hopeks. 

This is a monograph treating of the properties 
and methods used in the manufacture and measurement 
of the absolute activity of standard radioactive prepa- 
rations: y- standards of radium, radiothorium, actin- 
ium, Co”, ce", Na”, Mn™, B -standards of Ra 
(D + E), (U + UX, + UX,) and Ti™. o standards of 
uranium, Pu”? Po??, a special section is reserved 
for the manufacture and calibration of neutron standards 

The book is written for the wide readership of 
persons working with hot preparations, 

Nukleonika (Poland). Vol. IV, No, 3 (1959).Janik, 
E. Scattering and Slowing Down of Neutrons by Bound 
Nuclei; Baran, W. Dependence of the Reactor Period 
on the Multiplication Constant Under Subcritical Con- 
ditions; Pinkas, K, Method for Studying Uranium Raw 
Materials of Domestic Origin and Results of the 
Method; Radoszewski, T. Absolute Methods for Mea- 
suring B and y activity; Lesynski, J. GM Halide 
Counters; Pawlikevicz, J. Nuclear Facilities in the 
Merchant Fuel; Colloquium Sponsored by the Radio- 
chemical Division of the Institute for Nuclear Research, 
Chronicle, 
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Journal of Nuclear Materials, April 1, 1959 
saw the appearance of the first issue of " Journal of 
Nuclear Materials” (English title), The periodical 
will be devoted to problems of metallurgy, ceramics, 
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